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by 
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ABSTRACT 

This  report  describes  the  methods  and  computer  programs  employed 
in  the  processing  of  ambient  noise  data  resulting  from  the 
deployment  of  ACOustic  DAta  Capsules  during  the  period  of  1971  to 
1973.  It  summarizes  the  techniques  and  programs  used  to  analyze 
third  octave  ambient  noise,  and  narrow  band  frequence  spectra, 
as  well  as  high  speed  analog-to-digital  processing. A  detailed 
discussion  of  the  ACODAC  system  can  be  obtained  from  the~W7H .0717 
report  submitted  by  Daubin,  Berteaux,  Bitterman,  et  al,  reference: 
W.1I.0, 1 . -72-87,  dated  November,  1972. 

Volume  II  is  a  description  of  the  hardware,  interconnections, 
system  schematics,  program  listings,  and  program  flow  diagrams. 
This  volume  1-  scheduled  to  be  in  print  by  mid-October  1973. 
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1.  GliNHRAL  j 

1 . 1  The  ACODAC  System 

The  ACODAC  system,  figure  1,  consists  of  six  hydrophones  in  a  vertical 
array  suspended  in  the  ocean.  The  instrument  pressure  vessel,  IPV,  the 
heart  of  the  system,  contains  the  necessary  electronics  to  actuate,  main¬ 
tain,  and  suspend  the  operation  of  the  six  hydrophones.  Acoustic  signals 
are  recorded  on  n  magnetic  tape  at  a  speed  of  15/160  inches  per  second.  Hach 
hydrophone  is  recorded  on  a  separate  channel  of  the  recorder  for  a  total  of 
,  six  channels  using  direct  recording  methods.  A  seventh  channel  records  an 
IRIG-C  time  code  generator  track  using  a  50  hertz  carrier.  The  time  code  is 
generated  continuously  and  records  days,  hours,  and  minutes  along  with  control 
bits  representing  the  gain  state  of  each  of  the  six  hydrophone  data  amplifiers 
and  system  identification.  Figure  2  shows  the  ACODAC  system  electronics  used. 

Since  the  hydrophone  array  can  be  deployed  in  the  ocean  for  a  period  of 
a  few  days  up  to  six  months,  the  recorder  electronics  has  the  capability  of 
cycling  on  and  off  at  a  rate  dictated  by  the  duty  cycle  selected  prior  to 
launch.  To  date,  the  100  per  cent  duty  cycle  has  been  used  for  most  of  the 
deployments.  This  gives  10.66  days  continuously.  During  the  recording  phase, 
control  signals  are  injected  on  the  six  data  tracks  of  the  recorder  indicat¬ 
ing  hydrophone  overload  conditions. 

A  200  hertz  control  signal  is  modulated  on  the  50  Hz  time  code  carrier 
to  indicate  the  start  of  the  calibration  every  six  hours.  The  calibration 
frequencies,  50  and  200  liz,  are  then  recorded  on  the  separate  hydrophone 
data  channels  through  each  of  the  six  hydrophone  data  amplifiers  for  a  five 
minute  period.  Any  duty  cycle  other  than  100  per  cent  is  detected  on  play¬ 
back  by  control  signals  modulated  on  the  time  code  carrier  indicating  start 
of  data  and  end  of  data. 

Once  the  system  is  recovered,  the  magnetic  tape  is  removed  and  played  back 
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Figure  1.  Acodac  Mooring  System 
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to  a  Hewlett  Packard  2100  digital  computer  through  a  General  Radio  Type  1921 
Real-Time  Analyzer. 

1.2  Computer  Analysis  Programs 

Four  general  methods  are  employed  for  the  analysis  of  ambient  noise 
data  from  ACODAC  deployments  using  digital  computer  techriques: 

(1)  A  third  octave  analysis  using  a  Gene/al  Radio  analyzer  with  a 
Hewlett  Packard  2100  computer  to  obtain  average  spectrum  levels 
reference  a  microbar  at  the  desired  third  octave  bands. 

(2)  A  narrow  band  analysis  using  the  Sigma  7  anelog-to-digital 
convertc-  to  obtain  narrow  uand  frequency  spectra  of  a 
continuous  time  series. 

(3)  A  continuous  wave  analysis  system  using  one  hertz  bandwid*1' 
filters  to  determine  transmission  loss- 

(4)  A  high  speed  digitizing  process  on  the  Hewlett  Packard  2100 
computer  to  analyze  explosive  shots  to  obtain  transmission 
loss  data. 

Computer  programs  and  techniques  were  developed  to  reduce  the  raw  ACODAC 
ambient  data  to  meaningful  graphic  plots  and  statistical  information  which 
were  representative  of  the  ambient  noise  conditions  during  the  deployment 
phase.  Hie  type  of  output  data  was  carefuily  chosen  after  considerable 
s  tudy  to  provide  only  the  necessary  information  for  objective  acoustic 
analysis . 

It  should  be  emphasized  that  initial  attempts  of  automatic  data  process¬ 
ing  yielded  large  volumes  of  data  resulting  both  from  the  editing  phases 
and  from  the  final  stages  of  processing  which  were  thought  necessary'  in  the 
analysis  of  the  system.  From  this  experience,  the  interim  data  processes 
were  drastically  reduced  to  provide  only  that  information  which  was  thought 
absolutely  necessary  to  ascertain  the  quality  and  accuracy  of  data. 
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This  report  establishes  the  techniques  employed  in  converting  tne  raw 
ambient  noise  data  into  formats  that  can  be  used  with  the  appropriate 
statistical  subroutines  to  obtain  the  desired  acoustics  analysis. 

1 . 3  Summary  of  Output  Data 

The  following  represents  a  sunmary  of  the  output  data  available  from 
the  various  data  processing  systems  followed  by  a  brief  discussion  of  the 
methods  used  to  produce  each  of  the  data  outputs.  It  is  intended  that  the 
data  described  below  constitutes  a  ’’standard  data  package"  available  with 
the  processing  of  each  ACODAC  deployment  at  the  Woods  Hcie  Occ-anograpfic 
Institution. 

It  is  realized  that  there  can  be  varied  presentations  of  acoustic 
data  in  the  form  of  graphical  displays  and  tabular  statistical  summaries 
depending  upon  the  methods  of  analysis  and  results  desired.  However,  the 
formats  presented  have  been  chosen  to  aid  in  the  general  acoustic  analysis 
of  the  ACODAC  deployments  because  of  the  large  volumes  of  data  available. 

1.3.1  Third  Octave  Analysis  of  Ambient  Noise  Data  Using  a  Hewlett 

Packard  2100  Computer  with  a  General  Radio  RMS  Detector 

1.  Digital  magnetic  tape  of  ambient  noise  data  -  10  second 
average  spectrum  levels  corrected  to  microbar  levels. 

2.  Calibration  summary  tape. 

3.  Brush  analog  plots  of  ambient  noise  -  10  second  average 
spectrum  levels  corrected  to  microbar  levels. 

4.  Lrror  diagnostic  message  printout. 

1.3.2  CW  Analysis 

1.  Digital  magnetic  tape  of  ambient  noise  data  -  10  second 
averages  corrected  to  microbar  levels. 
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2.  Brush  analog  plots  of  narrow  band  levels. 

3.  Lrror  diagnostic  message  printout. 

1.3.3  bditing  Process  for  Third  0ct3ve  Digital  Magnetic  Tapes 

1.  Produces  edited  or  "clean"  tape  of  the  digital  input  ambient 
noise  data  from  the  third  octave  General  Radio  system.  All 
shots  and  spurious  signals  are  eliminated  during  the  editing 
process.  During  the  editing  process,  the  10  second  averages 
can  be  corrected  for  the  following: 

a.  change  relative  time  to  Greenwich  Mean  Time 

b.  insert  a  band  level  correction  for  all  10  second 
averages,  or  for  a  selected  time  window. 

c.  insert  a  level  correction  for  all  bands  for  all 

10  second  averages,  or  for  a  selected  time  window. 

d.  correct  the  10  second  averages  for  a  tiansfer  func¬ 
tion  shift  for  the  selected  six  hour  period  desired 
using  the  transfer  function  computed  and  stored  in 
the  data  record  for  that  period. 

e.  select  option  to  include  all  shot  intervals  on  the 
"clean"  tape. 

During  the  editing  process  (producing  a  "clean"  tape),  ary 
combination  of  15  averages  can  be  obtained  with  the  high  speed 
printer  for  the  entire  digital  tape.  The  averaging  base  is 
selectable  a:  run  time  and  the  editing  program  has  been  made 
flexible  to  include  any  averaging  period  desired.  Also,  during 
the  editing  process,  a  summary  is  printed  tagging  any  signal 
level  that  exceeds  the  ambient  noise  window  such  as  shots  or 
gain  changes  not  corrected  for  by  the  digitizing  program.  These 
diagnostic  messages  are  later  cross  checked  with  the  Brush  plots 
of  the  same  data  to  determine  the  nature  of  the  excessive  levels. 


-7- 


1.3.4  Processing  Ambient  Noise  Data  from  Clean  Tapes 

Once  the  editing  process  lias  been  completed,  several  options 
are  available  to  the  user  as  final  data  formats: 

a.  The  user  can  select  to  produce  two  and  ten  minute 
averages  of  the  average  spectrum  levels  for  the  en¬ 
tire  edited  tape  or  for  selected  time  windows.  The 
data  averaging  base  can  be  selected  it  process  time  and 
any  averaging  period  can  be  computed,  i.e.,  one  and  six 
hour  averages,  one  minute,  ten  minute,  and  sixty  min¬ 
ute  time  averages,  etc. 

b.  The  user  can  select  to  obtain  1  2  hour  statistical 
information  such  as  mean,  standard  deviation,  kurtosis, 
and  frequency  histograms.  These  statistics  are  avail¬ 
able  for  all  frequency  bands  of  interest.  See  fig¬ 
ure  71  for  sample  of  printout. 

c.  The  user  can  select  to  produce  a  digital  plot  of  the 
average  spectrum  levels  for  the  entire  period  of  the 
deployment  o-  for  a  selected  time  window.  Data  are 
plotted  using  a  one  minute  mean  of  tie  10  second 
averages.  Normally,  three  to  four  bands  are  plotted 
as  a  function  of  time;  however,  the  program  can  plot 
the  entire  15  bands.  Figure  23  shows  a  sample  plot. 

< .  7.5  Transmission  Loss  Data  From  CW  Sources 

A  transmission  loss  table  can  be  obtained  using  the  edited 
copy  of  the  narrow  band  system  data  tape  as  inr’it.  The  program 
computes  the  transmission  loss  of  the  narrow  hand  levels  of  a 
particular  source  along  with  the  ambient  noise  based  upon  a  two 
minute  mean.  The  table  provides  tiie  information  to  correlate 
transmission  lo«s  from  a  particular  CK  source  with  a  ship's  track. 


-8- 


1.3.6  Narrow  Band  Analysis  (Sigma  7  Computer) 

The  narrow  band  analysis  method  requires  the  following 
prog  rams : 

a.  An  analog-to-digital  conversion  program  to  input 
ambient  noise  data  from  the  Sangamo  recorder  via  the 
A/D  converter  to  the  Sigma  7  computer.  All  six  chan¬ 
nels  are  digitized  sequentially  at  a  54  KHz  sampling 

i ate,  transmitted  to  the  computer,  and  then  dumped 
onto  magnetic  tape. 

b.  A  media  conversion  program  which  inputs  digital 
tape  produced  in  a.  (above)  along  with  a  gain  tape 
produced  from  the  th.  rd  octave  digitizing  system, 
mei^es  and  re- formats  the  data  to  make  it  compatible 
with  existing  time  series  analysis  programs. 

c.  The  data  produced  from  the  media  conversion  is  then 
processed  through  a  time  series  analysis  program  to 
produce  standaid  plots  of  frequency  distributions, 
cross,  and  auto  correlation  studies. 


2.  DATA  ACQUISITION  SYSTTM 

The  ACODAC  system  consists  of  six  hydrophones  suspended  in  a  vertical  array 
with  a  data  capsule  (Instrument  Pressure  Vessel)  containing  all  the  required 
electronics  to  record  acoustic  signals  onto  magnetic  tape  from  the  hydrophone 
array.  Referring  to  figure  1,  the  electronics  contain  data  amplifiers,  over¬ 
load  detecting  circuits,  calibration  circuits,  and  a  master  clock  to  control 
the  internal  timing  of  the  system. 

The  automatic  gain  control  circuits  are  capable  of  injecting  four  levels 
of  gain  to  permit  recording  within  the  dynamic  range  between  the  hydrophone 
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and  the  instrument  recorder.  All  internal  timing  is  derived  from  the  master 
clock  which  generates  IfHG-C  time  code  information  with  a  50  hertz  carrier  for 
continuous  recording  of  days,  hours,  and  minutes  on  to  t he  time  code  track  of  the 
recorder. 

Gain  state  information  is  injects  on  the  time  code  via  the  control 
function  bits  in  the  IRIG-C  minute  frame.  Two  gain  bits  per  hydrophone  plus 
two  additional  redundant  gain  bits  are  injected.  The  gain  level  injection  occurs 
once  per  minute.  The  gain  level  of  each  data  amplifier  is  determined  by  a 
circuit  that  integrates  the  incoming  signal  for  one  minute.  A  decision  is  made 
to  determine  which  of  the  four  gain  levels  the  amplifier  gain  should  be  set  to. 

The  gain  state  is  then  switched  at  the  end  of  the  minute  and  the  gain  bits  in¬ 
jected  on  the  time  code  carrier. 

The  incoming  signal  once  amplified  is  recorded  on  a  magnetic  tape  re¬ 
corder  at  a  speed  of  15/160th  inches  per  second  using  the  direct  recording  method. 
Seven  tracks  are  recorded:  six  data  tracks  and  one  time  coJe  track.  The  elec¬ 
tronics  has  the  capability  of  discontinuous  recording  selectable  prior  to  launch, 
nils  method  can  be  used  to  increase  the  deployment  time.  Normally,  a  100  per  cent 
duty  cycle  is  employed,  but  a  5  per  cent  duty  cycle  has  been  usea  (1  minute 
on  and  19  minutes  off) . 

A  calibration  sequence  is  also  employed  in  the  IPV  electronics  to  allow 
injections  of  50  and  200  hertz  calibration  signals  once  every  six  hours.  Tnese 
calibration  sequences  are  used  later  in  the  data  processing  system  to  determine 
dynamic  changes  in  the  instrumentation  and  recording  media.  Table  1  shows  a  list 
of  the  control  signals  being  used. 
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TABLE  1 

ACODAC  CONTROL  SIGNALS 

INPUT  TO  COMPUTER 

SIGNAL  HOW  RECORDED  AND  ACTION  REQUIRED 

i 

Start-of-Data  I  50  Hz  +  150  Hz  for  4  sec.  Posit' ve  pulse  for  .05±  .01 

,  1  seconds  -  Integrations  shall 

begin  after  delay  of  .013 
seconds 

End-of-Data  75  Hz  150  llz  for  4  sec.  Positive  pulse  for  .05±  .01 

‘  seconds  -  stop  integrations; 

i  if  pulse  comes  during  integra- 

i  tion  period  data  shall  be  re¬ 

jected  and  noted  on  readout 
tape. 

Error  75  Hz  +  200  Hz  for  100  Positive  pulse  for  duration 

milliseconds  (real  time)  of  error  signal  (1.25  ms 
minimum.  min.)  -  Data  during  integra¬ 

tion  period  shall  be  rejected 
and  noted  on  readout  tape. 

Calibration  Signal  50  Hz  for  1  min.  at  4  Positive  pulse  for  duration 

amplitudes,  700  llz  for  of  calibration  signal  - 
1  min.  at  4  amplitudes  Process  normally  but  note  as 

calibration  signal  on  read¬ 
out  tape. 

Calibration  Sequence  200  Hz  modulated  on  Indicates  calibration 

Indication  time  code  for  5  min.  sequence. 

Time  Code  Continuously  on  Parallel  BCD  (21  bits)  in 

track  7  for  first  30s  days,  hours,  minutes  ••  Re- 

of  each  min  cord  on  readout  tape  once 

per  minute. 

Data  Amplifier  Settings  Included  in  time/state  Two  Bits  -  Record  on  readout 

code.  Follows  time  code  tape  once  per  minute, 
on  edge  track.  Occupies  00  corresponds  to  10  db  gain 

part  of  second  30s  of  01  corresponds  to  20  db  gain 

each  min.  10  corresponds  to  30  db  gain 

11  corresponds  to  40  db  gain 
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ANALVSIS  01-  AMU  HINT  N01SL  DATA  USING  THIRD  OCTAVli  FILTliRS  AND  10  SLCOND 
RMS  AVtiKAGhS 

3 . 1  The  Third  Octave  Analysis  System 

Figure  3  shows  the  instrument  connections  for  the  third  octave  analysis 
system.  The  selected  hydrophone  signal  from  the  Sangamo  direct  recording 
playback  unit  is  passed  through  the  control  function  indicator  (CFI)  panel 
where  calibration  and  overload  signals  are  detected  through  band  pass  fil¬ 
ters.  See  table  1  for  detection  frequencies.  From  the  CFI  panel,  the  wide 
band  signal  enters  the  General  Radio  1926  third  octave  Multifilter  through 
bands  31  to  44.  These  bands  correspond  to  the  frequencies  of  interest  of 
IS  to  300  hertz  at  an  80:1  playback  rate.  The  signals  are  then  fed  into 
the  General  Radio  1921  Analyzer  for  digitizirg.  The  wide  band  signal  is 
also  connected  to  band  45  but  bypasses  the  miltifilter  before  it  enters  the 
GR  analyzer.  The  output  accumulators  of  these  15  bands  are  then  connected 
directly  to  the  Hewlett  Packard  2100  computer  through  a  16  bit  microcir¬ 
cuit  interface  card. 

3.1.1  The  General  Radio  .Analyzer 

Lach  h.ydrophone  channel  of  the  analog  tape  is  played 
separately  through  the  General  Radio  1921  analyzer  to  con¬ 
vert  the  data  from  analog  to  digital.  In  the  analyzer,  or 
more  precisely,  the  rms  detector,  the  data  once  bandpassed 
through  the  series  of  third  octave  filters  are  -ampled  and  fed 
through  digital  circuits  that  are  time  shared  cn  all  channels. 
The  rms  detector  simultaneously  computes  the  rms  levels  for 
each  filter  band  taking  a  maximum  of  1024  samples  per  band 
and  converting  each  sample  to  a  digital  binary  number  and  squar¬ 
ing  it.  The  squared  values  are  accumulated  in  memory  registers 
until  the  end  of  the  selectable  integration  period.  At  the  end 


-13- 


of  the  integration  period  the  sum-of -squares  value  is  converted 
to  decibels  for  output  presentation  to  the  computer.  The  analyzer 
performs  the  above  operations  by  solving  the  equation: 

N  V 

rms  =  10  Log  £  [-^i]2}db  (3.1-1) 

i=  1  o 

where  Vq  is  the  voltage  at  the  reference  level,  and  N  is  the 
number  of  samples. 

At  this  time,  the  computer  is  ready  to  strobe  in  the  rms 
value  of  each  band.  See  Appendix  E  for  more  detailed  discus¬ 
sion  of  the  General  Radio. 

3.1.2  Time  Code  Reader 

The  ACODAC  processing  system  employs  a  Datametrics  model  344 
time  code  reader  that  converts  IRIG-C  serial  time  code  with  a 
modulated  50  hertz  carrier  to  parallel  BCD  days,  hours,  and 
minutes  t:me-of-year  information.  The  BCD  data  is  connected 
through  a  back  panel  connecter  to  the  computer  through  a  32  bit 
TTL  interface  card.  A  read  strobe  pulse  is  transmitted  to  the 
computer  on;e  per  minute  from  the  time  code  reader  so  that  the 
computer  knows  when  to  read  valid  BCD  information.  The  strobe 
occurs  at  the  on-time  reference  marker  of  the  minute  frame. 

Twenty -three  binary  digits  represent  days,  hours,  and  minutes. 

See  figure  4  for  IRIG-C  format. 

The  Datametrics  time  code  reader  accepts  IKIG-C  serial  code 
from  track  7  of  the  ACODAC  ambient  noise  analog  tape.  Time  of 
the  year  information  is  transmitted  continuously  in  patterns  of 
ones  and  zeros  using  amplitude  modulation.  The  time  code  reader 
decodes  the  bit  patterns  by  using  envelc.^c  detection,  a  one  bit 
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is  represented  by  a  series  of  five  cycles  of  the  carrier  having 
a  ratio  of  3:1  amplitude  from  a  zero  bit  pattern.  Three  of  the 
five  cycles  of  the  bit  pattern  must  be  detected  in  the  envelope 
with  an  amplitude  ratio  of  3: l  with  respect  to  a  zero  bit  pat¬ 
tern  in  order  to  decode  the  envelope  as  a  "one"  bit. 

The  BCD  time  of  year  code  can  be  transmitted  in  the  first 
21  seconds  of  the  minute  frame.  The  remainder  of  the  minute 
time  frame  from  25  seconds  to  60  seconds  is  used  to  transmit 
control  functions,  and  in  the  case  of  ACODAC  is  used  to  transmit 
the  gain  states  of  the  six  hydrophone  data  amplifiers. 

Computer  program  timing  is  obtained  by  a  one  pulse  per 
minute  signal  coming  from  the  time  code  reader.  T!ii3  signal  is 
derived  from  the  50  hertz  carrier  being  read  from  the  IRIG-C 
time  code  track.  Selectable  playback  rates  are  provided  by  a 
rotary  switch  on  the  front  panel.  Rates  of  80:1,  8:1,  10:1, 
and  1:1  are  selectable  although  present  methods  employ  an  80:1 
playback  rate  (the  Datametrics  reader  also  accepts  100  Hz 
carriers) . 

The  hydrophone  amplifier  gain  states  transmitted  in  the 
control  function  part  of  the  time  code  minute  frame  are 
detected  in  the  time  code  reader  and  output  through  a  rear  panel 
connector  directly  to  the  gain  state  decoder  panel. 

An  IRIG-D  slow  code  (one  hour  code  frame)  is  generated  in 
the  time  code  reader  and  available  as  an  output  through  a  8NC 
connector  on  the  rear  panel.  The  one-hour  time  frame  is  used 
as  a  timing  reference  marker  on  the  Brush  analog  recorder. 

Sec  figure  5  for  t Ire  IRIG-D  format. 
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3.1.3  Gain  State  Decoder 

As  was  stated  above,  the  gain  states  of  the  six  hydrophone 
data  amplifiers  are  transmitted  on  the  IR1G-C  carrier  in  the  con¬ 
trol  function  position,  Each  hydrophone  data  amplifier  gain  state 
is  represented  by  two  bits  indicating  the  gain  setting  of  a 
particular  hydrophone  during  that  minute  interval.  Since  the  IPV 
electronics  can  only  change  the  data  amplifier  gain  settings  once 
during  the  minute,  the  amplifier  gain  remains  constant  during  the 
entire  minute.  The  change  occurs  at  the  zero  second  reference  mark 
of  the  minute  frame.  Since  the  decoding  of  gain  bits  is  so  im¬ 
portant,  the  system  has  been  designed  to  transmit  gain  by  using 
two  sets  of  two  bits  for  each  hydrophone.  A  total  of  24  bits 
is  recorded  on  the  time  code  carrier  for  six  hydrophones.  Each 
hydrophone  has  a  range  of  40  db  of  programmable  gain  range  ar¬ 
ranged  as  follows: 

Bits  Gain  (db) 


00 

10 

01 

20 

10 

30 

11 

40 

The  gain  state  decoder  performs  three  functions:  1)  checks 
the  redundant  gain  bits  from  all  hydrophone  channels,  2)  displays 
the  gain  bits  and  the  gain  error  bit  from  one  hydrophone  channel, 
and  3)  reads  the  time  code  during  playback  into  the  computer. 

The  gain  bits  are  read  from  the  state  code  bits  according  to 
the  pin  assignments  of  table  2  and  table  3. 

1.  The  gain  error  bit  (E)  for  any  channel  is  "1"  when  A  i  AR  or 
B  i  BR  for  that  channel.  This  can  be  expressed  as 
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TABLE  2 

TIME  CODE  READER 


OUTPUT  PIN  ASSIGNMENTS  -  GAIN  INFORMATION 


(Connector  J4) 


Pin  #  Function 


Pin  #  Function 


1 

Ch 

1 

- 

A1 

26 

Ch 

1 

-  AIR 

2 

Ch 

1 

- 

B1 

27 

Ch 

1 

-  B1R 

3 

Ch 

2 

- 

A2 

28 

Ch 

2 

-  A2R 

4 

Ch 

2 

~ 

B2 

29 

Ch 

2 

-  B2R 

5 

Ch 

3 

- 

A3 

30 

Ch 

3 

-  A3R 

6 

Ch 

3 

- 

B3 

31 

Ch 

3 

-  B3R 

7 

Ch 

4 

- 

A4 

32 

Ch 

4 

-  A4R 

8 

Ch 

4 

- 

B4 

33 

Ch 

4 

-  B4R 

1) 

Ch 

5 

- 

AS 

34 

Ch 

5 

-  ASR 

10 

Ch 

5 

• 

B5 

35 

Ch 

5 

-  B5R 

11 

Ch 

6 

- 

A6 

36 

Ch 

6 

-  A6R 

12 

Ch 

6 

- 

B6 

37 

Ch 

6 

-  B6R 

TABLE  5 


Gain 

GAIN  BIT  ASSIGNMENTS 

+10  db  +20  db 

+  30  db 

+  40  db 

A 

0 

1 

0 

1 

B 

0 

0 

1 

1 

Computer  bit 

IS 

14 

13 

12 

»  11 

10 

9  8 

7 

6 

5 

4 

3 

2  1 

0 

First  word 

0 

E4 

B4 

A4 

0 

E3 

B3  A3 

0 

E2 

B2 

A2 

0 

El  B1 

A1 

DSI  pin  assignment 

10 

Z 

X 

R 

N 

F 

D  2 

8 

6 

V 

T 

L 

J  B 

4 

Second  word 

0 

0 

0 

0 

0 

0 

0  0 

0 

E6 

B6 

A6 

0 

E5  B5 

AS 

DSI  pin  assignment 

9 

7 

W 

U 

M 

K  C 

5 
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E  =  Gain  error  bit  =  A*AR  +  A*AR  +  B*BR  ♦  B*BR 
The  output  to  the  computer  for  each  channel  is  A,  B  and  E. 

The  18  bits  are  assigned  to  the  DSI  inputs  as  shown  in  table  3. 

2.  A  front  panel  switch  selects  AR,  BR  and  E  from  the  indicated 
channel  to  drive  three  panel  lights  (light  on  =  "1"). 

3.  Every  time  the  time  code  data  goes  from  59  minutes  to  00  min¬ 
utes  a  one-second- long  tick  is  generated  on  the  Brush  recorder 
event  marker. 

3.1.4  Control  Function  Indicator 

The  Control  Function  Indicator  Pane!  (CFI)  is  used  to  detect 
control  frequencies  injected  on  either  the  hydrophone  data  track 
or  the  IRIG-C  time  code  track  of  the  analog  tape  recorder  during 
deployment.  The  CFI  panel  contains  band-pass  filters  and  associated 
digital  logic  to  detect  these  control  frequencies  injected.  Any 
combination  of  four  bits  can  be  set  to  inform  the  computer  that  a 
control  frequency  has  been  injected.  These  art: 

Bit  0  -  Start  of  Data  (Duty  cycling  only) 

Bit  1  -  End  of  Data  (Duty  cycling  only) 

Bit  2  -  Error  -  overload  indication 

Bit  3  -  Calibration  sequence  (injected  once  per  six  hours) 
Figure  6  shows  a  block  diagram  of  the  CFI  panel.  A  control  in¬ 
terrupt  informs  the  computer  when  any  bit  has  been  set  by  the  CFI 
panel.  The  computer  then  strobes  all  four  bits  into  the  computer 
from  the  CFI  output  register,  stores  the  bit  configuration  in  the 
computer  program,  and  then  resets  the  output  register  so  that  a 
bit  set  can  only  be  read  once. 

During  the  deployment  phase,  the  error  or  overload  signal 
will  be  injected  through  the  hydrophone  data  amplifier  and 


20 


Figure  6.  Data  Reduction  Block  Diagram 
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recorded  on  the  analog  tape  when  the  input  hydrophone  signal  over¬ 
loads  the  data  amplifier.  The  remaining  three  signals,  end-of- 
data,  start -of-data,  and  calibration  are  injected  on  the  time  code 
carrier  since  they  are  common  to  all  six  hydrophones, 

3.1.5  Ihe  IIP  2100  Digital  Computer 

The  Hewlett  Packard  2100  digital  computer  is  a  standard 
general  purpose  digital  computer  used  to  control  and  process  the 
ambient  noise  data  recorded  on  analog  tape.  The  computer  has  16,000 
words  of  memory  and  employs  a  priority  interrupt  system  to  service 
external  devices.  Two  direct  memory  access  channels  are  available 
for  transfer  of  high  speed  data  to  and  from  the  computer.  One 
channel  is  used  to  strobe  in  data  from  the  General  Radio  analyzer. 
The  other  channel  is  used  to  output  ambient  noise  data  records  to 
digital  magnetic  tape.  Some  of  the  standard  peripheral  equipment 
attached  to  the  computer  are: 

Model  IIP  2748A  High  speed  reader 

Model  HP12925A  High  speed  paper  tape  punch 

Model  HP  7970B  9-track  digital  magnetic  tape  (800  BPI,  45  ips) 
Tektronix  4002  Input-output  CRT  display  and  graphics  console 
Model  365  Calcomp  Plotter  (digital) 

Model  ASR-35  Input-output  console  (backup) 

Model  486  Brush  analog  recorder 

3.1.6  The  Brush  Analog  Plotter 

The  eight  channel  Brush  analog  plotter,  Gould  Instruments, 
model  486,  is  connected  to  the  computer  through  digital-to-analog 
output  interface  cards.  Data  points  representing  10  second  averages 
of  corrected  ambient  noise  data  (re:  db /microbar)  for  selected 
third  octave  frequency  bands  are  converted  to  analog  signals  and 
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plotted  on  the  Brush  recorder.  The  data  points  are  plotted  con¬ 
tinuously  through  the  entire  hydrophone  process  period.  The 
Brush  plotter  has  eight  amplifiers,  model  134214-00,  capable 
of  operating  in  the  log  or  linear  mode.  The  present  configura¬ 
tion  requires  the  amplifiers  to  operate  in  the  linear  mode  with 
an  input  range  of  10  volts  full  scale  representing  an  output 
dynamic  range  of  SO  db. 

3.1.7  Sangamo  Playback  Analog  Tape  Recorder 

The  Sangamo  analog  recorder  playback  tape  unit,  model  3561, 
is  used  to  reproduce  the  acoustic  signals  recorded  on  the  Geotech 
tape  recorder.  The  Sangamo  tape  unit  is  equipped  with  7  data 
tracks  and  7  direct  reproduce  amplifiers  to  play  back  6  hydro¬ 
phone  data  channels  and  IRIG-C  time  code.  Hie  tape  unit  has  a 
selectable  speed  dial  capable  of  recording  and  playback  of  analog 
tapes  of  speeds  from  1  7/8  ips  .o  30  ips. 

3.2  Input  Data 

This  method  involves  the  analysis  of  a  single  hydrophone  data  channel  at 
an  80:1  playback  rate  (recorded  at  15/160  ips  and  played  back  at  7.5  ips). 

A  hydrophone  channel  is  played  through  a  General  Radio  1921  Analyzer  System  to 
obtain  10  second  averages  (1/8  second  in  playback)  every  12  seconds.  An 
IRIG-C  time  code  track,  recorded  during  deployment,  is  channeled  through  a 
gain  state  decoder  to  obtain  the  data  amplifier  gain  settings  once  per  minute. 
The  following  information  is  obtained  from  the  data  channels  recorded  during 
the  deployment  phase. 

1.  Wide  band  ambient  noise  data  from  15  to  300  Hz. 

2.  Overload  injection  from  data  capsule  (75  and  200  Ilz  signals)  indicating 
signal  saturation  above  minimum  gain  of  10  db. 

3. 3  Third  Octave  Data  Output 

The  following  forms  of  output  data  are  obtained  from  the  third  octave 
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analysis  process: 

(1)  Digital  magnetic  tape  record.  Each  record  contains  seven  minutes  of  data 
or  35  integrations  for  each  of  15  bands  (15-300  Hz).  The  sample  base  is  10 
second  averages.  Each  data  minute  contains  a  five  word  header  providing  the 
day  number,  hour  and  minute,  gain  setting,  data  type,  and  transfer  function. 

(2)  Analog  output  plot.  An  analog  plot  containing  8  channels  of  data  con¬ 
tinuously  recorded  for  the  entire  processing  period  of  the  hydrophone  through 
the  third  octave  system.  The  analog  channels  record  the  following  information: 

Channels  1-4:  any  four  out  of  15  bands  selected  for  plotting. 

The  average  spectrum  levels,  in  10  second  averages,  are  in 
microbars.  Generally,  the  selected  third  octave  bands  are: 

Channel  1:  SO  Hz 
Channel  2:  100  Hz 
Channel  3:  75  Hz 

Channel  4:  200  Hz 

Channel  5:  is  the  computer  output  gain  state  level  of  the  hydrophone 
being  processed.  Range  in  dc  level  represents  data  amplifier 
gains  of  10,  20,  30,  or  40  db. 

Channel  6:  Hydrophone  gain  as  read  directly  from  the  time  code  reader. 

This  redundancy  has  beer  useful  to  detect  incorrect  gain  state 
decoding  via  the  computer  processing.  If  channels  5  and  6  do  not 
agree,  problems  exist  either  in  the  computer  processing  program 
or  system  hardware. 

Channel  7:  Uncorrected  wide  band  signal  from  processing  hydrophone  through  the 
50  Hz  third  octave  filter  plotted  directly  on  the  Brush  plotter 
(not  through  computer).  Channel  7  is  used  to  check  validity  of 
data  being  processed  through  the  computer.  Visual  checks  are  made 
between  the  uncorrected  50  Hz  signal  level  of  channel  7  with  the 
corrected  50  Hz  signal  level  being  recorded  on  channel  1  from  the 
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computer.  In  crder  to  compare  the  levels  of  channel  7  with 
channel  1,  a  process  can  be  used  to  convert  the  un corrected  50  Hz 
signals  of  channel  7  and  then  to  compare  these  results  with  the 
signal  level  of  channel  1.  If  discrepancies  occur,  then  the  third 
octave  process  should  be  discontinued  until  the  problem  is  resolved. 
Incorrect  parameters  inserted  in  the  computer  program  during  initial¬ 
ization  are  indicative  of  discrepancies  between  these  two  channels. 

Channel  8:  This  channel  is  used  to  provide  a  suitable  time  base  when 

analyzing  the  Brush  plots.  An  IRIG-D  time  code  is  used  from  the 
time  code  reader  which  represents  a  one-hour  repetition  rate. 

Bach  DC  level  excursion  represents  one  minute  of  hydrophone  time, 
and  the  code  frame  repeats  every  hour. 

(3)  Post-Calibration  summary  tape.  This  calibration  summary  tape  is  used  to 
analyze  the  data  computed  during  the  calibration  cycle  once  every  six  hours. 

The  post  calibration  analysis  tape  provides  useful  information  to  determine 
the  quality  of  calibration  data  obtained  during  the  five  minute  calibration 
sequence  every  six  hours  of  deployment  time.  At  the  completion  of  processing 
a  hydrophone,  the  summary  tape  is  read  into  the  computer  and  then  printed 

out  on  the  high  speed  printer.  Figures  10  and  11  show  a  typical  output  from  this 
process . 

(4)  Hydrophone  gain  tape.  This  digital  output  tape  is  optional.  If  selected 
for  output,  the  program  writes  the  hydrophone  gain  states  for  all  the  six 
hydrophones  in  the  deployment  onto  magnetic  tape  once  per  hydrophone  minute. 

In  addition,  day  number,  hour  and  minute,  and  transfer  function  are  •ritten 
in  the  record.  This  output  tape  is  used  in  the  narrow  band  analysis  to  pro¬ 
vide  hydrophone  gain  information  along  with  digitized  ambient  noise  data  to  a 
fast  Fourier  transform  program. 

(5)  Summary  error  message  pr.ntout.  All  errors  detected  or  abnormal  indica¬ 


tions  encountered  during  the  processing  of  the  hydrophone  through  the  third 
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octave  system  are  recorded  on  hardcopy.  These  printouts  provide  information 
with  respect  to  the  state  of  processing  at  any  given  time.  For  example,  over¬ 
load  indications  are  recorded  as  messages  once  per  minute  of  clock  time  (80 
minutes  of  hydrophone  time)  to  inform  the  operator  as  to  the  state  of  process¬ 
ing.  It  also  proivdes  a  hard  copy  of  all  the  parameters  inserted  in  the  com¬ 
puter  program  during  processing,  At  present,  there  are  30  types  of  error 
messages  that  can  be  printed  as  a  result  of  the  third  octave  process. 

(6)  Third  octave  data  dump.  This  program  provides  a  sequential  listing  of  all 
ten  second  averages  obtained  during  the  hydrophone  third  octave  run.  It  is 
used  to  cross-check  the  results  of  the  third  octave  process  prior  to  re¬ 
leasing  the  digital  tape  for  further  processing.  The  ten  second  averages 
are  arranged  in  one  minute  blocks  and  tagged  with  day  number,  hours  and 
minutes,  gain  state,  and  transfer  function. 


4.  INTERNAL  COMPUTER  PROCESSING  OF  AMBIENT  NOISE  DATA 
4. 1  The  Acoustic  System  Equation 

The  computer  inputs  data  from  one  hydrophone  via  15  bands  (15  through 
300  Hz)  through  the  third  octave  filters  and  the  RMS  detector  to  form  10 
second  average  spectrum  levels  for  each  band.  The  spectrum  levels  are  then 
corrected  to  microbars  using  the  following  correction  equation: 

(LS) j  »  X.  -  Sh  -  Gc  -  (delta) .  -  (BWC)j  -  (T' 5Q)  (4.1-1) 

-  is  corrected  band  spectrum  levels  re: 
raicrobar/volt/IIz  for  band  j . 

-  is  uncorrccted  computer  ambient  noise 
level  for  band  j. 


where  (LS) . 

J 


X. 

J 


is  hydrophone  sensitivity. 
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G 


-  is  gain  setting  of  data  amplifier. 


c 

(delta) j  -  is  a  correction  factor  for  band  j  to 
account  for  the  frequency  response  of 
a  given  hydrophone  within  the  dynamic 
range  of  15  to  300  Hz. 

(BWC)  -  is  the  bandwidth  correction  for  band  j. 

J 

CT * ^q)  -  is  the  system  transfer  function  obtained 

during  six  hour  calibration  cycle. 

The  system  transfer  function  is  computed  from  the  50  and  200  Hz  signals  in¬ 
jected  every  six  hours  for  five  minutes.  The  program  stores  the  50  and  200  Hz 
10  second  averages  for  four  minutes  and  from  previous  knowledge  of  the  levels 
injected  along  with  the  appropriate  gain  settings  of  the  hydrophone,  determines 
the  following  equation: 


50 


<X' 


50 


-  (C1), 


Gc* 


(4.1-2) 


where  X'5Q  -  is  tr.e  ensemble  average  of  the  50  Hz  level 

obtained  during  the  5  minute  calibration 
sequence . 

C'  -  is  the  injected  50  Hz  level  applied  during 

minute  n. 

Gc'n  -  is  the  data  amplifier  gain  state  level  for 
minute  n. 

Typical  level  injection  and  gain  states  for  a  calibration  cycle  are: 


Minute 

2: 

-41 

db 

level. 

40 

db 

gain 

Minute 

3: 

-31 

db 

leve 1 , 

30 

db 

gain 

Minute 

4: 

-21 

db 

level. 

20 

db 

gain 

Minute 

5: 

-11 

dl> 

level , 

10 

db 

gain 
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Note  that  minute  1  of  the  calibration  cycle  is  not  used.  The  first  minute 
is  used  to  allow  the  data  amplifiers  to  seek  the  proper  gain  for  the  first 
level  injection. 

fhe  corrected  15  bands  of  data  in  the  form  of  10  second  averages  are  then 
written  out  onto  digital  magnetic  tape  arranged  in  seven  one-minute  blocks 
as  the  record  length.  Each  minute  represents  five  integrations  tagged  with 
day,  hour,  minute,  data  amplifier  gain  state,  and  transfer  function  obtained 
during  the  six  hour  caUbration  interval.  It  should  be  noted  that  one  trans¬ 
fer  function  is  used  during  the  entire  process  period  of  one  hydrophone.  The 
current  transfer  function  computed  during  any  six  hour  interval  is  stored  in 
the  one  minute  data  record  but  not  used.  Dynamic  changes  in  the  system  af¬ 
fecting  the  overall  transfer  function  should  be  examined  carefully  and,  if 
necessary,  corrections  then  would  be  applied  in  later  processes  (when  the 
clean  third  octave  tape  is  made). 

The  general  equation  is  determined  from  the  acoustic  system  below: 


Hydro¬ 

phone 


Hyd.  Pre- 
Pre  Amp  Amp 


Geotech 

Recorder 


Sangamo 

Playback 


General  Radio  8 
1/3  Octave  Analyzer 


Bandwidth  Corr. 
/  Computer 


o  o 

-  o  ■ 

- 

■> 

U 


i'j 


2')j  D'|  (BWC)j  Xj 


For  a  given  input  signal  (Ls)  from  the  hydrophone,  the  computer  stores  the 

average  spectrum  level,  X.,  in  db  for  each  third  octave  band  selected.  The 

3 

equation  representing  this  system  is: 

X  =  (L  )  ♦  S  ♦  (G  )  ♦  (T  )  ♦  (T,).  ♦  D  *  (BUG), 

j  s  j  n  ci  i  )  J  J 


(4.1-5) 
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where  =  computer  raw  value  for  band  j 

=  hydrophone  sensitivity  @  SO  Hz  (Temp,  and  pressure) 

Gc  =  gain  setting  of  data  amplifier 

=  transfer  function  of  hydrophone  preamp 
?2  =  transfer  function  of  recorder  and  playback  system 

D  =  transfer  function  of  third  octave  filter  bank  plus 

General  Radio  RMS  detector 

BKC  -  bandwidth  correction  factor  computed  from  10  log  BW, 
where  BW  =  bandwidth. 

For  purposes  of  the  playback  analysis,  the  transfer  functions  of  the  Geotech 
recorder,  Sangamo  playback,  and  General  Radio  detector  and  third  octave  filter 
are  combined  as  T. 

Tj  =  (1Vj  +  (Vj  "  Dj  (4.1-4) 

The  equation  becomes: 

X.  =  L.+S,  +  G.+T.*  (BWC)  .  (4.1-5) 

J  sj  h  cj  3  '3 

During  the  calibration  phase,  a  calibration  signal  at  a  fixed  frequency 
(50  and  200  Hz)  is  recorded  on  tape.  The  playback  system  is  then  represented  by: 

Geotech  Sangamo  General  Radio  & 


Combining  the  transfer  functions  of  the  Geotech  and  Sangamo  recorders  and  the 
General  Radio  Analyzer  yield  ar.  overall  transfer  function,  T'^,  during  the 
calibration  cycle.  This  can  be  expressed  as: 


T’.  =  (T/)  -  (T  ')■  ♦  D*  14.1-6) 

The  overall  equation  becomes: 

x' .  =  C* .  ♦  G'  ♦  T*. 
i  j  c  j 


(4.1-7) 
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where  C'  *  calibration  frequency  §  50  and  200  Hz  at  four  signal 

levels 

j  =  1/3  octave  bands  for  50  and  200  Hz  signals 

Solving  for  the  unknown  transfer  function  represented  by  the  Geotech,  Sa.igamo, 
and  General  Radio  yields 

T,50^X,50>  *  <C’n  -GCV  t4-1'8) 

where  n  =  1,  2,  3,  or  4  for  the  four  signal  levels  and  gain  settings. 

<  >  =  ensemble  average 

Since  the  calibration  signal  levels,  C1,  and  gain  settings,  G'c,  are 
known,  T*  represents  the  new  computed  transfer  function  of  the  recorder 
systems  for  a  given  calibration  period  at  the  reference  frequency  of  50  Hz. 

4.2  Pre-Process  Calibration 

Now,  to  obtain  a  relationship  between  the  transfer  function  of  the  over¬ 
all  system  at  the  50  Hz  reference  frequency  to  any  third  octave  center  fre¬ 
quency,  a  calibration  tape  must  be  prepared  and  played  through  the  system 
applying  a  calibration  signal  input  at  each  of  the  third  octave  center 
frequencies  from  10  to  300  Hz.  The  plot  of  figure  7  shows  the  transfer 
function  of  the  system  with  a  constant  signal  level  applied  to  the  input  of 
the  hydrophone  pre-amplifier  and  the  output  read  from  the  computer. 

Using  the  50  Hz  calibration  signal  as  a  reference,  it  can  be  seen  that 


at  any  given  frequency,  j,  the  overall  transfer  function  '  is: 

V  *  T'sc  * 

where  j  =  the  particular  third  octave  center  frequency.  Since  the  hydrophone 


sensitivity,  3^,  varies  as  a  function  cf  frequency,  the  A.  term  applys  a 
correction  factor  to  account  foT  the  frequency  response  of  a  given  hydrophone 


within  the  dynamic  range.  The  A^  term  is  then  an  array  of  values  for  each 
third  octave  center  frequency  applied  to  the  correction  equation. 
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Table  4  shows  a  typical  calibration  table  prepared  from  an  injection  of 
a  constant  voltage  at  the  center  frequencies  of  the  third  octave  bands  of 
interest. 

Returning  to  the  general  equation,  for  any  given  frequency  the  equation 
becomes 

(LS).  =  X.  -  S.  -  K  -  (T.).  -  T*.  -  (BWC).  (4.2-2) 

v  )  J  h  c  V)  3 

Substituting  equation  (4.2-1)  yields: 

(LS).  =  X.  -  Sh  -  Gc  -  T'so  -  A.  -  (BWC).  (4.2-3) 

Since  the  automatic  calibration  period  determines  the  transfer  function  at 
the  50  Hz  reference  frequency  every  six  hours,  equation  (4.1-8)  can  be  sub¬ 


stituted  for  T'so  to  give  the  general  equation: 


(LSJj  ■  X.  -  Sj,  -  Cc  -  4j  -  (BWC).  -  [X'so  -  C'n  -  C'cJ  (4.2-4) 
Equation  (4.2-4)  is  implemented  by  the  computer  program  and  converts  each  of 


the  third  octave  band  outputs  to  average  spectrum  levels  reference  to  a 


microbar. 


4.3  Automatic  Calibration 

The  program  performs  automatic  calibration  of  two  calibration  fre¬ 
quencies  once  every  six  hours.  The  two  mixed  frequencies,  50  and  200  hertz 
are  banc-passed  through  the  third  octave  filters  (bands  36  and  42  on  80:1 
playback)  and  detected  through  the  General  Radio  analyzer  using  a  10  second 
average. 

The  calibration  sequence  lasts  for  a  period  of  5  minutes.  The  program 
detects  the  cal  i brat  ion  phase  by  a  frequency  modulated  on  the  time  code 
carrier.  This  frequency  passes  through  a  narrow  band  filter  on  the  control 
function  indicator  board  and  into  the  computer  as  a  bit  indication.  Once 
this  bit  is  set,  the  program  enters  a  calibration  sequence  for  a  five  minute 
period,  the  duration  or  the  calibration  sequence.  The  mixed  50  and  20'*  her«.z 
frequencies  which  passed  through  tho  third  octave  filters  are  averaged  o.:  the 
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ACODAC  SYSTEM  "A"  CALIBRATION 


further  information 


NOT  ORIGINAL  DATA 

NOTE:  At  62,5  Hz  and  125  Hz  there  was  a  lot  of  60  Hz  pickup  and  the 
data  had  to  be  estimated  in  the  last  place. 

InPut  #4  was  disconnected  during  initial  portion  of  the  62.5  Hz  test 
NOTE:  Vj  =  input  to  hydrophone;  V'2  =  output  from  hydrophone;  V  =  output  from  amplifiers. 


ACODAC  SYSTEM  "B"  CALIBRATION 
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General  Radio  RMS  detector.  The  following  describes  the  sequence  of  events 
that  occur  for  the  five  minute  calibration  period. 

Minute  1  Idle  phase:  gain  settings  are  indeterminable. 

System  allows  one  minute  for  gain  to  be  settled 
in  preamplifier. 

Minute  2  Gain-1:  highest  gain  setting,  lowest  signal 
level.  Five  10-second  averages  collected  for 
one  minute. 

Minute  3  Gain-2:  next  lower  gain  setting,  next  higher  signal 
level.  Five  10-second  averages  collected  for  one 
minute . 

Minute  4  Gain-3:  next  lower  gain  setting,  next  higher 

signal  level.  Five  10-second  averages  collected 
for  one  minute. 

Minute  5  Gain-4:  lowest  gain  setting,  highest  signal 
level.  Five  10-second  averages  collected  for 
one  minute. 

At  the  start  of  the  calibration  period,  the  program  sets  up  a  one  minute 
wait  period  to  allow  for  the  settling  of  the  capsule  preamps  to  the  highest 
gain  setting  and  the  lowest  signal  level.  At  the  start  of  minute  2,  the 
program  starts  the  10  second  integration  cycle  of  the  General  Radio 
Analyzer.  At  the  completion  of  the  10  second  period,  the  program  tem¬ 
porarily  stores  the  integrated  values  of  the  50  and  200  llz  signals,  waits 
two  seconds  and  starts  the  next  integration  cycle.  This  process  continues 
for  four  minutes  storing  the  results  of  a  10-second  integration  period  for 
the  SO  and  200  llz  bands.  The  program  then  computes  the  mean  wG  for  each 
of  the  gain  settings  and  for  each  frequency.  Thus,  for  the  first  gain  setting 

of  minute  2  (gain  1),  the  mean,  y„,  is: 

0 
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5 

T, 

N=1 


(4.3-1) 


from  the  general  equation: 


S 

I 

N=1 


(4.3-2) 


where: 

B^  =  the  output  of  the  General  Radio  Analyzer  in  db 

j  =  Band  36  for  50  Hz  and  Band  42  for  200  Hz 

k  =  1,  2,  3,  4  representing  the  four  gain  settings  of 

the  preamplifier 

N  =5,  the  number  of  integrations 

Figures  8  and  9  illustrate  the  timing  of  the  integration  periods  for  the 
10  second  averages. 

The  purpose  of  the  calibration  is  to  determine  the  value  of  the 
transfer  function  of  the  Geotech  recorder  and  Sangamo  playback  unit  where  a 
given  input  signal  is  applied  to  the  data  amplifier  and  is  recorded  at  a 
known  gain.  If  any  dynamic  changes  have  occurred  between  the  recording  and 
playback  of  an  input  signal  from  the  hydrophone,  the  calibration  of  a  known 
signal  on  the  same  recording  shall  also  be  subjected  to  the  same  level 
change.  Sines  the  true  level  of  the  calibration  signal  is  known,  then  the 
ambient  noise  signal  coming  from  the  hydrophone  can  be  corrected  to  reflect 
this  dynamic  change.  Therefore,  the  true  hydrophone  signal  can  be  recon¬ 
structed  without  distortion  from  recording  and  playback. 

If  no  dynamic  changes  have  occurred  from  recording  to  playback,  or  if 
r.o  hardware  malfunctions  have  occurred,  the  four  calibration  levels  for  the 
two  frequencies  should  not  vary  more  than  ^  to  1^  db  through  the  entire 
range  from  the  expected  values  of  the  cal  levels  and  gain  settings.  If 
this  is  true,  then  all  20  values  (5  integrations  times  4  levels)  for  each 
frequency  will  be  averaged  together  to  obtain  one  value  which  will  then  be 


Figure  9.  Typical  Time  Spacing  for  50  &  200  Hz  Calibration  Signal 


-39 


used  in  the  calibration  equation  to  determine  the  transfer  function  of  the 
Geotech  and  Sangamo  recording  systems. 

However,  the  five  integration  results  for  each  level  are  first 
tested  to  observe  if  these  values  are  within  the  tolerance  range  of  the 
expected  value  (obtained  from  a  knowledge  of  the  gain  setting  and  the 
voltage  level).  If  the  five  integrated  values  for  any  given  level  fall 
beyond  this  tolerance  range,  then  the  particular  cal  level  will  be  suspect 
of  malfunction  and  excluded  fiom  the  running  average.  This  event  will  be 
flagged  in  the  program  and  the  values  printed  on  the  teletype.  This 
is  one  method  of  determining  hardware  electronic  malfunctions  in  the  ACO- 
DAC  capsule. 

Once  the  calibration  data  is  acquired  and  the  transfer  functions  com¬ 
puted,  a  calibration  summary  will  be  printed  on  the  teletype  for  purposes 
of  operator  examination.  See  Fig.  10  for  sample  printout. 

TYfo  tests  are  applied  to  the  calibration  data  to  insure  validity  of 
information  during  the  calibration  phase. 

(1)  The  results  of  the  five  integrations  within  the  minute  interval  are 
summed,  and  the  sum  is  compared  to  each  of  the  five  -  10  second  integrations 
of  that  interval.  If  any  of  the  10  second  averages  are  greater  than  the 
applied  tolerance,  in  this  case  1  db,  the  entire  one  minute  interval  is 
rejected.  If  the  interval  passes  the  tolerance  test,  the  mean  of  the  five 
integration  values  is  computed  and  is  used  in  the  computation  of  the  trans¬ 
fer  function. 

(2)  A  gain  check  is  then  made  to  insure  that  the  proper  gain  level  was 
applied  at  the  correct  minute  interval.  This  check  is  necessary  since  the 
cc’ibration  signal  level  injected  at  the  specific  minute  interval  forces  the 
gain  to  switch  to  the  proper  range  corresponding  to  the  level  of  the  signal. 
That  is,  an  injected  -41  db  cal  level  should  cause  the  amplifier  to  switch 


-40- 

to  the  highest  gain  range  of  40  db.  An  injected  level  of  -31  db  should 
correspondingly  cause  the  gain  to  switch  to  a  30  db  level.  If  the  computer 
program  were  not  in  sync  with  the  level  switching,  then  improper  gain 
settings  would  be  applied  to  the  cal  signal  levels  thus  resulting  in  an 
erroneous  transfer  function  value. 

During  the  early  development  stages  of  the  third  octave  analysis  pro¬ 
gram,  the  transfer  function  resulting  from  a  six  hour  calibration  sequence 
was  injected  into  the  system  equation  when  converting  raw  ambient  noise 
data  to  microbar  levels.  The  choice  was  left  to  the  operator  during  the 
playback  phase  to  decide  whether  to  allow  the  new  transfer  function  to  be 
inserted  into  the  system  equation,  or  to  keep  the  previously  computed  trans¬ 
fer  function.  The  operator's  decision  was  based  on  a  quick  analysis  of  the 
calibration  summary  printout  (now  called  Post  Analysis  Summary  Printout, 

Fig.  10).  It  was  felt  that  the  analysis  of  the  data  required  more  acousti¬ 
cal  knowledge  of  the  system  than  might  be  available  at  the  level  of  the 
operator  to  make  valid  decisions.  The  program  was  then  changed  to  print 
out  the  computed  transfer  function  during  the  calibration  phase,  store  the 
transfer  function  in  the  current  digital  output  record  and  then  proceed 
automatically  ahead  digitizing  the  raw  ambient  noise  data. 

The  current  program  prints  out  the  following  calibration  sequence 
message: 

1200  CALIBRATION  IN  SEQUENCE 
TRANSFER  FUNCT:  OLD  =  58.2  NEW  =  38.5 
TOL  USED  =1.0 
CALIBRATION  COMPLETE 

The  printout  shows  the  operator  the  previously  computed  six  hour  calibra¬ 
tion  transfer  function  as  well  as  the  current  computed  value  along  with  the 
value  of  the  tolerance,  in  db,  used  to  accept  or  reject  signal  levels  in 
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computing  the  mean  of  the  levels.  If  all  of  the  four  minute  intervals  (20 
data  points  total  for  each  frequency)  of  either  the  50  or  200  hertz  signal 
levels  are  unable  to  pass  the  tolerance  test,  then  the  program  will  print: 

NO  USABLE  50  HZ  DATA 
NO  USABLE  200  HZ  DATA 
CALIBRATION  COMPLETE 

If  the  above  condition  occurs,  the  program  will  store  the  last  good  computed 
transfer  function  in  the  output  data  record. 

4.3.1  Obtaining  the  Initial  Value  Transfer  Function 

The  initial  value  for  the  system  transfer  function  is  inserted 
in  the  program  during  the  initialization  phase  and  prior  to  the 
playback  phase.  This  value  is  obtained  by  actually  playing  the 
ambient  noise  tape  through  the  system  at  various  six  hour  cal¬ 
ibration  intervals  through  the  whole  deployment  to  obtain  a  mean 
value  for  the  transfer  function.  This  value  is  then  inserted  in 
the  program  and  used  through  the  entire  acquisition  phase. 

Plans  are  being  implemented  to  record  the  pre-calibration 
signal  levels  directly  on  the  deployment  data  tape  prior  to 
system  deployment  along  with  at  least  two  calibration  sequences. 
This  will  preclude  the  necessity  of  preparing  a  pre- calibration 
tape  on  a  separate  recording  which  has  been  done  in  the 
past.  There  will  still  be  a  requirement  to  generate  the  same 
calibration  data  on  a  recording  after  the  system  has  been  re¬ 
covered.  Usually,  there  is  insufficient  tape  remaining  on 
data  tapes  to  make  this  type  of  recording. 

4.3.2  Selecting  the  Transfer  Function  from  Either  the  50  or  200  Hz 
Signals 

As  was  previously  stated,  the  transfer  function  can  be 
computed  using  the  50  hertz  or  200  hertz  calibration  levels. 
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Normal  ly  the  50  llz  signal  is  \sed  if  at  least  one  of  the  four 
minute  intervals  passes  he  tolerance  test.  If  all  the  four 
minute  50  llz  intervals  fail  to  pass,  the  200  llz  signal  is  used 
for  the  transfer  function  computation  providing  again  th_t  at 
least  one  of  the  four  minute  200  Hz  intervals  passes  the  tol¬ 
erance  test. 

If  the  50  Hz  signal  is  used,  the  transfer  function,  TF, 
computed  is:  . 

£  C*50^n 

TF  =  T'50  =  -  (4. 3. 2-1) 


where  n  =  the  no.  of  minute  intervals  used  in  the  mean;  (n  <_  4) . 
When  the  200  Hz  signal  level  is  used,  the  transfer  function 
must  compensate  for  the  difference  in  frequency  response  be¬ 
tween  the  50  Hz  reference  and  the  200  Hz  signal.  This  is 
done  by  subtracting  the  delta  function  value  (A^)  at  200  Hz 
from  the  computed  mean  using  the  200  Hz  signal,  or 


TF  * 


T' 

200 


\  ^200^ 

n=l 


(4. 3.2-2) 


where  X  =  50  Hz  mean  of  the  minute  interval  n,  n  =  1  to  4 
^200  *  200  l!z  mean  the  m*nute  interval  n 
A^qo  =  Term  applied  in  db  for  the  difference  in  frequency 
response  between  the  50  Hz  reference  and  200  Hz 
signal  injected  through  system  during  pre-process 


calibration  run. 
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The  following  table  gives  the  calibration  signal  input 
levels  for  each  of  the  data  amplifier  gain  settings  during  the 
calibration  sequence. 


TABLE  5 


CALIBRATION  SIGNAL  INPUT 

LEVELS 

Minute  During 

Cal  Sequence 

Input  Level  to 
50  Hz 

Data  Amplifiers 

200  Hz 

Amplifier 

Gain 

1 

0.0  VRMS 

0.0  VRMS 

N/A 

2 

.009V  =  -41DB 

.009V  *  -41DB 

♦40  DB 

3 

.028V  -  -31DB 

.028V  *  -31DB 

+30  DB 

4 

.090V  =  -21DB 

.090V  =  -21DB 

+20  DB 

5 

.283V  =  -ilDB 

.283V  =  -11DB 

+10  DB 

(For  the  1973  field  operations,  pre-emphasis  at  the  higher 
frequencies  was  added  to  the  two  ACODAC  Mod  I  systems.  The 
new  values  for  the  200  Hz  calibration  ari  -46  db,  -36  db, 

-26  db,  and  -16  db  rather  than  the  above.) 

4.4  Post  Calibration  Summary 

The  post  calibration  summary  printout  is  used  to  evaluate  the  quality 
of  calibration  data  used  in  the  computation  of  the  system  transfer  function. 
The  evaluation  is  done  subsequent  to  i  hydrophone  third  octave  process  run 
md  is  used  primarily  to  assist  in  the  detection  and  analysis  of  hardware 
problems  that  may  have  occurred  in  the  system  either  during  the  deployment 
or  playback  phase.  Figures  10  and  11  show  typical  printouts.  Figure  11  shows 
obvious  problems  from  the  number  of  99. 9' s  that  appear  in  various  columns. 

A  considerable  amount  of  information  is  presented  to  assist  the  analyst 
in  determining  system  performance.  For  example,  the  mean  of  each  minute  in¬ 
terval  is  given  for  both  the  50  and  200  hertz  signals  along  with  the  computed 
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POST  CALIBRATION  NUMMARY  PRINTOUT 
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Figure  10.  Post  Calibration  Summary  Printout 
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POST  CALIBRATION  SUMMARY  PRINTOUT 
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Figure  11.  Post  Calibration  Summary  Printout 
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transfer  functions  for  each  minute  interval.  That  is: 

TPS  =  SO  Hz  transfer  function 

TP2  =  200  Hz  transfer  function 

The  injected  gain  for  each  minute  interval  is  given.  If  the  proper  gain 
level  is  rot  injected  or  detected  during  playback  for  the  right  minute 
interval,  the  program  prints  a  99.9  in  the  associated  minute  column  for 
both  the  gain  and  the  transfer  function,  since  the  transfer  function  for 
the  particular  minute  interval  cannot  be  computed  without  the  correct  gain. 
The  gain  error  printout  indicates  whether  the  system  had  a  gain  redundancy 
error  in  decoding  the  gain  bits.  If  this  occurs,  a  negative  gain  is  printed 
on  this  line  for  any  minute  interval  in  error. 

Along  with  the  mean  for  the  50  and  200  hertz  signals,  information  is 
provided  that  show  the  amount  each  integration  interval  varies  from  the 
tolerance  used.  Since  the  levels  for  each  of  the  five  10-second  inte¬ 
gration  intervals  should  be  the  same  for  each  minute,  with  only  slight 
variations,  any  large  variations  of  levels  around  the  mean  value  indicate 
possible  problem  areas  in  the  IPV  electronics  system.  Also,  it  is  noted 
that  the  system  transfer  function  can  be  computed  using  the  200  hertz 
signal  if  it  is  found  that  the  50  hertz  signal  levels  exceed  tolerance 
limits.  If  the  200  hertz  signal  is  used,  an  important  factor  in  the  compu¬ 
tation  of  the  transfer  function  is  the  -'election  of  the  delta  function  value 
for  the  200  hertz  signal.  Since  the  50  hertz  signal  is  the  reference, 
and  since  the  delta  function  value  accounts  for  the  system  frequency  re¬ 
sponse  at  200  hertz,  the  computed  transfer  function  would  be  in  error 
if  the  delta  function  were  in  error. 
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5.  THIRD  OCTAVE  ANALYSIS  PROGRAM  DESCRIPTION 
5.1  General  Specifications 

The  playback  system,  figure  12,  is  designed  to  be  controlled  by  the 
digital  computer  using  a  single  hydrophone  channel  input  connected  to  the 
system  for  processing  through  the  GR  1921  Real  Time  Analyzer  (1/3  octave). 
Analog  data  from  the  hydrophone  channel  is  digitized  through  the  analyzer 
and  provided  to  the  computer  for  further  processing.  The  computer  executes 
a  sequence  of  events  to  convert  the  raw  ambient  noise  data  from  15  third 
octave  bands  represented  by  10  second  average  spectrum  levels  to  ambient 
noise  levels  reference  to  a  microbar.  The  events  are  as  follows: 

(1)  Initiate  and  stop  integration  cycles  in  the  analyzer  at  specific 
intervals  of  time  dictated  by  the  analog  data  timing  sequence. 

(2)  Shift  15  channels  of  data  from  outputs  of  the  analyzer  accu¬ 
mulator  sequentially. 

(3)  Detect  and  act  up  in  control  signals  coming  from  the  playback 
channel  which  represent  (a)  start  of  data  signals,  (b)  overload  signals, 

(c)  calibration  signals,  and  (d)  end  of  data  signals. 

(4)  Read  tine  in  days,  hours,  and  minutes  each  minute  from  the  time 
code  read. 

(5)  Read  state  codes  from  the  time  code  reader. 

(6)  Convert  and  correct  the  incoming  data  to  average  spectrum  level 
reference  to  a  microbar. 

(7)  Identify  and  tag  all  values  from  any  integration  cycle  that  may 
contain  had  data. 

(8)  Perform  automatic  system  calibration  every  six  hours  and  store 
the  results  on  hard  copy  printout  and  on  digital  magnetic  tape. 

(9)  Provide  error  diagnostic  message  printout  on  system  performance. 


Figure  12.  Third  Octave  Ambient  Noise  Analysis  System 
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(10)  Plot  the  corrected  spectrum  levels  of  selected  third  octave  bands 
on  an  analog  plot  device. 

(11)  Dump  7-minute  block  records  of  data  onto  digital  magnetic  tape. 

These  events  are  controlled  by  a  program  employing  a  real  time  execu¬ 
tive  which  operates  in  a  foreground -background  mode  in  order  to  schedule 
the  various  internal  sequences.  The  foreground  section  of  the  progrz.7.  re¬ 
sponds  to  all  active  interrupts  from  the  devices.  The  priority  interrupt 
system  is  used  to  schedule  the  servicing  of  each  device  according  to  the 
rate  at  which  information  is  flowing  to  and  from  the  computer.  Higher  speed 
devices  are  given  higher  priority  servicing  over  slower  speed  devices. 

The  background  mode  performs  the  necessary  conversion  and  correction  of  data, 
initiates  or  terminates  various  program  modes  of  operation,  and  monitors 
the  operator  console  for  command  inputs. 

5.2  Computer  Timing 

A  one  pulse  per  second  (1  pps)  signal  is  supplied  to  the  computer  from 
the  IRIG-C  time  code  track  of  the  ACODAC  Analog  tape  by  the  time  code  genera¬ 
tor.  This  signal  is  the  basis  for  all  the  computer  timing  and  is  used  to 
schedule  all  events  such  as  starting  the  integration  cycle,  printing  out 
error  messages,  and  plotting  analog  d,.ta. 

rhe  1  pps  signal  to  the  computer  appears  in  playback  time  as  a  signal 
generated  each  l/'80th  of  a  second,  or  80  pps  (.0125  secs).  This  is  because 
the  analog  tape  recorded  at  15/160  inches  per  second  is  played  back  at  7.5 
inches  per  second  -  a  speed  up  of  80. 

5.3  Integration  Timing  Control 

Part  of  the  one  pulse  per  second  interrupt  routine  is  used  to  start 
successive  10  second  integrations  every  12  seconds  by  sending  out  a  control 
pulse  to  the  General  Radio  1921  Real  Time  Analyzer.  The  analyzer  in  turn 
uses  its  internal  clock  to  start  the  10  second  integration  rate.  During 


this  time,  the  analyzer  operates  by  sampling  each  input  channel  from  the 
GR  1925  Multifilter  (a  part  of  the  GR  1921  system)  and  supplies  these 
samples  to  the  digital  circuits  that  are  time  shared  on  all  channels.  The 
detector  simultaneously  computes  the  rras  level  for  each  filter  channel 
by  taking  128  samples  during  the  1/8  second  integration  period  (10  seconds 
in  hydrophone  time) .  The  samples  are  converted  to  digital  binary  numbers 
and  squared.  The  squared  values  are  accumulated  in  a  memory  register  where 
spaces  are  provided  for  all  channels.  At  the  end  of  the  integration  period, 
the  sum-of- squares  value  is  converted  to  decibels  for  output  presentation. 

At  this  time  the  computer  receives  an  interrupt  enabling  transfer  circuits 
to  strobe  the  data  into  the  computer.  This  process  is  repeated  continu¬ 
ously  while  the  program  is  in  the  acquisition  phase. 

During  this  process,  five  integrations  are  accumulated  every  minute 
and  stored  for  conversion  and  correction  to  microbar  levels.  Since  the 
data  amplifiers  can  shift  gain  with  respect  to  the  level  of  the  input 
ambient  noise  signal,  pr  >e  synchronization  must  be  maintained  between 
the  integration  process  .nd  the  IRIG-C  time  code  frame.  This  is  necessary 
since  the  input  signal  to  the  data  amplifiers,  during  the  deployment  phase, 
is  integrated  over  a  minute  period  to  determine  the  gain  range  for  the  am¬ 
plifier.  At  the  end  of  the  minute  frame,  if  the  signal  level  exceeds  the 
dynamic  range  of  the  present  gain  setting,  the  electronics  in  the  instru¬ 
ment  pressure  vessel  (IPV)  automatically  switches  the  amplifier  gain  range 
and  at  the  same  time  injects  the  new  gain  setting,  represented  by  two  bits 
per  hydrophone,  into  the  time  code  hit  stream  in  the  control  function  posi¬ 
tion  of  the  time  frame. 

Since  the  program  reads  the  time  code  and  gain  settings  at  the  start 
of  the  new  minute  frame,  the  program  now  applies  the  new  gain  settings  to 
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the  new  set  of  five  integrations  which  will  be  completed  at  the  end  of  the 
new  minute  frame. 

However,  if  synchronization  is  not  maintained  between  the  five  integra¬ 
tions  and  the  reading  of  the  gain  bits,  an  improper  gain  correction  will  be 
applied  to  the  raw  ambient  noise  data.  This  phenomenon  will  appear  as  a 
step  function  in  the  datr.  if  observed  on  the  Brush  analog  plot. 

Many  conditions  have  been  observed  which  cause  program  timing  to 
loose  synchronization.  Time  code  bit  drop-out  is  a  possible  occurrence 
resulting  in  the  computer  dropping  a  clock  pulse  and  thus  losing  the  60 
second  count  during  the  minute  frame.  Although  redundant  gain  bits  were 
added  to  the  electronics  for  each  hydrophone  amplifier  to  detect  gain  bit 
errors,  if  a  gain  bit  error  is  detected,  the  program  must  reject  the  present 
gain  setting  and  use  the  previous  minute  gain  setting.  This  of  course  can 
result  in  erroneous  spectrum  levels  from  the  five  integrations. 

Therefore,  it  should  be  noted  that  no  foolproof  method  has  been  de¬ 
vised  to  account  and  correct  for  all  possibilities  occurring,  although  the 
program  does  make  elaborate  synchronization  checks.  It  is  necessary  for  a 
trained  analyst  to  examine  the  data  output  to  verify  the  computations. 

5.4  Playback  Process  Ti  a  and  Throughput  Rate 

Most  ACODAC  moorings  have  been  deployed  for  periods  of  10  days 
(The  maximum  data  acquisition  time  is  actually  10  2/3  days  -  7200  feet  of 
magnetic  tape  at  15/160  inches  per  second).  Because  of  the  long  deploy¬ 
ment  duration  and  the  use  of  slow  recording  speeds  of  the  analog  recorder 
(15/160  ips) ,  an  80:1  time  compression  ratio  has  been  used  to  accelerate 
ambient  noise  processing. 

In  the  accelerated  process,  each  second  of  hydrophone  time  is  repre¬ 
sented  by  .0125  seconds.  The  10  second  (hydrophone  time)  integration  period 
would  then  be  completed  in  .125  seconds. 
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There  is  a  minimum  of  150  msecs  to  read  in  15  bands  (16  bits/band) 
of  data  and  convert  and  process  the  data  before  a  new  integration  can  be 
started.  This  period  is  broken  down  as  follows: 

Integration  period  .12S  seconds 

Time  between  integrations  .025  seconds 

Total  time  .150  seconds 


The  throughput  rate  then  becomes: 


15  words 
.15  secs 


100  words/sec. 


It  can  be  seen  that  at  this  slow  throughput  rate  there  remains  consider¬ 
able  time  in  the  background  to  accomplish  many  job  tasks  such  as  data  con¬ 
version,  correction,  analog  plotting,  and  error  detection. 

5.5  Digital  Output  Magne.ic  Tape  Data  Records 
5.5.1  Data  Recording 

Data  records  represent  a  7  minute  (hydrophone  time) 
period.  Each  record  j«=  1092  words  long  and  are  written  at 
800  bpi.  As  an  example,  assume  there  are  a  total  of  10  days 
of  recorded  data  on  the  analog  ambient  noise  tape.  Since  this 
period  represents  10  days  X  24  hours/day  X  60  min/hour  = 

14,400  min,  the  total  number  of  data  records  per  hydrophone 
channel  is: 


14400 


=  2057  records 


At  a  recording  density  of  800  bpi,  each  record  occupies: 
1092/800  =  1.36  inches 


plus  an  inter- record  gap  of  .75  inches,  or  a  total  record 
length  of  2.11  inches. 
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The  total  length  recorded  on  tape  per  hydrophone  will  be: 

2057  recs  X  2.1  inches  =  4320  inches 
or  approximately  400  feet. 

At  this  rate,  approximately  6  hydrophone  data  channels 
could  be  recorded  on  one  reel. 

In  practice,  one  hydrophone  is  recorded  on  a  1200  foot 
reel  since  the  data  tape  is  then  submitted  to  other  computer 
processes  for  editing. 

5.5.2  Data  Record  Format 

Data  records  are  arranged  with  a  general  header  label, 
table  6  of  120  words  written  at  the  beginning  of  the  tape 
representing  vital  statistics  of  the  hydrophone  being  processed 
and  containing  pertinent  information  with  regards  to  the  ACO- 
DAC  deployment. 

Each  succeeding  record  represents  actual  data  arranged  in 
a  1092  word  record.  See  table  7.  Each  minute  of  the  data 
record  contains  a  six  word  header  identifying  the  date-time 
group,  data  type  (0  for  ambient  noise  data,  1  for  calibration 
data),  gain  setting  (2  bits),  and  transfer  function  computed 
for  the  pu  .^ar  six  hour  interval.  The  header  is  then  fol¬ 
lowed  by  five  integrations  for  that  specific  minute  arranged  in 
successive  15  bands  for  each  integration.  Each  integration 
represents  a  10  second  average  for  each  of  15  bands  in  db 
re:  microbar.  This  procedure  is  repeated  for  the  next  successive 
six  minutes.  Therefore,  for  a  seven  minute  record,  there  are 
five  integrations  per  minute,  or  a  total  of  35  integrations. 


TABLE  6 


ACODAC  HEADER  LABEL 


Word 

Variable 

Name 

Function 

Field 

Type 

1,2 

MOOR 

Mooring  type:  A=0,  B=1 

2  alpha 

A  format 

3,4 

IDEPLO 

Deployment  number 

2  alpha 

A  format 

5,6 

ISIIIP 

Launch  ship 

2  alpha 

A  format 

7,8 

I  CRUZ 

Cruise  number 

2  alpha 

A  format 

9,10 

ISTA 

Hydrophone  station  (1-6) 

2  alpha 

real 

11,12 

HYSER 

Hydrophone  serial  number 

4  digit 

real 

13,14 

DEPTH 

Hydrophone  depth 

4  digit 

real 

15,16 

HYCAL 

Hydrophone  sensitivity 

(-xx.x) 

real 

17,18 

TPRIM 

Pre- calibration  transfer 
function 

(xx.x) 

real 

19,20 

DDAY 

Deployment  launch  day  number 

(xxx.) 

real 

21,22 

DYLAR 

Deployment  launch  year 

(xx.) 

real 

23,24 

SDAY 

Recorder  start  time:  day 
number  (GMT) 

(xxx.) 

real 

25,26 

SI10UR 

Recorder  start  time:  hours 
(GMT) 

(xx.) 

real 

27,28 

SMIN 

Recorder  start  time:  minutes 
(GMT) 

(xx.) 

real 

29,30 

GSYNC 

GSYNC*1  means  time  recorded 
with  data  output  which  has 
been  converted  to  GMT. 

GSYNC=0  means  time  is  in 
relative  time. 

(x.) 

real 

31,32 

EDAY 

Recovery  day  number 

(xxx.) 

real 

33,34 

EGWT 

Recovery  time  (GMT) 

(xxxx.) 

real 

35,36 

AN1NB 

Minimum  1/3  octave  band 
displayed 

(xx.) 

real 

37,38 

AMAXB 

Maximum  1/3  octave  band 
displayed 

(xx.) 

real 
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ACODAC  HEADF.k  LABEL  (cont.) 

Variable 


Word 

Name 

Function 

Field 

Type 

39,40 

ONTIM 

Capsule  recorder  on-time 

(xx.) 

-1.A.  - 

real 

41,42 

DUTY 

Recorder  duty  cycle 

(.10-1.0) 

real 

43,44 

i  ZONli 

Local  time  zone  (at 
deployment) 

±x. 

real 

45,46 

DRIFT 

Time  drift  constant 

1  E  -6 

real 

47,48 

DLAT 

Deployment  latitude:  Degrees 

xx. 

real 

49,50 

XL  AT 

Deployment  latitude:  minutes 

XX.  x 

real 

51,52 

DLON 

Deployment  longitude:  Degrees 

XX. 

real 

53,54 

XLON 

Deployment  longitude:  minutes 

XX.  X 

real 

55, S4 

BEC( 15) 

Band  equalization  correction 

15  bands 

XX.  X 

real 

85,114 

Di:LTA(15) 

P re- calibration  delta 
function  15  bands 

XX. X 

real 

115,116 

ATA'M. 

/Analog  tape  ID 

xxxxx. 

real 

117,118 

DTAPE 

Digital  tape  ID  (data) 

XXX. 

real 

] 19, 120  GTAPo  Digital  tapelD  (gains) 


XXX. 


real 
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DATA  RECORD  FORMAT  ACODAC  OUTPUT  FROM  ACDI 


Digital  Magnetic  Tape  Record  of  analog  to  digital  conversion  process  using 
General  Radio  1921  Analyzer. 

Word  Variable^  Type 

Minute  #1 

1  Day  Integer 

2  Hour  ft  Minute  (24  hour  clock)  Integer 

3  Data  Type  (0  or  1)  Integer 

4  Gain  setting  (0,  1,  2,  or  3)  Integer 

5,o  Transfer  Function  T'  Floating  Point 

7-150  Integrations  #1-5  Floating  Point 

(5  integrations  per  minute  for 
15  third  octave  bends),  75  values 
total.  Each  10  second  integration 
period  results  in  15  band  values 
in  db  re:  1  volt /microbar. 

Minute  #2  through  #7 
Same  format  as  above. 

Total  number  of  words  =  1092 


where: 

Day.  hour,  and  minute  are  read  from  the  time  code  reader. 

Data  type  equals  1  for  calibration  data  (5  minutes  per  6  hours)  or 
0  for  acoustic  data. 

Transfer  function  (T)  is  computed  from  the  calibration  data  during 
the  calibration  phase. 

Gain  settings  are  read  from  the  state  codes  by  the  Time  Code  Reader 
and  are  either  0,  1,  2,  or  3  representing  10,  20,  30,  or  40  db 
respectively.  A  negative  integer  (-1,  -2,  -3)  indicates  that 
the  previous  1  minute  gain  setting  was  used.  Since  a  "-0"  cannot 
be  represented,  a  gain  of  10  db  (both  b±*s  0)  will  be  represented 
by  a  -32767  which  is  equivalent  to  setting  bit  15,  sign  bit,  to  one. 

Integrations  1  through  5  represent  the  average  spectrum  level  in  db  of 
15  filter  bands.  This  represents  a  total  of  75  values  for  the 
5  integration  periods,  or  150  words  ir.  floating  point  since  the 
values  can  range  from  35  to  -95  db. 

The  above  format  is  repeated  for  each  minute  interval  for  7  minutes.  The 
entire  record  length  is  thus  representative  of  a  7  minute  period  and  is  re¬ 
corded  as  109*.  words  on  magnetic  tape  for  each  record.  A  double  end-of-filc  is 
written  at  the  end  of  t he  deployment. 
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If  signal  overload  occurs  during  any  integration  period, 
the  program  will  replace  all  actual  band  values  with  a  99.9 
for  all  15  bands. 

A  double  end  of  file  is  inserted  at  the  end  following  the 
last  data  record. 

5.6  Program  Modes 

The  program  modes  of  operation  are: 

IDLE  (ID) 

ACQU I S IT  I ON  (AQ) 

INITIALISATION  (IN) 

TEST  (TE) 

TERMINATION  (EN) 

5.6.1  IDLE  Mode 

The  IDLE  mode  is  used  as  an  interim  mode  between  acquisi¬ 
tion,  test,  or  initialization.  During  this  phase  all  active 
equipment  interrupts  are  turned  off,  pointers  and  counters  re¬ 
set,  and  the  program  enters  a  wait  loop  for  an  operator  command. 
No  clock  timing  is  available  during  IDLE  since  all  program  tim¬ 
ing  is  obtained  while  the  playback  tape  is  in  motion.  Normally 
the  tape  unit  is  turned  off  during  IDLE. 

5.6.2  ACQUISITION  Mode 

The  ACQUISITION  phase  is  selected  to  start  the  ambient  noise 
data  processing  of  one  hydrophone  through  the  1921  General  Radio 
Analyzer  system.  ACQUISITIGN  is  the  main  program  loop  which 
schedules  the  digitizing  and  processing  of  the  ambient  noise 
data,  services  all  active  interrupts  from  the  equipment,  pro¬ 
vide^  diagnostic  error  message  summaries  to  the  operator,  and 
plots  the  corrected  third  octave  frequency  bands.  While  in 
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this  mode  a  digital  magnetic  tape  is  written  containing  the  10 
second  averages  corrected  to  microbar  levels. 

Continuous  data  integrations  through  the  GR  analyzer  will 
be  sequericed  as  long  as  this  mode  is  selected  and  the  analog 
tape  is  in  forward  motion. 

5.6.3  INITIALIZATION  Mode 

This  mode  is  used  to  enter  various  initialization  parameters 
to  the  program  via  the  teletype  keyboard.  These  data  are  en¬ 
tered  and  stored  in  memory  and  also  printed  as  a  header  label 
record  on  the  magnetic  digital  data  tape. 

Table  12  shows  a  list  of  all  entries  required  to  complete 
the  initialization  phase.  The  operator  has  an  option  to  go 
back  and  correct  any  particular  line  if  the  entry  was  made  in 
error  prior  to  the  information  being  written  on  tape. 

5.6.3. 1  Data  Array  entries 

Two  important  data  arrays  are  required  as  input 
during  the  initialization:  (1)  bandwidth  corrections  and 
(2)  delta  functions. 

(1)  The  bandwidth  corrections  are  applied  as  an 
arr-  of  numbers  for  the  effective  frequencies  of  interest 
being  processed.  Table  8  shows  a  table  of  third  octave 
frequencies  being  used  along  with  the  bandwidth  and  the 
bandwidth  correction.  The  bandwidth  correction  (BWQ  was 
computed  using  the  equation: 

BWC  *  10  log  Bandwidth  (BN) 

Table  9  shows  the  center  frequencies  for  a  standard  set 


of  third  octave  bands. 
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TABLE  8 


SPECTRUM  LEVEL  CORRECTIONS  FOR  1/3  OCTAVE  FILTERS 


Frequency 

10 

12.5 

16 

20 

25 

32 

40 

50 

63 

80 

100 

125 

160 

200 

250 

320 

•.00 

500 

630 

800 

1,000 

1,250 

1,600 

2,000 

2,500 

3,150 

4,000 

5,000 

6,300 

8,000 

10,000 


BW 

2.31 

2.77 

3.69 

4.62 

5.76 

7.37 

9.23 

11.50 

14.50 

18.40 

23.05 

27.7 

36.9 

46.2 

57.6 

73.7 

92.3 
115.0 
145.0 
184.0 
230.5 
277.0 
369.0 
462.0 
576.0 
739.0 
923.0 

1,150.0 

1,450.0 

1,840.0 

2,305.0 


10  log  (BW) 

3.6 

4.6 

5.6 

6.6 

7.6 

8.6 

9.6 

10.6 

11.6 

12.6 

13.6 

14.5 

15.6 

16.6 

17.6 

18.6 

19.6 

20.6 
21.6 
22.6 

23.6 

24.6 

25.6 

26.6 

27.6 

28.6 

29.6 

30.6 

31.6 

32.6 

33.6 


V 


I 


(»(>. 


TABLE  9 


CENTER  FREQUENCIES  FOR  THE  1/3  OCTAVE  BANDS 


Standard  Band 

Center 

Frequency  Mzl 

Standard  Band 

Center 

5 

Frequency  (Hz) 

3.162 

28 

631.00 

6 

3.9?  1 

29 

796.30 

7 

S.022 

30 

1000.00 

8 

6.310 

31 

1259.00 

9 

7.944 

32 

1585.00 

10 

10.00 

33 

1995.00 

11 

12.59 

34 

2512.00 

12 

15.85 

35 

3162.00 

13 

19.95 

36 

3981.00 

14 

25.16 

37 

5010.00 

15 

31.63 

38 

6310.00 

15 

39.81 

39 

7944.00 

17 

50.22 

40 

10.00  KHz 

18 

63.10 

41 

12.59 

19 

79.44 

42 

15.85 

20 

100.00 

43 

19.95 

21 

125.90 

44 

25.22 

22 

158.50 

45 

31.62 

23 

19!).. M) 

4<> 

59.81 

24 

251.20 

47 

SO.  1() 

25 

317.00 

48 

63.10 

26 

398. 10 

49 

79.44 

SO l.oo 
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Actual  filter  characteristics  are  no',  ideal  in  the 
sense  of  completely  rejecting  signal  components  outs.de  the 
nominal  passband.  The  effective  bandwidth  of  an  actual  fil¬ 
ter  such  as  the  1000  hertz  third  octave  filter  being  used 
in  the  system  is  about  21%  of  the  center  frequency  and 
the  filter  extends  from  891  to  1122  hertz.  This  is  ap¬ 
proximately  equivalent  in  determining  the  effective  band¬ 
width  of  a  third  octave  filter  by  injecting  a  white  noise 
signal  in  the  input.  The  white  noise  signal  is  uniform 
over  a  power-per-hertz  bandwidth  over  a  very  wide  fre¬ 
quency  range.  The  effective  bandwidth  is  then  the  total 
output  power  divided  by  the  output  power  for  a  one-hertz 
bandwidth  at  the  frequency  of  maximum  response. 

(2)  The  delta  function  is  an  array  of  numbers 
obtained  during  a  pre-process  calibration  run  (see  para¬ 
graph  4.2)  which  is  a  correction  factor  to  account  for  the 
frequency  response  of  a  given  hydrophone  within  the  dy¬ 
namic  range.  Table  10  shows  the  delta  function  array 
used  for  a  typical  ACODAC  deployment. 

5.6. 3.2  Hydrophone  Sensitivity 

The  hydrophone  sensitivity  value  used  in  the 
acoustic  equation,  equation  4.1-3,  and  entered  during 
initialization  was  obtained  from  data  derived  from  a  cali¬ 
bration  performed  by  the  Underwater  Sound  Reference  Divi¬ 
sion,  Naval  Research  Laboratory,  Orlando,  Florida.  liach 
hydrophone,  ITC  Model  8004,  used  in  the  ACODAC  system  was 
calibrated  at  temperatures  and  pressures  shown  in  table  11. 
The  hydrophone  sensitivity  value  used  in  the  initialization 
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TABLE  10 

DELTA  FUNCTION  ARRAY  FOR  SYSTEM  1A1  (November  1972) 


Third 

Octave  Hydrophone 

Band 


No. 

f/80 

#1 

#2 

#3 

#4 

#5 

#6 

5 1 

15.7 

-7.0 

-2.5 

-5.0 

-5.0 

-8.0 

-5.0 

32 

19.8 

-3.5 

0.0 

-2.0 

-3.0 

-5.0 

-2.0 

33 

24.9 

-2.5 

0.0 

-1.0 

-1.5 

-4.0 

-1.0 

34 

31.4 

-1.5 

.5 

0.0 

-2.0 

-2.0 

-.5 

35 

39.5 

-.5 

0.0 

.5 

0.0 

-1.0 

o 

o 

36 

49.8 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

37 

62.6 

0.0 

-.5 

0.0 

0.0 

0.0 

0.0 

38 

78.9 

0.0 

-1.0 

-.5 

-.5 

0.0 

-.5 

39 

99.3 

0.0 

-1.0 

-.5 

-.5 

0.0 

0.0 

40 

125.0 

-.5 

-.5 

-.5 

0.0 

0.0 

0.G 

41 

157.4 

-1.0 

0.0 

-.5 

0.0 

0.0 

0.0 

42 

198.1 

-3.5 

-.5 

-2.0 

.5 

-1.0 

1.0 

43 

249.4 

-5.5 

0.0 

-3.5 

0.0 

-2.0 

1.0 

44 

315.25 

-10.0 

-2.0 

-8.5 

-1.0 

-5.5 

-.5 

NOTE:  The  Delta  function  arrays  listed  above  for  hydrophones  1  through  6, 
system  1A1,  represent  a  system  calibration  during  a  run  conducted 
in  November  1972  and  is  only  shown  for  illustrative  purposes.  These 
values  will  vaiy  with  respect  to  adjustments  made  to  recorder/ 
playback  amplifier  settings  and  head  re-alignments . 
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TABLE  11 

FREE-FIELD  VOLTAGE  SENSITIVITY 
(Decibels  re  cne  volt  per  micropascal) 

(Sensitivity  values  negative)  -  ITC  Model  <*00*  Hydrophones  -  Unbalance  ' 
Voltage  across  470  ohm  resistor  at  output  of  amplifier 


Serial  1 


Serial  2 


Serial  3 


142.7 

143.4 

144.7 

341.4 

141.7 

142.5 

140.6 

140.0 

141.6 

138.9 

139.1 

139.2 

139.7 

138.6 

139.8 

137.8 

138.0 

137.6 

139.5 

138.1 

139.  1 

137.6 

137.7 

137.3 

139.5 

138.3 

138.7 

137.6 

137.6 

137.2 

139.5 

138.3 

138.4 

137.5 

137.6 

137.0 

139.5 

138.3 

138.4 

137. S 

137.6 

136.4 

139.5 

138.3 

138.4 

137.5 

137.6 

135.5 

.5 

.5 


136.8  137.5 


22°C 

3.5 

6J.95 

kPa 

k°a 

142.? 

142. 9 

139.2 

139.4 

137.7 

138.0 

137.4 

137.5 

137.2 

137.4 

137.1 

137.4 

137.0 

137.1 

137.0 

137.1 

Serial  5 


Serial  7 


Serial  8 


um 

142.8 

mm 

141.9 

142. 1 

142.7 

EES 

139.3 

IQS 

139.4 

139.6 

139.4 

138.0 

137.7 

138.3 

133.2 

137.8 

137.9 

137.3 

139.0 

137.8 

138.0 

137.5 

137.7 

137.2 

138.9 

137.8 

137.9 

137.5 

137.7 

137.1 

138.5 

137.8 

137.7 

137.5 

13V. 7 

137.1 

138.4 

137.7 

137.6 

137.2 

137.7 

137. 1 

138.4 

137.7 

137.6 

137.2 

142.5 

145.3 

138.9 

140.3 

137.7 

139.2 

137.6 

158.8 

137.5 

138.7 

137.5 

138.4 

137.1 

137.8 

137.1 

137.8 

Serial  9 


Serial  10 


42.8 
39.3 
38.1 
.9 


Serial  11 


141. 

8 

139. 

*7 

/ 

138. 

7 

138. 

5 

138. 

,4 

158. 

,4 

138. 

,4 

i38. 

.4 

142.2 

143.3 

140.7 

139.5 

139.7 

138.0 

138.8 

137.7 

138.6 

137.5 

138.6 

137.5 

138.4 

137.5 

138.4 

137.5 

(continued) 
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FREE-FIELD  VOLTAGE  SENSITIVITY 
(Decibels  re  one  volt  per  micropascal) 


Serial  4 

Serial  6 

3 

°C 

22° 

C 

3‘ 

’C 

22° 

C 

Freq 

3.5 

6895 

6895 

3.5 

6895 

ra 

6895 

(Hz) 

kPa 

kPa 

kPa  I 

kPa 

kPa 

kPa 

EBB 

kPa 

HSi 

141  8 

141.6 

142.1 

L41.8 

141.8 

TBKI 

HP 

139.4 

139.3 

139.0 

139.1 

i  Saw 

139.4 

HT~ 

138.8 

138.7 

137.8 

137.8 

[ 

138.3 

137.9 

■!  1 

138.  e 

137.6 

13.  .6 

Wm 

138. 1 

137.6 

137.6 

Wt&l 

138.2 

138.0 

137.6 

137.6 

!  H 

137.9 

137.6 

137.6 

137.  {j 

138.0 

137.6 

137.6 

'  laplsl 

137.9 

137.6 

137.6 

800 

137.8 

138.0 

137.6 

137.6 

'  'Mvl® 

137.9 

137.6 

137.6 

1000 

137.8 

138.0 

137.6 

137.6 

[137.9 

137.9 

137.6 

137.6 

ITC  Model  8004  Hydrophones 

Voltage  across  470  ohm  resistor  at  output  of  amplifier 

Unbalanced 


(Sensitivity  values  negative) 
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SAMPLE  PRINTOUT  OF  INITIALIZATION  DATA  ENTERED 
IN  1/3  OCTAVE  PROCESSING  PROGRAI^^  *  IZ.j'Ji, 
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0 

CC  ? 
ID 

0 

IDLE 

0 

CC  ? 
IN 


3£Ai-l 

f??o<.esse 

k"  . 

X4  J./ 


l)*J?3 


MT  UR :  (  DA  TA  ,  UA I  N)  : 

ENTER  0  FOR  NO  GAIN  TAPE  OUTPUT:  7,0 
LP  UR  #=  2 


CAL  DUMP  UR=? 

8 

JUMP  TO  LINE  CHANGE? 

NO 

1.  MOORING  TYPE, DEPLOYMENT  NOJ  2A4  001D 

2.  SHIP, CRUISE  NO:  NS  7*0 

3.  HYDRO  STA,EER  NU,  4  DtPlHJ  1  0UU1  Ob^ 

4.  HYDhuPhONc  SENSII1V1TY:  -t>0.0 

5.  PRE-CAL  THANSF.  FUNCT.  (TPhlM):  4:>.U 
CAL  TUL=1.0 

4  CAL  LEVELS  ARE: 


-41.0  -51.0  -21.0  -11.0 

ENTER  CHANGE  1-4  (X,XX> 

OH  0  FOR  NO  CHANGE 

0 

6.  DEPLOYMENT  LAUNCH  DAY  NO,  YEmr;  554,72 

7.  RECORDER  START  TIME:  (DAY  #,nH,M.N:  5-1-.,! \"J 

PROCESS  IN  RELATIVE  TIME  On  gHT 
RT=0  AND  GMT  =  l:  0 

8.  RECOVERY  DAY  NO.,  TIME (GMT ) :  54^,1552 

$.  MIN  4  MAX  SAND  DISPLAY:  5l,4I> 

10.  LOC  TIME  ZONE,  TIME  DRIFT  CONST: 

11.  DEPLOYMENT  LAT  (DeG,MIN>:  18, 4*.0 


12.  DEPLOYMENT  LosU  (DEU,MIN):  -79, 52.7 

13.  HYDRO  MFOR  <6  ALPHA):  ITC 
1  A .  HYDhu  MODEL  No.  (8  alPaa):  il'C~BOEO 

is.  playback  time  compression  ts  ua  bo.-d: 
lb.  analou  tape  id:  1272 

17.  PRL-C)  L  DELTA  *  UNO  flow 

ENTEK  ACTIVE  BmNdS  ONLY;  TEhMi.  NA  Ie  LILT  wiTrt  BhNu  N0=0 
HAND  No.  (XX),  VALUE  (XX. X  DB) 

0 

ib.  bandwidth  correction:  oo  bands  vs  5o.u  db  t«x 
ENTER  ACTIVE  BANDS  ONLY;  TEHM  Li  SI  wilh  uaNd  #si. 

xx,xx.x 

31.5.62 
32,5. 

3.’, a. 62 
3  >, 

3  j  ij .  62 

35.9.62 

36.10.62 
37, 1  1  .62 
38, 12  .62 

39.13.62 

40.14.62 

41.15.62 

42.16.62 

43.17.62 

44.18.62 
45,24.50 
0 

19.  PARAMETER  INSERTION  SECTION 
ENTER  MAX  15  VALUES. 

TerM  WITH  LINE=0 


PRINTOUT  CONTINUED 

TABLE  12 
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DeFaULT  VALUES  USED: 

3,15  -  F'i  N  FHAME  SYNC 
4,36  50  REF  BAND  NU 
14,  1  HYD  STA  NO 

0 

20.  OBSERVATION  ON-TIME  4  UFF-TIMe 
0,0  FUR  CONTINUOUS  (IN  MIN):  0,0 

DUTY  CYCLE: 100«—  CONTINUOUS 


21.  DlulTAL  TAPE  ID: (Da Ta  6  omIN):  XAeI 

22.  WANT  BRUSH  PLOT  OUTPUI? 

YES 


ENT-  SdLtClED  BAND  #  Kir  EaCH  CHaN 
aN d  Lower  PLof  LIMIT  (do);  4  baNdS  Max 

FORMaT  :  XX, YY  WHERE  X:oANd  No  4  Y=  RaNoe  LOWea  LxM(dd) 


ortUSH  CHAN: 

1 

BANS  #,LOV 

LI  Pi 

36,-3u,o 

PRINTOUT  CONTINUED 

SHUSH  CHAN: 

2 

BANS  #,LOw 

LI  Pi 

39,-30.0 

TABLE  12 

SHUSH  CHAN: 

3 

BANS  f,LQ» 

Lx  Pi 

38,-30.0 

1 

o 

00 

1 

SHUSH  CHANi 

4 

BANS  #,L0W 

LI  Pi 

46,-30.0 

SELECTIONS 

MALE: 

BANS  it  36 

lim: 

-30.0 

BAND  it  39 

LIPU 

-30.0 

BANS  it  38 

LIM: 

-30.0 

BAND  it  42 

LI  PL 

-30.0 

CAL  LEVEL  LOW  LIP! 


0.0 


ACCEPT  ALL  INPUT  VALUES? 
YES 

INPUT  DATA  COMPLETE 


0 

<9  0 
IULE 
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Phase  of  the  program  was  that  value  taken  at  a  temperature 
of  3°C,  at  a  pressure  of  689S  XPa,  and  at  a  frequency  of 
50  hertz  for  the  particular  hydrophone  being  used.  This 
value,  as  can  be  seen  from  the  table,  does  not  compensate 
fully  for  the  effects  of  actual  temperature  and  pressure 
of  the  hydrophone  at  a  specific  depth  and  for  the  frequency 
response  through  the  third  octave  filters. 

The  actual  hydrophone  sensitivity  value,  N. 
required  by  the  acoustic  equation  can  be  expressed  as: 

J  h  j  j  C5  b.3,2-x) 

where  j  =  tnird  octave  band  number 

Sh  *  hydrophone  sensitivity  obtained  from  the 
table 

=  thp  correction  applied  to  account  for  the 
frequency  response  of  a  given  hydrophone 
within  the  dynamic  range  for  band  j 

Kj  -  correction  factor  applied  to  compensate  for 

actual  temperature  and  pressure  of  a  hydro- 

pu.one  at  a  given  depth  and  at  the  center 

frequency  of  the  third  octave  band. 

The  acoustic  equation  accounts  for  the  parameters  SL  and 

h 

which  are  entered  in  the  program  during  initialization. 
However,  the  "K"  factor,  K.,  is  introduced  as  a  correction 
during  the  editing  phase  of  the  third  octave  analysis  pro¬ 
cess  rather  than  during  initialization.  This  *as  done  to 
expedite  third  octave  processing  because  of  the  computations 


involved. 
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5.6. 4  TEST  Mode 

The  TEST  mode  is  used  to  check  system  operation  prior  to 
entering  the  data  acquisition  Dnase.  The  main  difference  between 
this  and  the  ACQUISITION  mode  is  that  no  output  data  records  are 
writ  jn  on  the  digital  magnetic  tape. 

5.6.5  TERMINATION  Mode 

It RUINATION  is  used  to  complete  the  acquisition  phase  by 
terminating  the  digital  magnetic  tape  records  with  a  double  end- 
of-file.  Once  a  double  end-of-file  is  written,  the  program 
backs  over  the  last  end-of-file  in  preparation  for  processing 
another  hydrophone. 

5 . 7  Command  Keyboard  Entries 

The  operator  can  at  all  times  communicate  with  the  program  regardless 
of  the  program  mode  currently  in  execution.  The  operator  interacts  with 
the  program  by  initiating  keyboard  commands  via  the  teletype  or  CRT  dis¬ 
play  unit.  The  commands  allow  the  operator  to  change  modes  of  operation, 
request  information  to  be  printed  onto  hard  copy  and  dynamically  change  or 
read  out  specific  memory  locations.  Program  operation  is  not  interrupted 
Curing  this  process  •-ir.ee  the  entire  conversation  is  accomplished  in  back¬ 
ground.  A  list  of  keyboard  commands  is  shown  in  table  13. 

5 . 8  Program  Error  Diagnostics 

The  program  prints  out  eiror  diagnostics  once  every  minute,  80  minutes 
in  hydrophone  time.  These  messages  provide  dynamic  status  with  respect  to 
the  state  of  processing  during  acquisition.  The  messages  inform  the  opera¬ 
tor  when  unusual  conditions  exist  such  as  time  synchronisation  errors  oc¬ 
curring,  time  and  gain  code  errors  detected,  and  plot  limits  exceeded. 

The  messages  take  the  form  of  a  three  or  four  character  mnemonic  to  identify 
the  exact  error  diagnostic.  The  format  printed  is: 
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TAltl.l.  15 

COMMAND  khYISOARb  l-.NIHII.S 
Program  Adi)  1  -  Ambient  Noise  Processing 


Keyboard  Command 
*D@X 
*C@A,3A 
IN 
AQ 
UN 
All 
'll: 

ID 

TS 

RD 

CD 

UN 

CG 

SI) 

SY 

IC 

CA 

DC 

DU 

RC 

SB 

SP 

IT 

iC 


Descrioti on 
- .. 

Dump  memory  location  at  octal  location  X 

Change  memory  location  X  to  octal  value  A 

Start  initialization 

Start  acquisition 

Terminate  deployment 

Abort  deployment 

Start  test  mode 

Start  idle 

Resynchronize  24  hour  clock  to  time  code 
Raw  data  printout  -  Tl;  mode 
Corrected  data  printout  -  TL  mode 
Select  third  octave  band  for  printout 
Print  out  current  value  of  hydrophone  gair. 

Stop  test  data  printout 

Re-synch  days,  hours,  and  minutes  to  time  code 
Compare  computer  time  with  time  code  -eading 
Start  calibration  (special  test) 

Dump  current  IS  bands  of  corrected  data 

Dump  current  15  bands  of  ur.corrected  data 

Change  range  of  brusn  plot 

Start  brush  plotting 

Stop  brush  plotting 

Dump  calibration  data 

Print  integration  count 
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1220 

3 

src  2 

HRBD  5 


meaning  that  at  relative  time  1220,  day  3,  two  types  of  errors  were  detected: 
two  gain  code  errors  (STC)  and  five  hydrophone  amplifier  overloads  (ERBD) . 
These  errors  were  detected  within  the  time  period  of  80  minutes  since  the 
last  error  message  printout.  Table  14  shows  a  list  of  error  messages  cur¬ 
rently  being  used.  If  no  errors  occurred  during  the  80  minute  hydrophone 
period,  only  the  time  will  be  printed. 

S.9  Analog  Plots  of  Third  Octave  Bands 

Figure  13  shows  a  sample  Brush  analog  plot  of  ambient  noise  data.  The 
analog  plot  contains  6  channels  of  data  continuously  recorded  for  the  en¬ 
tire  period  of  the  hydrophone  through  the  third  octave  system.  The  analog 
channels  contain  the  following  information: 

Channels  1-4:  any  four  out  of  15  bands  selected  for  plotting.  The 

average  spectrum  levels,  in  10  second  averages,  are  in 
microbars.  Generally,  the  selected  third  octave  bands 
are: 


Channel  1:  50  Hz 

Channel  2:  100  Hz 

Channel  3:  75  Hz 

Channel  4:  200  Hz 

Channel.  S:  is  the  computer  output  gain  state  level  of  the  hydro¬ 

phone  being  processed.  Range  in  dc  level  represents 
data  amplifier  g3ins  of  10,  20,  30,  or  40  db. 

Channel  6:  Hydrophone  gain  as  read  directly  from  the  time  code 

reader.  This  redundancy  has  been  useful  to  detect 
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TAB  LL  14 

ERROR  MESSAGE  SUMMARY  FOR  PROGRAM  ACUI 


Error  Mnemonic 

Description 

MINI 

A  minute  interrupt  was  not  detected 
after  5  integrations 

I  ROD 

Day  code  error 

I  ROM 

Minute  code  error 

IROH 

Hour  code  error 

OVSI  i 

Integrations  not  in  sync  with  minute 
time  frame 

ERBP 

Hydrophone  overload  detected 

1  0 

Input  exceeding  output 

SIC 

Gain  code  error 

RSYC 

Synchronization  error 

UBX 

Brush  channel  X  upper  iimit  exceeded 
where  X  =  Biush  channel  1  thru  4 

LbX 

Brush  channel  >  lower  limit  exceeded 
where  X  =  Brush  channel  3  thru  4 

Figure  13.  Sample  Ambient  Noise  Brush  Analog  Plot  Output 
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Channel  7: 


Channel  8: 


incorrect  gain  state  decoding  via  the  computer  process¬ 
ing.  If  channels  5  and  6  do  not  agree,  problems  exist 
either  in  the  computer  processing  program  or  system 
hardware . 

The  uncorrected  signal  from  the  processing  hydrophone, 
through  a  50Hz  third  octave  filter,  is  plotted  directly 
on  the  Brush  plotter  (not  through  computer) .  Channel  7 
is  used  to  check  validity  of  data  being  processed  through 
the  computer.  ual  checks  are  made  between  the  uncor¬ 
rected  50  Hz  signal  level  of  channel  7  with  the  corrected 
50  Hz  signal  level  being  recorded  in  one  of  the  first 
4  channels  from  the  computer.  In  order  to  compare  the 
levels  of  the  two  SO  Hz  signals,  a  process  can  be  used 
to  correct  the  uncorrected  50  Hz  signals  of  channel  7 
by  comparing  it  with  the  signal  level  of  channel  1.  If 
discrepancies  occur,  then  the  third  octave  process  should 
be  discontinued  until  the  problem  is  resolved.  Incor¬ 
rect  parameters  inserted  in  the  computer  program  during 
initialization  are  indicative  of  discrepancies  between 
these  two  channels. 

This  channel  is  used  to  provide  a  suitable  time  base 
when  analyzing  the  Brush  plots.  An  IRIG-0  time  code 
is  used  from  the  time  code  reader  which  represents  a 
one-hour  repetition  rate.  Each  DC  level  excursion 
represents  one  minute  of  hydrophone  time,  and  the  code 
frame  repeats  every  hour. 


-7o- 


Referring  back  to  the  initialization  printout,  table  1?,  the  program  asks 
the  operator  to  input  scale  factors  for  each  of  five  channels  to  be  used 
for  plotting  data  on  the  Brush  recorder.  The  program  stores  the  lower  limit 
of  each  band  selected  and  automatically  sets  the  upper  range  of  the  Brush 
plot  50  db  higher  than  the  lower  limit.  An  example  of  bands  to  be  plotted 
are  as  follows: 


56,  -40. 


37,  -40. 
39,  -40. 


15,  -40. 


The  above  entry  indicates  that  the  operator  wishes  to  plot  third  octave 
band«  36,  37,  39,  and  45  all  using  -40  db  as  the  lower  plot  limit. 

Since  the  output  to  the  Brush  recorder  through  the  digital -to-analog 
converter  is  a  10  volt  dc  level  with  full  scale  being  represented  in  256 

bits,  the  scale  factor  applied  would  be: 

.  SO.jb  fdygamc  rangej  ,  b/ 

256  bits 

Then  scaled  data  point  XC,  represented  by  a  10  second  average  spectrum 
level  (BJjj  in  db  re:  microbar/Hz,  would  be  scaled  as  follow^: 


XC. 

i 


(BJ)j  -  LIMj 
SFAC 


(5.9-1) 


(5.9-2) 


where : 

XC 

i 


LIM 


1 


SFAC. 


-  the  scaled  data  point  to  be  plotted  on  channel  i, 
where  i  -  1  to  4 . 

-  average  spectrum  level  for  band  j,  where  j  -  31  to  45. 
=  lower  limit  of  a  50  db  dynamic  range  of  the  Brush 

recorder , 

=  Brush  recorder  scale  factor  applied  through  a  D-A 
output  converter  wnerc  10  volts  represents  full  scale. 
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The  program  has  the  ability  to  change  the  lower  limit  of  any  selected 
band  in  the  event  that  an  error  message  indicates  that  out  of  bound  values 
are  being  recorded.  This  change  is  done  dynamically  while  the  program  is  in 
process  of  plotting  data  points. 

5.9.1  determining  Validity  of  Ambient  Noise  Data  Derived  from  Brush 
Recordings 

The  analog  Brush  recorder  output  generated  from  the  third 
octave  analysis  program  contains  plots  of  selected  third  octave 
ambient  noise  data  referenced  to  microbar  levels.  Eight  brush 
v.annels  are  recorded:  six  channels  plotted  directly  from  the 
computer,  one  represents  the  gain  of  the  hydrophone  from  the 
gain  state  decoder,  one  from  the  IRIG-D  slow  code  (hour  code) 
from  the  time  code  reader,  and  the  eighth  from  the  recorded 
broadband  signal  through  a  third  octave  50  Hz  filter  to  the 
brush  recorder. 

A  selection  of  1/3  octave  bands  for  a  particular  hydrophone 
is  shown  below.  This  is  typical  output  plot  sequence  on  a 
Brush  recorder  of  the  amb’ent  noise  data  from  an  ACODAC  deploy¬ 
ment. 


Brush 

Channel 

Signal 

Origination 

Band  Number  (80:1) 

1 

50  Hz 

computer 

36 

2 

100  Hz 

computer 

39 

5 

75  Hz 

computer 

37 

4 

200  Hz 

computer 

42 

5 

Hydro  gain 

computer 

6 

Gain 

gain  state  decoder 

•9 

50  Hz 

broadband  through  50  Hz  filter 

8 

time  code 

IRIG-D  slow  code 
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Note  that  two  amplifier  gain  state  oh:.. i. 'els  are  displayed. 
This  redundancy  has  been  used  as  a  cross-check  on  the  hydrophone 
amplifier  gain  setting  on  playback  to  insure  that  the  computer 
program  is  detecting  the  gain  states  properly  with  respect  to 
the  gain  state  being  decoded  from  the  time  code  envelope. 

One  purpose  of  the  Brush  recorder  output  data  is  to  assist 
in  determining  the  validity  of  data  as  the  data  is  being  digi¬ 
tized,  integrated  by  10  second  averages,  corrected  to  microbar 
levels,  and  written  onto  digital  magnetic  tape  for  further 
processing. 

The  following  procedure  should  be  used  to  verify  that  all 
system  constants  have  been  inserted  during  the  initialization 
phase  of  program  execution,  and  that  all  hardware  connections, 
including  proper  hydrophone  selections,  have  been  made. 

Calibration  Phase 

Figure  14  represents  a  Brush  recording  of  a  typical  cali¬ 
bration  sequence  which  was  recorded  during  the  deployment  phase 
of  ACODAC  once  every  six  hours.  During  the  calibration  phase,  a 
50  Hz  and  200  Hz  signal  is  injected  at  the  input  to  the  hydro¬ 
phone  data  amplifiei  with  the  hydrophone  input  disconnected. 

Four  levels  are  gonerated  for  each  signal  in  the  capsule,  each 
lasting  one  minute.  These  levels  are:  -41,  -31,  -21,  -11  db. 

The  levels  are  selected  to  force  the  amplifier  gain  states  to 
switch  to  the  four  dynamic  ranges  of  the  amplifier,  thus  allowing 
a  test  of  the  amplifier  gain  ranges  while  at  the  same  time  pro¬ 
viding  a  calibration  level  for  ambient  noise  analysis. 
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CALIBRATION  SEQUENCE -BRUSH  PLOTTER 


j-«* calibration  cycle*-|  tj 


FIGURE  14 

The  four  calibration  levels  with  the  corresponding  gain 


ranges  are  shown  below. 

Signal  level  (db) 

Gain  State 

-41 

40 

-31 

30 

-21 

20 

-11 

10 

Referring  to  figure  14,  observe  the  signal  through  a  50  Hz  third 
octave  filter.  If  the  system  is  behaving  normally,  the  resultant 
level  that  will  be  recorded  on  the  Brush  recorder  should  repre¬ 
sent  a  relative  level  of  (-1  db)  recorded  for  the  four  minute 
period.  The  procedure  requires  converting  this  raw  50  Hz  signal 
to  a  sound  pressure  level  -e:  1  »icrobar/Hz,  and  then  to  compare 
this  level  with  the  output  of  the  50  Hz  third  octave  signal  being 
plotted  from  the  computer  re:  1  microbar/Hz.  If  the  ambient  noise 
levels  of  the  two  signals  agree,  within  1  or  2  db,  then  it  can  be 
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assumed  that  the  computer  process  is  normal  and  that  all  initiali¬ 
zation  constants  and  hardware  interconnections  have  been  made 
properly. 

To  Convert  the  Raw  50  Hz  Signal  to  a  Sound  Pressure  Level  Re:MBar/Hz 

1.  Observe  the  calibration  level  for  the  four  minute  period.  The 
level  s'-.^uld  be  smooth  during  this  time.  If  uneven  levels  appear 
for  the  1  minute  period,  it  may  be  difficult  to  pick  the  correct 
level,  sirce  the  uneveness  was  probably  caused  by  the  amplifier 
not  switching  to  the  proper  gain  state  when  the  signal  level  was 
switched.  However,  if  the  correct  level  can  be  observed,  then 
this  level  represents  -1  db  re:  1  volt. 

2.  To  convert  the  calibration  level  to  microbars/Hz,  subtract 

the  hydrophone  sensitivity  (S^)  from  this  level  (CL):  -  S^. 

This  is  the  equivalence  of  generating  an  RMS  sound  pressure 
level  of  one  microbar  equal  to  1  dyne  per  square  centimeter. 

In  this  example,  the  hydrophone  sensitivity  is  -38.7  db/mbar. 

The  converted  level  then  becomes: 

-1  -  (-38.7)  =  37.7  db 

3.  Since  the  broad  band  s;gnal  was  passed  through  a  50  Hz 
third  octave  filter,  the  signal  must  be  corrected  for  band¬ 
width.  The  bandwidth  correction  for  the  50  Hz  signal  is  10.62 
db.  The  reference  level  then  becomes: 

37.7  -  10.62  =  27.08  db 

If  Xj  represents  the  ambient  noise  level  at  time  t^.  after  the 
calibration  cycle  has  been  completed,  then  the  true  ambient 
noise  level  re:  mbar/cycle  is: 

27.08  -  X. 
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The  level  X  therefore  can  be  used  for  the  next  six  hour  interval 
c 

to  determine  the  ambinet  noise  level  at  any  time  t. 

This  53  Hz  level  can  then  be  used  to  compare  as  a  reference 
the  50  Hz  signal  obtained  through  the  computer  process.  If  the 
two  levels  do  not  agree  within  one  or  two  db,  then  it  is  as¬ 
sumed  that  certain  initialization  parameters  have  not  been  in¬ 
serted  properly  in  the  computer  program  or  that  instrumentation 
malfunctions  may  exist. 

6.  ANALYSIS  OF  CONTINUOUS  WAVE  SIGNALS  USING  THIRD  OCTAVE  FILTERS  AND  1  HZ 
NOTCHED  FILTERS  TO  OBTAIN  TRANSMISSION  LOSS  DATA 
6.1  Input  Data 

The  continuous  wave  (CW)  or  narrow  band  analysis  is  similar  to  the 
third  octave  ambient  noise  analysis  with  the  following  exceptions: 

(11  An  array  of  six  1  Hz  wide  combed  filters  were  added  to  bandpass 
source  frequencies  of  20  and  130  Hz.  (Thc'e  frequencies  were  from  a  1972 
field  operation;  other  frequencies  are  possible  with  appropriate  filters). 

A  patch  panel  was  used  to  ar 'ange  th?  desired  signals  into  15  bands  which 
were  then  passed  through  the  General  Radio  1926  RMS  Detector.  From  then  on, 
the  process  wa-.  similar  in  nature  to  the  third  octave  system.  Figure  15 
shows  the  hardware  interconnections  of  the  analysis  setup.  For  certain 
applications,  filters  of  less  thac  one  hertz  bandwidth  may  be  used. 

(2)  The  patch  panel  was  used  to  connect  the  combed  filter  and  third  octave 
filter  outputs  to  selected  band  inputs  to  the  General  Radio  Real  Time  Ana¬ 
lyzer  Although  in  practice,  bands  31  to  45  represent  the  true  third  octave 
bands  of  the  frequency  inputs,  in  the  narrow  band  system,  bands  31  to  45  have 
no  meaning  to  the  band  frequency.  The  connections  through  the  patch  panel  must 
be  known  to  ascertain  which  band  is  associated  with  a  particular  frequency. 
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6 • 2  Calibrating  the  Narrow  Band  System 

The  transfer  function,  T'  ,  used  in  equation  4.1-1  represents  the  trans¬ 
fer  function  of  the  third  octave  analysis  without  the  narrow  band  filters 
connected.  The  true  system  transfer  function  of  the  narrow  band  system  was 
obtained  by  injecting  a  known  signal  level  into  the  system  of  figure  15 
at  the  center  frequency  of  the  20  and  130  hertz  filters.  The  correction  to 
the  overall  transfer  function  was  made  in  the  correction  equation,  4.2-4, 
by  using  the  BWC  array  since  the  bandwidth  correction  for  narrow  band  fil¬ 
ters  was  not  used.  Figuie  16  shows  a  more  detailed  flow  of  filters  used 
for  the  130  hertz  analysis. 

6 . 3  Computing  Transmission  Loss 

Once  the  c.w  energy  is  digitized  through  the  GR  Analyzer,  corrected  to 
microbar  levels  and  recorded  on  digital  tape,  the  digital  tape  is  sub¬ 
mitted  for  the  editing  process  to  eliminate  spurious  signals.  During  this 
process,  a  transmission  loss  table  is  generated  for  the  cw  source  frequency 
selected.  Table  15  shows  a  typical  set  of  transmission  loss  values  compu¬ 
ted  for  the  13C  Hz  source. 

6.3.1  System  Equations 

Consider  a  set  of  simultaneous  equations  which  are  used 
to  define  a  cw  source  system  through  narrow  band  and  third 
octave  filters.  The  average  power  of  the  total  energy  received 
and  recorded  during  the  time  the  cw  source  is  transmitted  is: 

Narrow  band  filter  system: 

The  total  average  power.  A,  obtained  through  the  narrow 
band  system  is: 

PA 

g(BW) 


P130S  ♦ 


(6.3-1) 


Genera! 


Figure  15.  Acodac  System  Diagram  for  CW  Analysis 


Hgure  16.  Acodac  Continuous  Wave  Analysis 
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where  P130S  *  average  power  of  130  Hz  source  through  1  Hz 

wide  narrow  band  filter  system 
PA  «  average  power  of  ambient  noise  energy  through 
narrow  band  filter  system 
g  =  weighting  factor  to  compensate  for  the  dif¬ 
ference  in  the  average  power  between  the 
center  frequency  of  the  third  octave  band 
(125  Hz)  and  the  signal  source  (130  Hz); 
g  =  1.22 

Third  Octave  filter  system 

The  total  average  power,  B,  obtained  through  the  third 
octave  125  Hz  band  is: 


P130S  PA 
BW  BW 


(6.3-2) 


where  P130S  =  Average  power  of  130  Hz  energy  through 

third  octave  filter  system 
PA  =  Average  power  of  ambient  noise  energy  through 
third  octave  filter  system  at  a  center 
frequency  of  125  Hz. 

BW  =  Bandwidth  of  125  Hz  filter;  BW  =  28.8  db. 
Using  the  simultaneous  equation  we  solve  for  the  ambient  noise 
energy,  PA: 


P130S  + 


PA 

g(BW)  " 


A 


(6.3-3) 


P130S  ♦  PA  =  (BW)B 


(6.3-4) 


solving  for  PA 


pa  s  JkifflQJktffP 

l-g(BW) 


(6.3-5) 
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The  average  power,  PA,  computed  for  ambient  noise  energy  in 
equation  6.3-5  must  be  corrected  for  the  effective  bandwidth 
for  the  actual  filter  in  order  to  obtain  the  power-per-hertz 
value  of  ambient  noise. 


PA 

Btf 


l-g(BW) 


Re-arranging  the  terms  yields: 


(6.3-6) 


Equation  6.3-7  and  6.3-8  can  now  be  used  to  obtain  the  equivalent 
of  source  plus  ambient  noise  from  the  total  power  received  dur¬ 
ing  the  period  when  the  source  is  transmitting. 

Substitutions  are  made  in  equation  6.3-8  to  obtain  the  power  re¬ 
ceived  through  the  129  Hz  narrow  band  filter  since  the  results 
showed  that  the  peak  source  energy  varied  between  the  129  and 
130  Hz  frequency  bands. 

Figure  17  shows  a  typical  Brush  plot  of  narrow  band  spectra 
of  the  130  Hz  recorded  by  hydrophone  no.  3.  Channels  1  and  2 
are  outputs  from  the  narrow  band  filters  while  channels  3  and  4 
represent  outputs  of  the  wide  band  signal  through  the  GR  1925 
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Channels  1  through  4  represent  the  ambient  noise  plus  cw 
energy  processed  through  the  computer  and  corrected  to  db  re: 
microbar/Hz.  Channel  S  is  a  plot  of  the  gain  of  the  hydrophone 
as  detected  via  the  cr tputer,  and  channel  6  is  the  same  gain 
levels  plotted  from  the  output  of  the  gain  state  decoder  and 
wired  directly  to  the  Brush  plotter.  Channel  7  represents  the 
uncorrected  wide  band  signal  through  a  separate  50  Hz  third  oc¬ 
tave  filter  connected  directly  to  the  Brush  plotter.  Channel  8 
is  a  plot  of  the  IRIG-D  slow  code  used  as  a  time  reference 
track . 

7.  EDITING  PROCESS  FOR  THIRD  OCTAVE  AMBIENT  NOISE  TAPES 

7. 1  Purpose 

The  editing  process  for  ambient  noise  data  is  required  to  eliminate  all 
spurious  signals  such  as  shots,  instrument  injected  noise,  and  improper  cor¬ 
rections  made  in  gain  ranging.  Figure  18  shows  a  flow  diagram  of  the  edit¬ 
ing  process  using  the  Sigma  7  as  the  processing  computer.  The  unedited  am¬ 
bient  noise  tape,  from  the  HP  third  octave  analysis  system,  is  mounted  as 
an  input  tape  to  the  system.  The  editing  program  provides  various  process¬ 
ing  options  which  are  selectable  at  run  time  by  input  cards  to  the  program. 

7.2  Editing  Options 

Editing  options  are  provided  as  card  inputs  to  the  program.  These 
options  are: 

Time  Window  Selection 

An  option  is  allowed  to  select  any  time  window  for  editing.  Either  a 
series  of  time  windows  can  be  used,  or  the  entire  tape  can  be  selected  for 
editing. 


Band  Level  Correction 


Corrections  can  he  made  to  any  or  all  third  octave  bands.  Corrections 
ure  made  in  dh  and  are  applied  as  algebraic  sums.  The  "kM  factors  discussed 
in  paragraph  5.6.3  2  are  typical  of  the  corrections  applied  to  individual 
bands  in  this  category 

Transfer  function  Corrections 

This  option  allows  for  a  transfer  function  correction  of  any  six  hour 
period.  If  selected,  the  editing  program  will  pick  up  the  new  transfer 
function  from  the  data  header  record  (word  5  and  6  of  table  7)  and  apply 
the  correction  to  the  ambient  noise  spectrum  levels  for  each  10  second  aver¬ 
age  and  for  each  band.  The  correction  applied  is  actually  the  difference 
between  the  new  transfer  function  and  the  one  used  during  the  processing. 

7.3  Non-Optional  Editing 

In  addition  to  the  selectable  options,  the  editing  program  will  elimin¬ 
ate  all  calibration  data  records  from  the  editing  tape  and  reject  all  10 
second  integration  intervals  that  have  been  tagged  with  a  99.9. 

7.4  Data  Validity  Test  Criteria 

Several  tests  are  made  to  validate  the  data.  Tests  are  applied  to  in¬ 
dividual  bands  to  check  for  spikes  in  data  which  may  be  induced  by  instrumen¬ 
tation  Checks  are  also  made  to  detect  shot  arrivals  and  erroneous  gain 
ranging. 

7.4.1  Band  Level  Test 

Each  band  is  checked  to  see  if  it  exceeds  a  certain  ratio 
of  the  wide  band  signal.  The  ratio  is  set  at  run  time  and  should 
be  a  number  between  .5  and  .85  of  the  wide  band  signal.  The 
default  value  is  .5.  Since  we  are  dealing  with  ambient  noise 
levels  rather  than  pure  cw,  it  can  be  assumed  that  the  third 
octave  energy  levels  are  distributed  throughout  the  IS  bands 


wis^ywwi'iw,!)- 


SI-  .  '"f 


-92- 


with  peak  levels  occurring  at  approximately  50  Hz  and  dropping 
off  at  the  higher  and  lower  frequencies.  It  should  be  remembered 
that  the  band  level  test  is  used  to  eliminate  spurious  signals 
caused  by  instrument  electronics  and  the  like.  Further  tests 
are  employed  to  check  for  shot  energy  returns  and  improper  gain 
ranging. 

7.4.2  Erroneous  Gain  Level  Shifting  and  Shot  Energy  Detection 

A  test  is  made  to  determine  whether  the  data  has  undergone 
a  shot  arrival  transition  or  an  erroneous  gain  level  shift. 

Since  the  purpose  of  the  program  is  to  eliminate  all  but  pure 
ambient  noise  levels,  these  phenomena  are  edited  out  of  the 
data. 


An  erroneous  gain  level  shift  is  defined  as  the  condition 
existing  when  the  third  octave  processing  program  introduces 
a  10  db  gain  step  in  the  average  spectrum  levels  for  each  of 
15  bands  during  the  process  of  converting  raw  ambient  noise  to 
microbar  levels.  This  condition  occurs  when  the  program  gets 
out  of  synchronization  with  the  IRIG-C  time  code  frame  thus 
causing  a  step  function  to  appear  in  the  data.  More  precisely, 
a  gain  adjustment  is  being  applied  to  data  in  the  wrong  minute 
time  frame. 

The  editing  program  first  forms  the  difference  between 
the  sum  of  the  first  ten  bands  of  the  ith  and  ith  ♦  1 
10  second  interval.  The  sum,  S,  of  the  ith  interval  is: 

10 


(7.4-1) 
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where  *  the  average  spectrum  level  of  the  ith  10 

second  integration  interval  in  db. 
j  *  the  third  octave  band  number  representing 
the  frequency  range  of  15  to  125  Hz. 

The  sum  of  the  ith  +  1  interval  ir: 

ID 

Si-1  *  l  (Bi+1)i  (7*4“2) 

j-1  J 

The  difference  is  then 

AS  *  Si+1  -  Si  (7.4-3) 

This  difference  is  then  tested  to  see  if  it  exceeds  a  level  of 
26  db. 

It  can  be  assumed  that  for  prediction  purposes,  average 
representative  ambient  noise  spectra  will  not  usually  exceed 
26  db  level  shifts  in  short  time  intervals  with  respect  to  fre¬ 
quency  ranges  of  10  to  300  Hz,  and  the  test  appears  to  be  a 
valid  threshold,  although  selected  rather  heuristically,  when 
attempting  to  separate  ambient  noise  from  shot  arrivals. 

Nevertheless,  from  the  data  processing  results  obtained, 
the  26  db  threshold  appears  to  work  in  isolating  ambient 
noise  from  shot  arrivals.  Since  most  of  the  shot  energy  appears 
in  the  frequency  range  of  50  to  100  Hz,  the  sum  of  the  first 
ten  third  octave  frequency  bands,  between  15  and  125  Hz,  were 
selected  as  the  frequency  ensemble 

Now,  if  the  difference  of  the  suns  represented  by  AS 
exceeds  the  26  db  level,  further  tests  are  made  to  isolate  be¬ 
tween  a  true  shot  arrival  and  an  erroneous  gain  level  shift. 

This  isolation  can  roadily  be  accomplished  by  testing  AS  to  see 


if  it  exceeds  90  dh.  Recall  that  an  erroneous  gain  level  shift 
will  introduce  a  10  db  spectrum  level  error  in  each  band.  If 
the  absolute  value  of  the  difference  between  the  ith  and  the 
ith  ♦  I  spectrum  level  intervals  exceed  90  db,  then  the  ith  ♦  1 
interval  contains  an  erroneous  gain  level  shift,  or 


XU  XU 

45  -  I »“  -  .1 0»i>j 

J=1  1  J=1  J 


(7.4-4) 


AS  >  90  db 


If  the  10  second  interval  contains  a  shot  arrival,  then 
26  <  AS  <  90  db 

If  any  of  The  five  10  second  intervals  fails  to  pass  the  data 


(7.4-5) 


(7.4-6) 


validity  test,  the  minute  interval  is  rejected  from  further  pro¬ 
cessing  and  is  omitted  from  the  edited  tape  being  generated  and 
from  the  computation  of  the  mean. 

7.5  Data  Averaging 

A  maximum  of  15  data  means  can  be  computed  while  undergoing  the  editing 
process.  These  means  are  obtained  by  converting  the  10  second  average  spec¬ 
trum  levels  to  power,  forming  the  sum  and  then  converting  back  to  decibels. 
These  means  are  printed  on  the  high  speed  printer  for  all  15  bands  along 
with  the  time  of  d*y  and  the  mean  time  base,  i.e.,  2  minute,  10  minute, 

60  minute  means,  etc. 

7.6  Error  S'jcmary  Printout 

A  hardcopy  printout  is  provided  showing  a  suranary  of  all  shot  intervals 
and  erroneous  gain  level  shifts  encountered  during  the  editing  process. 

See  figure  19  for  sample  printout.  This  summary  is  helpful  in  listing  the 
minute  intervals  omitted  during  the  editing  process  and  provides  an  accurate 
record  of  where  the  shot  sequences  were  detected  during  the  deployment  phase. 
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Figure  19.  Error  Summary  Printout 


7.7  Processing  from  the  Edited  Tape 

Once  the  edited  or  clean  tape  is  produced,  further  processing  is  re¬ 
quired  to  obtain  meaningful  !:atistical  data  and  plots  of  the  average  spec¬ 
trum  levels  foi  the  ambient  noise.  These  data,  in  the  form  of  12  hour 
statistics,  frequency  histograms,  and  plots  provide  the  final  output  stages 
of  data  processing  presently  available  for  ambienc  noise  representation. 


8.  STATISTICAL  SUMMAP.Y 

The  program  ACOSTAT  is  used  to  compute  statistics  on  ambient  noise  data  in 
the  form  of  10  second  average  spectrum  levels.  The  program  uses  the  10  second 
interval  as  the  data  base,  and  all  statistical  functions  are  computed  in  decibels. 

The  statistics  shown  in  figure  21  are  .  ^dard  moment  statistics  which  can 
be  found  in  any  statistics  text,  e.g..  Statistical  and  Experimental  Design  by 
N.L.  Johnson  and  F.C.  Leone.  The  chi-squared  test  (last  two  lines,  figure  21) 
is  a  test  for  normalcy  of  the  decibel  distribution  function.  Figure  20  shows  the 
hydrophone  station  identification  for  the  statistics  being  computed. 

Program  ACOSTAT  also  has  an  option  to  output  a  digital  plot  of  ambient  noise 
spectrum  levels  for  selected  frequency  bands  vs  time  with  each  data  point  plotted 
representing  a  seven  minute  mean.  This  plot  is  not  widely  used  and  has  been  super¬ 
seded  by  the  one  minute  mean  plot  to  be  described  below. 

If  selected,  a  frequency  histogram,  figure  22,  of  each  third  octave  band  is 
printed  on  the  high  speed  printer  giving  the  class  intervals  and  population  density 
of  the  distribution.  See  Appendix  B  for  a  program  description  of  ACOSTAT. 

9 .  ACODAC  AMBIENT  NOISt:  PLOT 

Tli is  program  is  used  to  plot  ambient  noise  data  using  the  10  second  intervals 
as  the  data  base  and  computing  a  one  minute  mean  for  each  data  point.  All  10  sec¬ 
ond  data  points  are  converted  to  power,  the  mean  is  computed  and  then  converted 
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Figure  20.  Station  Identification 


Figure  21.  Statistical  Summary  of  Ambient  Noise  for 
Selected  Third  Octave  Bands 
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FREQUENCY  HISTOGRAM  FOR  A  THIRD  OCTAVE  BAND 
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Figure  22.  Frequency  Histogram  for  a  Third  Octave  Band 
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back  to  decibels  for  plotting.  Figure  23  shows  a  representative  plot  of  three 
frequency  bands  plotted  against  time.  The  data  points  are  converted  from  micro¬ 
bars  to  micropascals  for  plotting.  The  number  of  points  plotted  is  shown  for 
each  band  along  with  a  histogram  of  the  distribution. 


PROGRAM  SPECIFICATIONS  FOR  THE  NARROW  BAND  ANALYSIS  PROGRAM  (SIGMA  7) 
(ANALOG-TO-DIGITAL  CONVERSION  PROGRAM) 

10.1  Program  Specifications 

This  program  has  the  ability  to  input  raw  ambient  noise  data  from  a 
Sangamo  7-channel  recorder/play-back  unit.  Each  hydrophone  channel  (6 
total)  are  connected  directly  to  analog  input  channels  of  the  MD-51  digi¬ 
tizer  in  sequential  order. 

The  program  automatically  starts  and  stops  the  digitizing  process 
under  computer  control.  The  start  time  is  deteimined  by  decoding  time  in¬ 
formation  from  the  Time  Code  reader  in  days,  hours,  minutes.  The  time  code 
reader  is  interfaced  to  the  computer  via  the  A/D  controller  unit.  The  pro¬ 
gram  starts  by  sending  the  appropriate  control  word  out  to  the  A/D  controller 
to  start  the  time  reading  process.  A  series  of  punched  cards  are  an  input 
to  the  computer  which  have  the  start  times  and  the  duration  of  each  digi¬ 
tizing  cycle.  Once  coincidence  with  the  start  time  is  made,  the  program 
starts  the  digitizing  cycle  for  the  duration  designated  on  the  card.  This 
process  continues  until  the  entire  tape  has  been  digitized  or  until  speci¬ 
fied  series  of  cycles  has  been  accomplished.  Figure  24  shows  the  hardware 


interconnections  to  the  Sigma  7  via  the  high  speed  multiplexer  digitizer. 

An  initialization  parameter  is  used  to  provide  the  program  with  the 
program  with  the  duty  cycle  being  used  on  the  ACODAC  tape.  Two  duty  cycles 
are  expected:  50%  (7  minutes  on  and  7  minutes  off)  and  100%  continuous. 


Figure  24.  Sigma  7  -  Narrow  Band  Analysis  System  for 
Ambient  Noise 
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The  program  has  the  ability  to  set  the  full  scale  gain  of  the  digitizer 
to  obtain  maximum  resolution  of  the  digitized  data  point.  Resolution  of  the 
MD-51  is  14  bits  plus  sign.  Automatic  gain  can  be  provided  by  the  computer 
allowing  for  full  scale  resolution  of  input  voltages  of  1.2S,  2.5,  5.0,  and 
10  volts.  The  input  voltage  level  from  the  tape  playback  will  normally  not 
exceed  a  maximum  of  2.5  volts  full  scale.  This  -represents  a  gain  of  4:1. 

The  computer  program  sets  the  A/D  converter  to  digitize  6  channels  in 
the  sequential  mode  with  clock,  meaning  that  the  zero  axis  crossing  detec¬ 
tor  triggers  the  convert  cycle  at  the  sampling  rate  set  by  the  time  code 
reader  signal  connected  to  the  detector. 

The  computer  control  word  sets  the  digitizer  to  the  double  word  trans¬ 
fer  mode  such  that  two  15-bit  data  points  are  accumulated  in  the  controller 
before  the  32  bit  word  is  transferred  to  the  computer.  This  method  will 
effectively  reduce  the  A/D-to-computer  throughput  rate  by  one  half. 

The  program  will  detect  one  full  scan  of  six  hydrophone  channels 
(A/D  channels  0  through  5)  by  testing  the  32  bit  transfer  for  bits  15  or 
31  set  to  "1".  Since  two  dcta  points  are  packed  per  32  bit  word,  the  first 
conversion  will  always  be  in  the  right  half  of  the  32  bit  word  (bits  16-31) 
while  fci  ;  second  conversion  will  show  up  in  the  left  half  of  the  32  bit  word 
(bits  0-15).  Since  an  even  number  of  channels  are  being  digitized  (6  hydro¬ 
phone  channels  total),  bit  15  of  the  left  half  when  set  to  "1"  will  indicate 
that  the  right  half  of  the  next  32  bit  transfer  will  be  the  start  of  a  new 
scan,  or  channel  0.  The  output  format  for  a  dual  halfword  transfer  is  shown 
on  the  following  page. 


LEFT  HALF 


RIGHT  HALF 


0 

L 

S 

B 

X1 

S 

L 

S 

B 

0 

0 

14 

15 

16 

30 

31 

S  *  SIGN  BIT 

Xj,  X2  When  bit  is  set  indicates  that  one  complete  scan  has  elapsed. 
10.2  Digitizer  Sampling  Rate 

The  sampling  rate  of  the  digitizer  is  controlled  by  the  Zero  Axis 
Crossing  Detector.  This  box  is  connected  between  tnc  output  of  the  time 
code  signal  recorded  on  the  ACODAC  tape  and  the  "sampling  rate  .input"  to 
the  digitizer. 

A  294.12  Hz  signal  is  modulated  on  the  time  code  signal.  This  frequency 
is  the  maximum  frequency  of  interest  (in  hydrophone  time)  for  analysis.  In 
sampling  a  time  function,  this  implies  that  the  system  must  be  able  to  sample 
slightly  more  than  twice  per  period  of  the  highest  frequency  one  wishes  to 
resolve.  This  makes  reference  to  Shannon's  sampling  theorem  where 


F  < 
max 


1 

2At 


(10.2-1) 


where:  At  =  sampling  period 

F  =  maximum  frequency  of  interest 

The  zero  axis  crossing  detector  is  designed  to  generate  the  Nyquist  fre¬ 
quency  of  588.24  Hz  and  to  multiply  this  signal  by  six,  the  number  of  hydro¬ 
phone  channels  being  sampled. 

A  playback  rate  of  10  to  1  speedup  is  used  to  digitize  the  hydrophone 
data.  Therefore,  the  overall  sampling  rate  as  seen  by  the  digitizer  is: 


or 


where 
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Sampling  hate  =2  f  X  N£  x  S 


=  2  (294. 12) (6) (10) 


(10.2-2) 


Sampling  Rate  =  35,294.4  samples/sec. 

Fjaax  =  Maximum  frequency  of  interest  (hydrophone  time) 
N£  =  Number  of  hydrophone  channels 
S  =  Playback  rate 


10.3  Frequency  Resolution 

Frequency  resolution  (Af)  is  defined  as  the  reciprocal  of  the  total 
record  length,  T,  or 

Af  =  i  (10.2-3) 

For  ACODAC  Analysis,  a  one  Hertz  resolution  has  been  desirable.  To  obtain 
this  resolution,  the  length  of  time  required  to  digitize  one  continuous 
analog  data  section  is  one  second  in  hydrophone  time.  Since  a  playback  rate 
of  10:1  is  employed,  this  time  is  reduced  t<  0.1  seconds  for  each  data  sec¬ 
tion. 


The  maximum  number  of  data  points,  N,  required  to  obtain  a  one  hertz 
frequency  resolution  based  upon  the  maximum  frequency  of  interest  of  294 
hertz  (hydrophone  time)  is  therefore 

N  =  *  583  (10.2-4) 

10.4  Analog  Filtering  and  Aliasing 

ACODAC  ambient  noise  data  is  digitized  so  as  to  obtain  a  spectral 
resolution  of  one  hertz  bands  to  a  maximum  frequency  of  300  Hz  in  hydro¬ 
phone  time.  To  avoid  aliasing  below  the  maximum  frequency  selection  point 
requires  the  use  of  a  low  pass  analog  filter  and  a  deteimination  of  the  fold¬ 
ing  frequency,  a  factor  characteristic  of  the  type  of  filter  used. 

Aliasing  is  a  phenomenon,  inherent  in  analog  input  devices  which  results 
in  higher  frequencies  being  folded  back  onto  lower  frequencies,  making  the 


f 
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user  believe  that  there  are  frequencies  present  below  the  maximum  frequency 
selected  which  in  fact  may  not  be  there. 

To  circumvent  this  problem,  analog  filters  such  as  the  Allison  type 
can  be  used  between  the  data  output  channels  of  the  Sangamo  playback  unit 
and  the  input  to  the  analog  channels  of  the  digitizer.  The  cutoff  frequency 
of  the  filter  is  set  to  have  a  flat  response  to  the  maximum  frequency  of 
interest,  in  this  case,  300  Kz,  and  to  roll  off  beyond  this  frequency  at  some 
rate  of  attenuation  to  prevent  "folding  back"  of  the  higher  frequencies. 

Since  analog  filters  do  not  have  a  sharp  drop  off,  or  attenuation,  at 
the  cutoff  frequency  selected,  frequency  and  amplitude  information  around 
the  300  Hz  band  will  be  somewhat  distorted  and  misleading.  To  prevent  this, 
the  filter  cutoff  point  is  taken  higher  than  the  expected  maximum  frequency 
of  interest. 

10. S  Initializing  Parameters 

The  program  uses  the  following  parameters  entered  via  data  input  cards 
to  initialize  the  digitizing  program.  Th  jse  parameters  are  used  to  set  the 
proper  control  functions  of  the  digitizer  as  well  as  to  label  the  digitized 
data  which  will  be  written  on  the  output  magnetic  tape  m  the  first  file. 

The  following  parameters  will  be  used  to  set  the  control  functions  of 
the  digitizer  and  time  code  reader. 

1.  Maximum  full  scale  gain  setting:  2.5  volts  full  scale. 

2.  Type  of  scan:  sequential  with  clock 

3.  Digitizer  Control  Indicator:  1  -  start  the  digitizing  process 

on  time  code 

0  -  start  digitizing  process  on 
manual  toggle. 

4.  Data  points  per  transfer:  single/double 

5.  Digitizing  rate:  rate  per  channel  (hydrophone  time) 

l 

% 

I 


6.  Number  of  channels  being  digitized. 

7.  Playback  ratio. 

8.  Duty  Cycle. 

9.  Duty  cycle  reference  time. 

Label  Data 

In  addition  to  the  above  initialization  data,  the  following  parameters 
will  be  written  on  the  output  magnetic  tape  and  used  to  identify  the  digitized 


data: 

1.  Cruise  No. 

2.  Deployment  number 

3.  Mooring  type 

4.  Ship  name 

5.  Launch  day  number 

6.  Start  time  of  data  recording  (Greenwich  time) 

7.  Stop  time  of  data  recording  (Greenwich  time) 

8.  Recovery  day  number 

9.  Duty  cycle  (percent) 

10.  Date-time  group  of  start  of  on  time,  if  duty  cycle  is  not  equal  to  100% 

11.  Time  compression  (example:  10:1  playback  rate) 

12.  Playback  tape  speed 

13.  Mooring  latitude  and  longitude 

14.  Hydrophone  serial  numbers  (6  total). 

15.  Analog  tape  identification  number 

16.  Comments:  comments  inserted  at  run  time. 


Card  inputs  are  used  to  determine  what  portions  of  the  data  will  be 


digitized  in  order  to  (1)  eliminate  all  shot  intervals  so  that  only  ambient 
noise  data  is  digitized,  or  (2)  eliminate  all  umhient  noise  portions  in  order 
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to  digitize  only  shot  intervals.  The  Brush  recorder  analysis  may  be  used 
to  select  those  intervals  from  the  ambient  noise  tape  which  will  be  digitized. 

A  series  of  cards  are  keypunched  prior  to  the  digitizing  process  de¬ 
fining  which  sections  of  the  analog  record  are  to  be  digitized.  Hach  card 
contains  the  following  information: 

Tq  -  (Start  time):  in  days,  hours,  and  minutes  of  the  start  of  the 
record  to  be  digitized. 

Ts  -  (Stop  time):  in  days,  hours,  and  minutes  defining  the  time  on  the 
analog  record  to  stop  digitizing. 

T  -  (Burst  length) :  in  seconds  defining  the  total  time  of  the  sample 
record,  where  T  is  the  total  record  length  =  number  of  samples  X 
sampling  interval.  T  is  a  measure  of  the  frequency  resolution  de¬ 
sired  (seconds  in  hydrophone  time). 

Ti  -  (Burst  interval):  in  seconds  defining  the  interval  between  start 
time  of  each  burst  (in  hydrophone  time) . 

Figure  2S  illustrates  an  example  where  a  section  of  the  analog  record 
is  to  be  digitized  from  times  1230  to  1410  on  day  288.  Since  a  one  hertz 
resolution  is  desired,  T  =  1  second.  The  burst  interval  T^  is  picked  at 
10  seconds.  The  number  of  data  sets,  S,  generated  from  digitizing  in  this 
mode  is  given  by: 

c  -  VT0  -  Cl  10)  (60) 

T.  10 

l 

S  =  660 

Therefore,  660  data  sets  will  be  generated  and  stored  on  output  tape.  This 
would  provide  660  data  blocks  for  time  series  analysis  where  each  data  block 
contained  1166  data  points  per  hydrophone. 
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BURST  MODE  (Ambient  Noise  Anclysis)  June  1975 

Determining  the  Ntmber  of  Data  Sets  Generated  in  Burst  Mode 

Day  288  Day  288 


Digitizing  the  analog  record  in  burst  mode 
Total  data  sets,  S,  accumulated. 

r  _  Ts"T0  110(60)  _ 

S  *  -y- —  =  — — —  =  660  data  sets 

where  Tg  =  stop  time  -  1410 
Tq  =  start  time  =  1230 

*  burst  interval  =  10  seconds  (hydrophone  time) 
T  *  burst  length  -  1  second  (hydrophone  time) 


Figure  25.  8urst  Node  (Ambient  Noise  Analysis) 


The  output  tape  generated  as  a  result  of  digitizing  contains  a  tape 
label  record  plus  sets  of  data  records,  bach  burst  is  represented  on  the 
output  magnetic  tape  by  a  set  of  records.  Hie  first  part  of  each  data  record 
consists  of  a  label  containing: 

Start  t^me  of  burst 
Record  number  within  burst 
End-of-burst  indicator 

For  a  given  burst,  the  record  number  within  the  burst  will  start  at  1 ,  and 
increase  moootonically  thus  enabling  the  program  to  determine  the  number  of 
records  within  the  burst.  The  last  record  in  the  burst  will  have  the  end-of- 
burst  indicator  set  to  one  (1);  otherwise  this  indicator  will  be  zero.  The 
next  burst  continues  in  the  same  manner,  with  no  end-of-files  between  bursts, 
except  that  the  record  count  starts  again  at  1  and  increases  monotonically. 

An  end-of-file  mark  will  be  written  at  the  end  of  all  bursts  on  the  tape. 

10.8  Duty  Cycle  Considerations 

Since  either  a  continuous  or  a  discontinuous  duty  cycle  is  expected 
for  ACODAC  moorings,  the  program  has  the  ability  to  determine  whether  the 
intervals  to  be  digitized  in  the  burst  mode  represent  a  continuous  time 
series  during  the  on-tirae  of  the  ACODAC  recording  process.  Referring  to 
Figure  26,  it  can  be  seen  that  if  the  time,  T,  the  burst  length,  includes  the 
off-time  of  the  ACODAC  tape  recorder,  a  discontinuous  time  series  would  re¬ 
sult.  To  circumvent  the  problem,  the  program  computes  the  valid  on-times 
of  a  series  by  using  two  parameters  provided  during  the  initialization  phase. 

1)  Duty  Cycle 

2)  Date-time  group  of  start  of  on-time  of  recorder 
For  example,  if  the  duty  cycle  for  an  ACODAC  mooring  is  50%  (7  minutes  on, 

7  minutes  off)  and  if  the  start  time  of  the  first  on-time  minute  is  given. 
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DIGITiZiNG  WITH  DISCONTINUOUS  DUTY  CYCLE 


k  Recorder  Off 


Discontinuous  tine  scries  resulting  from  system 
digitizing  analog  record  through  ACODAC  recorder 
turn  off  period. 

Duty  Cycle  =  50*.  (7  minutes  on  -  7  minutes  off) 

T  =  Burst  Length 


Figure  26.  Digitizing  with  Discontinuous  Duty  Cycle 
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the  program  can  compute  any  valid  on-time  period  in  the  analog  record. 

A  test  is  made  using  the  start  and  stop  times  of  the  burst  in  question 
to  deteimine  if  any  ACODAC  recorder  off  time  is  present  during  this  period. 

If  any  off  time  is  found,  the  entire  burst  shall  be  rejected. 

10.9  Digitizing  Data  from  One  to  Six  Hydrophones 

The  program  has  the  ability  to  digitize  data  from  one  to  six  hydro¬ 
phones  connected  to  the  analog  inputs.  The  zero  axis  crossing  detector  has 
a  hydrophone  selection  switch  as  part  of  the  unit  which  is  used  to  select 
the  nunber  of  hydrophones  connected  for  digitizing.  Selecting  the  total 
number  of  hydrophone  channels  will  effectively  change  the  digitizer  sampling 
rate  by  the  number  of  hydrophones  connected. 

1 1 •  TIME  series  analysis  or  NARROW  band  spectra 

Once  the  digital  tape  is  generated  from  the  process  described  in  paragraph  10 
above,  the  output  tape  can  then  be  entered  in  the  noimal  SIGMA  7  batch  processing 
environment  for  further  processing.  The  Sigma  7  library  contains  standard  time 
series  analysis  programs  that  can  be  used  to  obtain  narrow  band  frequency  spec¬ 
tra  using  fast  fourier  transform  techniques  to  produce  standard  plots  of  frequency 
distribution,  cross  andauto  correlation  studies.  Figure  27  shows  a  software  flow 
diagram  of  the  system. 

12.  DIGITAL  PROCESSING  OF  TRANSMISSION  LOSS  (SHOT  ANALYSIS) 

12. 1  Background 

Digital  processing  of  shot  arrivals  was  first  developed  for  project 
NEAT  for  the  R/V  KNORR  cruise  24.  Transmission  loss  data  was  successfully 
acquired  and  processed  using  a  Hewlett  Packard  2116  digital  computer.  Since 
then  the  program  has  been  revised  to  process  ACODAC  shot  data  using  basically 
the  same  methods  employed  as  before  with  the  exception  that  shot  arrivals 


Cross 
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for  ACODAC  are  recorded  on  the  ambient  noise  analog  tape  whereas  the  shots 
from  Project  NEAT  were  processed  in  real  time. 

This  report  sunmarizes  the  two  processing  phases  required  to  obtain 
transmission  loss  data  from  shot  arrivals.  Phase  I  describes  the  digitiz¬ 
ing  of  shot  energy,  and  phase  II  descripes  the  processing  of  the  digital 
data  to  obtain  transmission  loss.  Since  the  program  is  not  fully  operational 
as  of  this  writing,  a  more  detailed  program  report  will  follow  at  a  later 
date. 

12.2  Phase  I  -  Analog-to-digitai  Conversion 

Phase  I  involves  the  digitizing  of  signals  from  one  hydro¬ 
phone  through  two  third  octave  filters  to  obtain  energy  levels  of  the  50 
and  100  Hz  components.  The  digitized  voltage  levels  from  the  filtered 
channels  are  squared  and  partially  integrated  in  the  computer  and  the 
results  written  on  magnetic  tape.  Time  code  and  gain  state  information  is 
also  brought  in  to  the  computer  and  used  in  the  output  magnetic  tape  records. 

The  digitizing  and  recording  process  produces  a  series  of  magnetic 
tape  records  where  each  record  consists  of  1024  words  plus  a  trailer  of 
15  words.  The  summing  and  squaring  process  produces  a  double  word  integer 
representing  the  sum  of  squares  of  n  data  points,  where  n  is  an  input  para¬ 
meter  and  is  typically  set  at  80.  Each  record  thus  represents  128  data 
points.  A  12  KHz  sampling  rate  is  used  for  digitizing  two  input  channels 
at  a  300  Hz  nyquist  sampling  rate  with  20:1  playback  time  compression. 

Each  record  contains  an  identifying  trailer  consisting  of  the  following 
parameters: 

1.  Sta’-t  time  at  beginning  of  record. 

2.  Gain  setting  of  hydrophone  data  amplifier  at  beginning  of 

record. 

3.  Gain  setting  of  hydrophone  amplifier  at  end  of  record. 
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4.  Deep/shallow  shot  indicator. 

5.  The  time  of  the  shot. 

6.  Two  indicators  showing  whether  an  overload  condition  or  a 
calibration  cycle  was  encountered. 

7.  The  bit  configuration  of  a  switch  register  which  is  used 
as  a  sense  switch  input. 

During  the  tape  playback  process,  an  oscilloscope  and  loud  speaker 
system  is  used  to  locate  each  shot  arrival  by  listening  to  the  signal  be¬ 
ing  digitized.  At  a  20:1  playback,  a  shot  arrival  through  the  system 
sounds  like  a  clicking  noise  generally  distinguishable  from  the  background 
noise.  When  a  shot  is  recognized,  a  push-button  switch  is  used  to  gen¬ 
erate  an  interrupt  to  the  computer  which  flags  the  record  being  accumulated 
as  containing  a  shot. 

The  time  at  which  the  button  was  pushed  is  then  stored  on  digital 
tape  along  with  the  other  time  and  gain  information. 

12.3  Phase  II  -  Data  Analysis 

In  Phase  II  the  analysis  program  reads  the  digital  data  records  from 
Phase  I  and  searches  for  records  which  have  been  flagged  as  containing  a 
shot  arrival.  The  flagged  record  is  then  displayed  for  the  operator  on  an 
oscilloscope  in  a  format  which  is  essentially  acoustic  energy  as  a  function 
of  time  for  the  period  encompassed  by  the  record.  A  typical  display  is 
shown  in  Figure  28.  Ambient  noise  produces  a  gentle  positive  slope  to  the 
line  of  points.  The  arrival  of  a  shot  produces  a  steeper  positive  slope, 
the  steepness  depending  on  the  signal  to  noise  ratio.  Manually  controlled 
cursors  are  also  displayed  and  are  used  to  identify  an  interval  of  the 
record  which  may  be  used  to  determine  an  ambient  noise  power  level  which 
is  subsequent ly  displayed.  Cursors  on  the  ’’noise  free"  display  are  then 
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aligned  to  encompass  the  shot,  and  the  energy  content  of  the  shot  is  deter¬ 
mined,  and  the  program  computes  the  transmission  loss  for  the  shot.  One 
of  the  cursors  is  also  used  to  get  a  more  accurate  shot  arrival  time  by 
aligning  it  on  the  beginning  of  the  shot  in  the  display. 
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APPEND1X  A 

INSTRUCTIONS  FOR  USING  PROGRAM  ANPLOT  TO  PLOT 
ACODAC  AMBIENT  NOISE 

These  instructions  show  how  to  make  a  standard  time  series  plot  using 
the  XD  series  ACODAC  ambient  noise  edited  tapes. 

The  following  deck  of  instruction,  control,  and  data  cards  is  used: 

1.  !JOB  card 

Computer  control  card  specifying  user  nunber,  priority  and  account 

number. 

2.  ! MESSAGE  cards 

iMESSAGE  PLOT  JOB  ID  _ 

!  MESSAGE  INPUT  FROM  xn  _ -9T 

3.  ! ASSIGN  CARDS  and  IFTC'jPY 

! ASSIGN  F: PLOT,  (DEVICE,  7T) ,  (SN,  PLT1) ,  (OUT),  (BIN), 
(UNPACK),  (TRIES,  10) 

! ASSIGN  M: El,  (DEVICE,  9T),  (SN,  XL  _ ) 

! ASSIGN  M:E0,  (FILE,  TEMP) 

! FTCOPY 

! ASSIGN  F:30,  (FILE,  TEMP) 

4.  !OLAY  card 

(Assuming  that  anplot  and  its  subroutines  are  stored  on  the 
disk  in  the  ACODAC  account.) 

!OLAY  (EF,  (ANPLOT)),  (UNSAT,  (3)),  (MAP) 

5 .  ! RUN  card 
!  RUN 


6.  !DATA  card 
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7.  Run  time  comment  card 

72  characters  on  this  card  will  be  reproduced  at  the  top  of  the  plot. 
For  example:  PLOTTED  6  FBB  73  FROM  XD11-9T  L.  BAXTER,  II 
7a.  Run  time  header  card 

72  characters  reproduced  below  run  time  comment 

8.  Unit  reference  number  card 

In  format  (13)  punch  unit  reference  number  for  input  file. 

Example:  30. 

9.  Printout  parameters  card 

In  format  (5X,  SIS)  punch  LXFMX,  LXPMX,  LXNMX,  LXQMX,  LFLMX. 

Where  parameters  control  print-out  as  follows: 

LXFMX  «  Maximum  count  for  LXF  where  output  on the  printer  will  occur  lor  data 
off  the  high  end  of  the  time  scale  with  LXF  *  1  to  LXFMX. 

LXPMX  =  Maximum  count  for  LXP  reading  higher  or  lower  than  the  extremes  of 
the  decibel  scale  will  not  be  plotted.  Output  to  the  printer  will 
occur  for  LXP  =  1  to  LXPMX. 

LXNMX  *  Maximum  count  for  unconditional  printing  of  values  from  input  tape. 

Values  from  input  tape  will  be  printed  for  LXN  =  1  to  LXNMX. 

LXQMX  =  Maximum  count  for  attempts  to  read  header  from  input  tape.  Attempts 
to  read  header  will  be  made  until  it  is  read  correctly  or  until 
LXQ  =  LXQMX. 

LFLMX  =  Maximum  count  for  output  checking  DATRD,  This  output  will  occur 
for  LFL  *  1  to  LFLMX. 

The  card  used  at  present  has  the  following: 

20  100  20  5  10 

10.  Factor  card 

Scales  of  the  basic  plot  are  multiplied  by  this  factor  in  format  (F4.2) 

2.50  fills  the  plotter  paper 

1.50  gives  a  smaller  plot  and  reduces  PDP5  plotting  time. 


_  i»VJ' >^»  <*  '■"•*^**  **  •t*f>k:' 
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11.  Selector  card 

Usual  lower  db  limits  T  or  F 

If  this  card  is  T,  leave  out  the  next  card  and  minimum  db  plotted 
in  db  re  1  ybar  will  be  as  follows: 

(db  re  1  y  Pascal) 


XMIN(l)  =  -20.0 

80 

XMIN(2)  =  -20.0 

80 

XMIN(3)  *  -20.0 

80 

XMIN(4)  =  -25. 0 

75 

XMIN(5)  =  -25.0 

75 

XMIN(6)  =  -25.0 

50  Hz 

75 

XMIN(7)  =  -30.0 

70 

XMIN(8)  =  -30.0 

70 

XMIN (9)  =  -30.0 

70 

XMIN(10)  =  -31.0 

125  Hz 

65 

XMIN (11)  =  -35.0 

65 

XMIN (12)  =  -40.0 

200  Hz 

60 

XMIN (13)  =  -40.0 

60 

XMIN (14)  =  -45.0 

55 

XMIN (IS)  =  -50.0 

50 

XMIN(16)  =  -50.0 

50 

12.  Lower  DB  limit  card  in  format  (16F  4.1) 

If  one  wishes  to  change  any  of  the  above  values,  card  section  11 
above  must  be  F  and  card  13  must  have  a  value  entered  for  each  channel  that 
will  be  plotted.  Channels  not  plotted  will  be  changed  to  zero  by  blanks 
but  this  will  not  have  any  effect  on  the  plot.  Note:  in  present  version 
of  ANPLOT  all  values  except  the  final  plot  are  in  db  re  1  ybar  -  the  final 
plot  is  in  db  re  1  y  Pascal. 

13.  Selector  card 

Usual  scales  -  T  or  F 

Note:  the  scales  are  determined  by  several  parameters: 

Factor  has  already  been  set  above  and  determines  over-all  size  of  plot. 
The  dimensions  given  below  are  those  for  factor  *  1.0. 
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Variable 

Value 

name 

for  t 

Remarks 

SCALY 

36.0 

hours  per  iach 

SCALX 

30.0/2.0 

decibels  per  inch 

YAXL 

8.5 

length  of  time  axis  in  inches 

XAXL 

2.0 

length  of  each  decibel  axis  in  inches 

XDISP 

20.0/30.0 

distance  between  time  series  plots 
in  inches 

XOFS 

1.5 

distance  from  origin  to  first  time 
axis  in  inches 

ISCAL 

24 

hours  per  division 

DBTNT 

1.0 

decibel  interval  for  histogram 

HIMAG 

1.0 

HIOFS 

1.5 

14.  Scales  card  (first  in  format  (6  F  10.2) 

SCALY,  SCALX,  YAXL,  XAXL,  XDISP,  XOFS 
for  24  hour  plot  this  card  will  read 

300  1500  850  200  67  150 

15.  Scales  card  (second) in  format  (15,  3F10.2) 

ISCAL,  DBINT,  HIMAG,  HI0FS 

for  13  day  plot  this  card  will  read 
3  100  100  100 

16.  Selector  card 

Channel  order  to  plot  USUAL  CHANNEL  ORDER  T  or  F  USE  "F"  Note: 
this  gives  control  of  channels  to  plot  to  the  following  cards  otherwise 
(in  current  version)  prog  would  attempt  to  plot  1  thru  16  and  there  are 
only  15  channels  of  data  returned  by  DATRD  on  ACODAC  TAPES  the  16th  is  set 
to  zero. 

17.  Channel  plot  control  cards  -  In  format  (318) 
punch  ICHPL  (I),  IFRK  (I),  ID?  (I) 

where  ICHPL  (I)  is  integer  1  to  15  (or  999) 

IFRK  (I)  is  frequency  in  Hz 
I DP  (I)  is  depth  in  meters 


R 


H  i 


I 
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Note:  where  all  channels  are  the  same  depth  as  they  are  when  plotting  one 
hydrophone  at  different  frequencies,  leave  IDP  blank  and  it  will  not  appear 
on  the  separate  time  series  plots. 

Note:  if  ICHPL  (I)  *  999  that  card  terminates  the  series  of  channel  plot 
control  cards 

Example  of  channel  plot  series: 


6 

10 

12 

999 


50} 

125} 

200} 


4  cards 


3  channels  will  be  plotted  -  one  each  for  50,  125,  and  200  Hz 

18.  Origin  card  -  In  format  (215) 

IDABG,  IHRBG 

where  IDABG  is  an  integer  0  to  365 
and  IHRBG  is  an  integer  0  to  23 

(Giving  the  day  of  the  year  and  hour  of  the  day  for  the  origin  of  the  time 
axis) 

19.  Selector  card  T  or  F 
T  =  NO  DELETIONS 

F  *  YES  THERE  ARE  DELETIONS 
at  run  time  -  if  T  skip  (i.e.  leave  out)  21.  and  22. 

20.  Card  for  mrober  of  deletions  -  In  format  (12) 

punch  an  integer  1  to  20  for  number  of  run  time  deletion  intervals 
desired 

21.  Parameter  card  for  each  deletion  interval 
In  format  (24,  IX,  13,  IX,  14,  IX,  14) 
punch  ICHWD,  IDAY,  IIIRM1,  IHRM2 

ICHWD  is  a  hexidecimal  number  representing  the  16  bit  binary  with  1  bit  for 
a  deleted  channel  0  bit  for  non-deleted  channel  where  channel  in  this  case 
refers  to  the  time  series  in  order  plotted.  Thus  to  delete  an  interval  in 
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the  first  and  third  channels  plotted  would  call  for  ICHWD  =  5.  To  delete 
in  the  first,  fourth,  fifth,  and  nineth  ICHWD  =  119.  To  delete  in  the  first, 
third,  fourth,  and  fifth  channels  plotted.  ICHWD  =  ID 
I DAY  is  an  integer  0  to  365  representing  the  day  of  the  year  on  which  the 
deletion  begins. 

IHRM1  and  IHRM2  are  4  digit  nunbers  representing  hour  and  minute  the  dele¬ 
tion  begins  and  hour  and  minute  it  ends.  If  IHRM2  is  less  than  IHRM1  the 
program  assumes  that  the  deletion  ends  on  the  following  day  otherwise  both 
hour  and  minute  nunbers  are  assumed  to  be  on  the  same  day.  There  will  be 
one  parameter  card  for  each  deletion  interval  required.  That  is  the  number 
of  parameter  cards  equals  the  number  of  deletions  given  in  card  21. 

22.  Selector  card 

F  to  terminate  run  with  this  plot  and  file 
T  to  plot  another  graph  on  another  file  of  plot  tape, 
normally  use  F 

T  would  be  used  only  if  one  desired  in  the  same  run  with  the  same 
input  tape  to  re-enter  the  program  tu  plot  for  example  a  different  set  of 
frequencies  on  another  plot  entirely.  In  this  case  additional  data  cards 
must  be  provided  like  those  from  8  through  23,  etc.  finally  terminating  the 
last  file  and  the  run  when  a  23  selector  card  is  punched  with  an  F. 
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APPENDIX  B 

PROGRAM  DESCRIPTION:  ACOSTAT 


NAME:  ACOSTAT 

TYPE:  Main  Program  and  subroutines 

PURPOSE:  To  provide  statistical  and  graphical  output  of  the  ACODAC  data. 

MACHINE:  Xerox  Sigma  7 

SOURCE  LANGUAGE:  FORTRAN  IV  (Extended) 

PROGRAM  CATEGORY:  Data  processing  -  analysis 

DESCRIPTION:  The  program  system  ACOSTAT  is  used  to  process  the  Sigma  7  digital 
tape  of  the  1/3  octave  ACODAC  data.  The  user  has  the  option 
of  outputting  the  statistical  summary  of  selected  frequency 
bands,  db  distribution  of  selected  frequency  bands,  and/or 
the  plots  of  all  frequency  bands.  All  summaries  are  broken 
down  into  12  hour  intervals.  All  edited  points  (|db|^95.0) 
are  excluded  from  the  output  summaries. 

INPUT:  The  input  is  from  the  Sigma  7  ACODAC  digital  data  tape  and  from  data 
cards . 

Card  Input 

Card  1  -  Input  tape  record  structure 


LHD 

Colunns 

1-10 

length  of  header  identification 
record  (number  of  Sigma  7  words) . 

LID 

Columns 

11-20 

length,  in  Signa  7  words,  of  hour/ 
minute  data  calibration  identifier. 

IPM 

Columns 

21-30 

number  of  integration  periods  per 
minute. 

MPR 

Columns 

3:  -40 

the  duty  cycle,  or  nunbar  of  minutes 
of  data  per  record. 

NBN 

Columns 

41-se 

number  of  frequency  bends  of  data. 

All  of  the  above  are  integer  values,  right  justified  in  the  field. 

Card  2  -  Unit  device  code  for  digital  tape  mount. 

NT  Columns  1-10  unit  device  code  of  unit  on  which 

Sigma  7  digital  data  tape  is  to  be 
mounted.  This  must  be  the  same  as 
the  DCB  assignment  on  the  J ASSIGN 
card  for  the  tape.  Integer  value, 
right  justified. 
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Card  3  -  Plot  option  specification 

IPLT  Columns  1-3  YES  for  plotting;  blank  or  NO  to 

delete  plotting 

DBIN  Columns  11-20  the  plot  scale,  in  db's  per  inch. 

Card  4  -  Statistics  option  specification 

NST  Columns  1-10  right  justified  integer  specifying 

the  number  of  bands  for  which 
statistical  summaries  are  desired. 

If  blank  or  0  (zero)  no  statistics 
are  output. 

Card  4  ♦  . . . . (NST/8) 

IBAN1-NST  If  NST  >  0,  the  next  cards  contain 

the  frequency  band  numbers  for  which 
statistical  summaries  are  desired. 

Each  band  number  is  entered  as  an 
integer,  right  justified  in  fields 
ten  columns  wide.  The  value  NST 
specifies  the  number  of  band  numbers 
to  input.  If  more  than  8  were  speci¬ 
fied,  a  second  card  is  necessary. 

Card  5  ♦  ...  +  (NST/8)  -  Histogram  option  specification 

IHIST  Columns  1-10  right  justified  integer  specifying 

the  number  of  bands  for  which 
histogram  summaries  are  desired.  If 
blank  or  0  (zero)  no  histograms  are 
output . 

Remaining  Cards 

IHBAN1_IHIST  If  IHIST  >  0,  the  next  cards  contain 

the  frequency  band  numbers  for  which 
histogram  summaries  are  desired. 

Each  band  number  is  entered  as  an 
integer,  right  justified  in  fields  ten 
columns  wide.  The  value  IHIST  speci¬ 
fies  the  number  of  band  numbers  to  be 
input.  If  more  than  S  were  specified, 
a  second  card  is  necessary. 

Tape  Input 

The  input  tap*  is  the  Sigma  7  ACODAC  digital  data  tape  which 
was  created  from  program  AC0EDIT.  The  record  structure  of  the  binary  tape 
is  described  below. 
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Record  1  is  the  header  identification  record.  Its  word  length  is 
specified  by  the  input  variable  LHD. 

Records  2  to  the  end  of  file. 

Word  1  packed  word  containing  the  integer  day  number  in  the  leftmost 
16  bits  and  the  integer  hour/minute  identifier  in  the  right¬ 
most  16  bits. 

Word  2  packed  word  containing  the  integer  value  of  the  data  type 

identifier  in  the  leftmost  16  bits  and  the  integer  value  of 
the  gain  setting  integer  value  in  the  rightmost  16  bits. 

Word  3  real  value  of  the  transfer  function. 

Words  4-18  edited  db  data  valuer  for  each  of  the  frequency  bands  during 
the  first  integration  period  of  the  hour/minute  specified  in 
word  1. 

Words  19-33  same  as  above  for  second  integration  period. 

Words  34-48  same  as  above  for  third  integration  period. 

Words  49-63  same  as  above  for  fourth  integration  period. 

Words  64-78  same  as  above  for  fifth  integration  period. 

The  format  of  words  1-78  is  then  repeated  for  the  next  hour/minute 
identifier  in  the  record.  There  are  currently  seven  such  cycles  per 
record. 

Note  that  the  record  structure  above  is  determined  by  the  card  input  data. 

OUTPUT :  The  output  is  determined  by  the  user.  Run  information,  statistics, 
and  histograms  are  done  by  the  line  printer;  plotting  is  done  by 
a  Cal comp  plotter. 

Automatically  output  to  the  line  printer  is  a  run  identification. 

An  example  is  shown  in  Figure  20.  If  the  statistics  option  was 
chosen,  every  twelve  hours  the  statistics  for  those  frequency 
bands  selected  are  output  as  shown  in  Figure  21.  For  the  histogram 
option,  every  twelve  hours  the  histogram  of  db  distribution  is 
output  for  each  of  the  selected  frequency  bands.  An  example  of 
the  histogram  is  shown  in  Figure  22, 

A  twelve  hour  summary  of  all  frequency  bands  may  also  be  output 
to  the  plotter.  Each  curve  is  determined  by  points  representing 
the  means  of  each  record.  Thus,  a  curve  will  have  one  point  for 
the  mean  power  of  all  valid  points  in  a  seven  minute  cycle. 
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where  P^  is  the  average  power  in  db 

n  is  the  number  of  accepted  points  in  the  duty  cycle 
(seven  minutes) 

USAGE :  The  program  is  treated  as  an  ordinary  FORTRAN  IV  main  program.  The 
job  setup  is  described  in  Figure  4.  If  the  plotting  option  is 
specified,  the  user  must  submit  the  appropriate  special  action  card 
with  the  job  deck.  The  first  message  card  requests  that  the  plot 
tape,  if  used,  be  saved.  The  user  may  wish  to  change  this  feature 
and  plot  directly.  The  Computer  Users'  Handbook  or  a  member  of 
the  IPC  staff  may  be  consulted  if  the  procedure  for  changing  is  not 
understood.  If  plotting  is  not  used,  the  F:PLOT  device  assignment 
should  be  changed  to  (DEVICE,  NO). 

Figure  4 

2J0B  aaaa,  uuuu,  p 

IMESSAGE  LABEL  AND  SAVE  7T,  2400  FT  TAPE  — —  FOR  PLOTTING 
IMESSAGE  MOUNT  — — ,  ENABLE  WRITE 

'ASSIGN  F:PLOT,  (DEVICE,  7T),  (SN,  — ),  (OUT),  (BIN),  (UNPACK) 
IMESSAGE  MOUNT  TAPE  — — ,  800  BPI  (DIGITAL  DATA  TAPE) 

1 ASSIGN  F: 1,  (DEVICE,  9T),  (SN,  — ) ,  (IN),  (BIN) 

I  FORTRAN  GO 

FORTRAN  program  deck  goes  in  here  with  subroutines 
IOLAY  (MAP),  (UNSAT, (3)),  (GO) 

I  RUN 
I  DATA 

card  data  goes  in  here 

IE0D 

RESTRICTIONS:  Provision  has  been  made  to  handle  only  one  input  data  tape 

per  run.  Currently  dimensioned  to  handle  20  frequency  bands, 

200  frequency  db  intervals,  40  points  per  band  per  cycle, 

10  minutes  per  cycle  (record) ,  200  plot  points  per  band  per 
twelve  hour  period. 

STORAGE  REQUIREMENTS:  18,4S010  words. 

SUBPROGRAMS  REQUIRED: 

PL0TID  -  used  to  set  up  all  the  titles  and  identifiers  used  on  plots 


STORP  (IND,  NBN,  IPR)  -  used  to  calculate  the  mean  power  of  a  cycle 


and  store.  Also  determines  frequency  dis¬ 
tribution. 

IND 

plot  point  index 

NBN 

nunber  of  frequency  bands 

IPR 

number  of  integrations  per  record. 

BPI.0T  (BEGHR, 

FINHR,  IND)  -  used  to  plot  the  twelve  hour  sunmaries. 

BEGI1R 

hour  at  beginning  of  plot 

FINHR 

hour  at  end  of  plot 

IND 

number  of  points  in  plot 
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CHKIN  -  used  to  ensure  that  plot  point  falls  within  boundary. 

BSTAT  (IC,  NCLS,  CLSINT,  XMIN,  XMAX,  CUL,  NX)  -  used  to  calculate 

the  statistics  for  the  twelve  hour  summaries. 
IC  array  index  of  current  frequency  band 

NCLS  number  of  intervals  in  distribution  range 

CLSINT  size  of  interval 

XMIN  minimum  frequency  db  value 

XMAX  maximum  frequency  db  value 

CUL  upper  limit  of  each  interval  in  range 

NX  population  of  each  interval 

CHISQ  (NT,  RM1,  SD6,  HI,  X2)  -  chi  square  test  of  fit  for  normal 

distribution. 

NT  total  number  of  points  for  band  in  summary 

RM1  mean  of  points 

SD6  standard  deviation  of  points 

NI  degrees  of  freedom  minus  3 

X2  chi  squared 

CDF  (NI,  RM1,  SD)  -  calculates  the  degrees  of  freedom  of  the  db 

distribution  curve. 

NI  number  of  degrees  of  freedom 

RM1  mean  db  of  point  in  frequency  distribution 

for  band 

SD  db  standard  deviation  of  points  in  frequency 

band 

ARNRM  (Z,  TOL)  -  calculates  the  area  under  a  normal  curve  for  the 

current  frequency  band  distribution  of 
db  values. 

Z  standard  deviation  of  enlarged  intervals 

TOL  tolerance 

STOUT  (MINB,  BEGHR,  FINHR)  -  outputs  the  statistics  calculated  in 

BSTAT. 

MINB  minimum  frequency  band  number 

BEGHR  hour  at  beginning  of  statistical  summary 

period 

FINHR  hour  at  end  of  statistical  svmmary  period 

HISTO  (MINB)  -  calculates  and  outputs  the  histograms 
MINB  minimum  frequency  band  number 

Additional  library  routines:  PLOTS,  PLOT,  SYMBOL,  NUMBER,  AXDRAW, 
ISA. 


OPERATIONAL  ENVIRONMENT: 


Device 


Function 


Special  requirements 


card  reader 
line  printer 
Sigma  7  ACODAC  digital 
tape 
plots 


input 

output 


F:10S 
F:  108 


input 

output 


user  specified 
F:PL0T 
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OPERATIONAL  CHARACTERISTICS: 

To  start,  the  user  must  include  the  appropriate  device  assignment 
cards  for  plotting.  The  card  data  is  then  read  in,  followed  by 
the  input  of  the  digital  tape  header  record.  After  initialization 
of  variables  and  determination  of  output  options,  a  run  identifica¬ 
tion  record  is  output  to  the  line  printer. 

A  data  record  is  then  input  from  the  digital  tape.  Its  time  is 
extracted  and  compared  to  the  current  summary  interval  limits.  If 
it  falls  within  the  limits,  its  location  in  the  frequency  disl 
tribution  table  for  the  band  is  determined,  along  with  the  mean  of 
the  valid  points  within  the  record.  This  mean  (and  its  standard 
deviation)  are  stored  for  plotting.  When  the  summary  period  (12 
hours)  is  complete,  the  options  requested  are  executed.  The  ap¬ 
propriate  variables  are  cleared  and  the  loading  procedure  is 
initiated  for  the  next  sequence. 

Only  those  points  with  absolute  db  values  less  than  95.0  are  con¬ 
sidered  as  contributing  toward  any  summaries. 

TIMING:  The  program  run  time  will  vary  depending  on  the  number  of  output 
options  requested  by  the  user.  The  following  estimates  may  be 
used  for  computing  run  times: 


program  compilation 
input  and  interval 
assignment 

plotting  (Sigma  7  time) 
statistics 

histograms _ 


plotting 


4  minutes 

1.25  minute/day’s  data 
0.25  minute/day's  data 
1.00  minute/day’s  data 
negligible _ 

approx.  25  minutes /day  (2  plots) 


ERRORS  §  DIAGNOSTICS: 


HEADER  RECORD  ERROR  -  STOP  152 

There  is  something  amiss  with  the  input  tape.  Either  the 
header  record  is  in  error  or  does  not  exist. 

READ  ERROR  -  LAST  DAY  III  JJJ  STOP  210 

An  input  error  was  encountered  in  reading  the  input 
digital  tape 

PLOTTING  COMPLETED 

The  plotting  for  the  current  stannary  period  was  success¬ 
fully  completed. 

STATISTICS  COMPLETED 


The  statistics  for  the  current  stanmaiy  period  were  success¬ 
fully  completed. 
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APPENDIX  C 

ACODAC  SIGMA  7  THIRD  OCTAVE  DATA  EDIT  PROGRAM 

REPLACES  TAPEWRITE 
Control  Cards 


Control  card  type  #  right  justified  in  columns  1-5  (15) 
Followed  by  15  5-bit  column  data  fields  (15F5.0) 

TYPE  0  Tine  window  and  processing  control  (must  be  last  card  of  set) 


Field  1-3 

start 

time:  days,  hours,  minutes. 

4-6 

stop 

time. 

7 

=  0: 
i  0: 

use  original  transfer  function, 
use  new  transfer  function. 

8 

<  0: 

>  0: 

use  time  window  to  control  processing, 
count  of  number  of  records  to  process 
(ignore  time  window). 

9 

<  0: 

>  0: 

do  not  skip  any  records  before  processing 
count  of  number  of  records  to  skip  before 

processing. 

TYPF.  1  Band  corrections 

Field  1-15  hand  corrections  for  individual  bands 
TYPE  2  Universal  correction 

Field  1  correction  to  all  bands 


TYPE  3  Control  information 


Field  1 


<_  0:  no  output  tape  dump. 

>  0:  number  of  records  of  output  tape  to  dump  at 
end  of  run. 


2 

3 


4-6 


7 


>  0:  new  fraction  to  be  used  in  heuristic  check. 

<_  0:  no  change  to  fraction  (default  *  .5). 

>  0:  maximum  day  number  to  believe  from  time 

code  reader. 

*  0:  no  change  to  maximum  day  number  (default* 12) . 

time  to  be  added  to  tape  time  for  output:  days, 
hours,  minutes. 

*  0:  pass  through  editing. 

i  0:  do  not  pass  through  editing. 


^  ■  -'-it.-  f-T*wr"hri» 
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8 

9 

10 

11 


=  0:  do  not  change  maximum  error  message  counter. 
i  0:  new  maximum  error  message  counter. 

=  0:  do  not  dump  input  records, 
y  0:  dump  input  records. 

*  0:  write  output  tape, 
y  0:  do  not  write  output  tape. 

=  0:  use  tape  time  for  time  windows, 
y  0:  use  corrected  time  for  time  windows. 


12 


=  0:  normal  processing. 

y  0:  give  user  subroutine  "PROC"  control  for  each 
data  minute. 


TYPE  4  Fields  1-15  number  of  minutes  to  be  used  in  each  of  1  to  15 

averages. 


TYPE  5  Fields  1-15  =0:  normal  processing 

y  0:  gives  user  routine  "TIMER"  control  when  cor¬ 
responding  average  is  computed. 


TYPE  6  Field  1 


logical  unit  nunber  from  which  to  read  subsequent 
control  cards. 


WHU.  Ml' » U'w  im 
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APPENDIX  D 

STAND  ALONE  GENERAL  PURPOSE  ANALOG  TO  DIGITAL 
CONVERSION  PROGRAM  FOR  AMBIENT  NOISE  ANALYSIS 
USING  A  HIGH  SPEED  MULTIPLEXER  DIGITIZER 


NAME:  SAGPAD  (Stand  Alone  General  Purpose  Analog  to  Digital) 

TYPE:  Stand  Alone  Main  Program 

PURPOSE:  The  program  is  designed  to  provide  A/D  conversion  on  the  Sigma  7 

for  ACODAC  analog  data  tapes.  Since  most  parameters  are  dynamically 
set,  it  can  be  used  as  a  general  purpose  A/D  program 

MACHINE:  Sigma  7 

SOURCE  LANGUAGE:  Metasymbol 

PROGRAM  CATEGORY:  Utility 

DESCRIPTION: 

SAGPAD  is  a  stand  alone  program  designed  to  be  loaded  and  run  with 
the  Xerox  stand  alone  loader  with  I/O  handlers.  The  program  reads 
analog  data  via  the  Sigma  7  analog  to  digital  converter  and  writes 
it  to  a  9-track  magnetic  tape.  The  operation  of  the  SAGPAD  may  be 
divided  into  two  phases: 

I.  Parameter  Input  and  Intitialization 

During  this  phase  the  program  reads  parameter  cards  which  de¬ 
fine  the  environment  in  which  it  is  to  operate,  the  modes  of  opera¬ 
tion,  and  identification  data  for  the  output  tape  header  record. 
Using  this  data,  the  program  calculates  its  operational  parameters 
and  sets  up  an  internal  clock  which  will  operate  at  a  rate  which 
reflects  real  time  with  respect  to  the  actual  collection  of  the 
data.  This  is  referred  to  as  hydrophone  time.  For  example,  if 
the  playback  ratio  is  10:1,  the  clock  will  be  set  to  operate  at 
10  times  wall  clock  time. 

If  the  user  is  digitizing  by  a  predetermined  schedule,  based  on 
time,  the  program  will  next  read  a  series  of  shot  cards  and  set 
up  a  schedule  of  operations. 

II.  Digitizing 

This  phase  digitizes  the  data  and  writes  it  out  on  magnetic  tape. 

If  the  manual  mode  is  being  used,  the  program  digitizes  a  burst 
every  time  the  user  presses  the  interrupt  button  on  the  A/D  unit. 

If  in  automatic  mode,  the  program  compares  the  digitizing  schedule, 
set  up  in  phase  1,  against  its  internal  clock  and  digitizes  as 
the  schedule  directs. 
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If  the  analog  tape  contains  an  I RIG  "CH  time  code  channel,  the 
user  may  direct  the  program  to  read  it  between  shots  for  the  pur¬ 
pose  of  updating  the  internal  clock.  This  minimizes  differences 
between  the  internal  clock  and  hydrophone  time  which  are  bound 
to  accumulate  after  a  period  of  time  due  to  playback  rate  inac¬ 
curacies  . 

The  actual  digitizing  is  done  via  a  read  direct  loop  with  no 
interrupt  being  used.  This  assures  the  fastest  possible  rate, 
up  to  100,000  15-bit  data  points  a  second.  Since  this  rate  is 
far  too  fast  to  allow  direct  transfer  to  magnetic  tape,  the  data 
is  temporarily  stored  on  disk  during  the  digitizing.  Two  core 
buffers  are  used.  While  the  data  is  being  stox-ed  into  one  buffer, 
the  other  is  being  written  to  disk.  When  the  input  buffer  is 
full,  the  buffers  are  switched.  This  process  continues  until  the 
digitizing  burst  is  done. 

'When  the  burst  is  finished,  the  data  stored  on  the  disk  is  read 
back  into  core,  split  up  into  tape  record  size,  and  written  to 
magnetic  tape.  When  this  is  done  the  program  is  ready  to  digitize 
another  burst.  When  all  bursts  call  for  are  done,  the  program 
writes  a  double  end-of-file  on  the  output  tape,  backspaces  over 
the  second  one  and  terminates.  This  allows  more  than  one  data 
file  to  be  stored  on  one  tape,  if  desired. 

INPUT;  I.  Card  Reader 

A.  General  Parameter  Cards 


These  cards  define  general  operating  parameters.  Each  card 
is  identified  by  a  key  which  is  always  four  characters  long. 
The  card  foxmat  is  as  follows,  starting  in  colunn  1: 

Key,  Pj  [ , . . . ,  Pjq] 

Imbedded  blanks  are  not  allowed  except  in  character  strings. 
The  brackets  denote  optional  values  on  multi-parameter 
cards.  They  do  not  appear  on  the  card. 

There  are  three  types  of  parameters. 

1.  Character  string  -  there  may  be  any  EBCDIC  characters 
up  to  the  maximum  call  for  by  the  particular  item. 

They  are  stored  as  is. 

2.  Integer  -  these  are  represented  by  decimal  numeric 
characters  only  and  are  stored  as  32-bit  binary  integers. 
A  sign  may  precede  the  number.  Examples:  6 

-100 

5025 

3.  Real  -  these  are  represented  as  a  decimal  numeric  string 
which  may  contain  a  decimal  point.  The  numbers  may  be 
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signed.  They  are  stored  as  Sigma  7  single  precision 
floating  point  numbers.  1.5 

-2.6 

35. 

0.51 

.75 

25. 

The  following  table  lists  the  allowable  parameters. 

The  column  labelled  KEY  gives  the  4-character  identifier  which 
is  punched  in  columns  1-4  of  the  card. 

TYPE  refers  to  parameter  vype  where  C  stands  for  EBCDIC 
characters,  I  for  integer  and  R  for  real. 

The  column  labeled  NUMBER  refers  to  maximum  number  of  char¬ 
acters  if  it  is  a  C  parameter  or  the  maximum  number  of  parameter 
values  if  it  is  I  or  R. 

The  DEFAULT  column  shows  the  initial  default  value  which  will 
be  used  if  the  item  is  not  specified.  Input  values  for  a  run 
will  become  the  default  values  for  the  next  run  unless  the 
program  is  reloaded. 

The  DESCRIPTION  column  describes  the  parameter.  Descriptions 
preceded  by  an  asterisk  are  not  used  by  SAGPAD  but  are  output 
as  part  of  the  header  record  for  subsequent  use  or  identifica¬ 
tion.  The  description  of  these  items  refers  to  their  use  by 
subsequent  ACODAC  processing  programs.  Non-ACODAC  users  may 
use  them  for  any  purpose  they  wish. 

B.  Comment  Cards 


The  user  is  allowed  up  to  five  cards  (400  characters)  of 
comments  which  will  be  placed  at  the  end  of  the  output  tape 
header  record.  The  set  up  of  the  comment  deck  is  as  follows: 

:C0M  comment  Key  card  up  to  five  cards  of  comments 
:END  signals  end  of  comments. 

C.  Shot  Specification  Cards 

Shot  cards  are  used  to  set  up  a  schedule  of  digitizing  when 
operating  in  the  automatic  mode.  There  may  be  up  to  100  shot 
cards.  When  operating  in  the  manual  mode  oily  one  (1)  shot 
card  is  input  and  only  the  Burst  Length  is  specified. 

The  set  up  of  the  shot  deck  is  as  follows: 

:SU0T5  signals  shot  cards  follow  deck  of  shot 
specification  cards 
:END  signals  end  of  shot  cards. 
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GENERAL  OPERATING  PARAMETERS 


KEY 

TYPE 

NO. 

DEFAULT 

DESCRIPTION 

TPID 

C 

9 

Blanks 

•Output  tape  ID. 

RLEN 

I 

1 

2048 

Output  data  record  length  in  bytes, 
exclusive  of  header  words. 

NPPW 

I 

1 

2 

Number  of  data  points  from  A/D  con¬ 
verter  per  transfer  (1  or  2). 

GAIN 

R 

1 

2.5 

Gain  setting  for  digitizer.  Only 

1.25,  2.5,  5.0,  or  10.0  allowed. 

TSCN 

I 

1 

0 

Type  of  scan.  0  =  Sequential,  1  = 

Random. 

CTRL 

I 

1 

1 

Digitizing  control.  0  *  Manual, 

1  =  Automatic  using  preset  time  schedule 
for  digitizing  bursts. 

RATE 

I 

1 

900 

Digitizing  idte  per  channel  in  hydro¬ 
phone  time. 

NCHN 

I 

1 

6 

Number  of  channels  being  digitized  (1-6). 

PBCK 

R 

1 

10.0 

Playback  rate  of  analog  tape. 

DPNO 

C 

4 

Blanks 

•DEPLOYMENT  NO. 

MOOR 

C 

4 

Blanks 

•MOORING  TYPE. 

LDAY 

C 

4 

Blanks 

•LAUNCH  DAY. 

YEAR 

C 

4 

Bi°r.ks 

•LAUNCH  YEAR. 

SGMT 

I 

3 

0,  0,  0 

Recording  start  time  (GMT) ,  Day,  Hour, 
Minute. 

EGMT 

I 

3 

O 

o 

o 

Recording  end  time  (GMT) ,  Day,  Hour, 
Minute. 

ETON 

I 

1 

0 

•Duty  cycle  on  time. 

DTOF 

T 

A 

1 

0 

•Duty  cycle  off  time. 

RTSP 

R 

1 

0.0 

•Recording  tape  speed. 

LATI 

C 

12 

Blanks 

•Mooring  Latitude 

LONG 

C 

12 

Blanks 

•Mooring  Longitude 
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KEY 

TYPE 

NO. 

DEFAULT 

DESCRIPTION 

HSTA 

1 

6 

0  (all) 

’‘Hydrophone  station  numbers.  For 

ACODAC  users  there  should  be  entered 
as  they  are  ordered  on  the  analog 
tape,  i.e.,  the  first  number  should 
be  the  station  no.  on  channel  0, 
the  2nd  on  channel  1,  etc.  Unused 
channels  should  be  set  to  0. 

HSEN 

R 

6 

0.0  (al.) 

‘Hydrophone  sensitivity  numbers.  These 
are  ordered  in  the  same  manner  as  the 
station  numbers. 

GMT  I 

I 

1 

0 

Indicate  that  time  entered  on  shot  cards 
are  either  relative  tape  times  (0) 
or  GMT  (1).  If  GMT,  item  SGMT  must 
be  entered. 

TCMX 

I 

1 

10 

Maximum  allowable  number  of  unsuc¬ 
cessful  time  code  reads  allowed 
before  job  aborts. 

OVRX 

I 

1 

1 

Maximum  number  of  missed  digitizing 
bursts  allowed.  This  occurs  when 
burst  intervals  are  too  short  to 
allow  writing  the  data  to  mag  tape. 

TCWX 

A 

1 

62 

Maximum  number  of  seconds  (hydrophone 
time)  program  will  wait  for  a  time 
code  interrupt  before  assuming  read 
was  unsuccessful. 

TLAG 

I 

1 

55 

Time  lag  (in  seconds  of  hydrophone 
time)  between  actual  time  and  the 
finish  of  the  time  code  generation. 

TCDE 

I 

1 

1 

Time  code  indicator.  1  Time  code 
is  present  and  is  to  be  used  for  up¬ 
dating  internal  clock.  0  -*■  Time  code 
does  not  exist  and/or  is  not  to  be 
used. 
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Each  shot  specification  card  describes  a  digitizing  burst  or 
series  of  uniformly  spaced  bur-ts  in  terms  of  hydrophone  time. 
The  card  format  is  as  folio. (all  values  are  decimal  integer 
and  are  right-justified  in  the  field): 

Start  time  of  shot  or  shot  series 


column  1-3 
colunn  5-8 

colunn  10-11 


Days 

Hours  and  minutes  24  hour 
clock  format 
Seconds 


End  time  of  shot  series  (not  punched  for 
single  shot  card) 


column  13-15 
column  17-20 
column  22-23 


Days 

Hours  and  minutes 
Seconds 


3urst  information 


colunn  25-30 
column  32-35 


Burst  length  in  milliseconds 
Interval  between  start  of 
bursts  in  seconds  (not  punched 
for  single  shot  card) 


where  a  burst  is  defined  as  a  period  of  continuous  digitizing. 
Examples: 

To  specify  one  burst  of  1.5  seconds  duration  to  begin  at  time 
1425:35  of  day  1  (V  are  spaces): 

00 1 VI 425 V35  V  VV77V777W  WOO 1500 

To  specify  a  series  of  2.25  second  shots,  spaced  five  minutes 
apart  starting  at  2230:00  of  day  2  anu  ending  at  0130:00  of 
day  3: 

0027223070070037013070070022570300 

The  shot  cards  must  be  input  in  increasing  chronological  order 
and  there  must  be  no  intersection  of  time  intervals. 

II.  Analog  Input 

SA6PAD  reads  digitized  input  and  time  code  input  via  the  A/D  con¬ 
verter  attached  to  the  Sigma  7.  The  user  must  interface  his 
analog  tape  with  the  converter  and  provide  a  suitable  pulse  source 
to  determine  his  digitizing  rate.  He  should  consult  with  appro¬ 
priate  IPC  technicians  for  help  in  accomplishing  this  task. 

III.  Sigma  7  Console  Interrupt  Button 


This  is  used  to  signal  the  program  to  continue  operation  after  it 
has  entered  the  wait  state  to  allow  operator  intervention. 
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IV.  A/D  Converter  Manual  Input  Button 

This  is  used  when  operating  in  the  manual  mode.  It  signals  SAGPAD 
to  start  a  digitizing  burst. 

V.  Sense  Switches 

SAGPAD  interrogates  the  Sigma  7  console  sense  switches  periodically. 
Sense  switch  settings  are  interpreted  as  follows: 

Switch  #4 

If  switch  4  is  on,  the  program  will,  each  time  it  reads  a 
time  code,  log  out  its  value  and  the  internal  clock  setting 
on  the  line  printer.  Time  code  read  failures  will  also  be 
logged  out.  This  is  a  diagnostic  tool  and  should  not  be  used 
when  actually  digitizing. 

Switch  #3  and  #2 


Not  presently  used. 

Switch  Hi 

Setting  switch  1  on  will  force  the  program  to  terminate.  All 
normal  end  of  run  operations  are  executed.  This  is  the  normal 
way  to  terminate  a  manually  controlled  run.  It  may  also  be 
used  to  cut  short  an  automatically  controlled  run. 

OUTPUT:  1.  Line  Printer 

The  line  printer  is  used  for  listing  input  cards  and  logging  out 
error  and  diagnostic  messages.  If  an  error  is  specific  to  a 
particular  card,  the  error  message  is  logged  inmediately  following 
the  card  image. 

When  shot  cards  are  listed,  a  sequence  number  is  added.  If  an  error 
occurs  during  the  digitizing  process,  the  sequence  nunber  of  the 
currently  active  shot  specification  is  printed  along  with  the  error 
message.  This  gives  the  user  a  general  idea  of  where  the  program 
is  in  the  digitizer  "script. " 

I I .  Console  Teletype 

The  teletype  is  used  primarily  to  issue  instructions  to  the  opera¬ 
tor.  Generalized  error  messages  are  output  to  the  teletype  along 
with  recovery  instructions,  if  recovery  is  possible.  The  line 
printer  output  should  always  be  checked  if  an  error  occurs,  since 
it  will  generally  give  more  complete  information. 

III.  Magnetic  Tape 

The  digitized  data  is  output  to  9-track  800  BPI  magnetic  tape. 

Each  digitizing  run  results  in  one  data  file.  The  first  record 
contains  identification,  parameter  values  for  the  run  and  user 
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comments.  This  is  followed  by  data  records  containing  the  digitized 
data.  The  last  data  record  is  followed  by  a  double  end-of-file  mark. 

If  more  than  one  run  is  made  using  the  same  output  tape,  one  end- 
of-file  mark  separates  the  data  files  with  a  double  end-of-file 
following  the  last  file  on  the  tape. 

A  detailed  description  of  the  tape  format  follows.  In  the  header 
record  section,  the  column  labeled  "Input  Key"  identifies  the  in¬ 
put  parameter  which  produces  the  item  (see  INPUT,  Section  I-A). 

The  data  section  consists  of  digitized  points.  They  are  written 
just  as  received  from  the  digitizer. 

The  digitized  output  for  a  single  point  consists  of  16  bits  ar¬ 
ranged  as  follows: 


|S|  14  data  bits  |xj 

01  14  15 

If  single  point  per  transfer  is  specified  (NPPW=1),  the  16  bits 
will  be  right  justified  in  32  bit  Sigma  7  words  on  the  output 
tape.  If  double  point  transfer  is  specified  (NPPW=2),  two  16-bit 
points  will  be  packed  in  each  Sigma  7  word.  The  order  of  digitiz¬ 
ing  will  be  rightmost  point  first,  leftmost  second. 

The  interpretation  of  each  16-bit  point  is: 

Bits  0-14  represent  a  signed  binary  fraction  in  2*s  complement 
form.  The  actual  value  of  the  input  voltage  to  the  digitizer 
is  obtained  by  multiplying  this  fraction  by  the  gain  setting  being 
used  (set  by  user  via  parameter  GAIN).  In  other  words,  the  gain 
setting  determines  the  full  scale  representation,  plus  or  minus, 
of  the  digital  fraction. 

In  practice  the  full  scale  value  is  not  used,  since  a  maximum 
plus  or  minus  fraction  may  represent  an  overload  condition.  The 
ideal  situation  is  to  adjust  the  playback  equipment  so  that  the 
maximun  input  to  the  digitizer  is  slightly  less  than  one  of  the 
four  possible  gain  settings  (1.25,  2.5,  5.0,  or  10.0  volts),  thus 
obtaining  the  greatest  possible  resolution  without  causing  an 
overload  condition. 

Bit  15,  labeled  X,  is  used  as  an  end  of  scan  Indicator  when  sequen¬ 
tial  scanning  of  several  channels  is  specified.  If  this  bit  equals 
1,  it  implies  that  this  data  point  is  the  last  in  the  scan  and  that 
the  next  point  represents  the  start  of  a  new  scan. 

USAGE:  The  program  is  designed  as  a  stand  alone  system  and  must  be  loaded 
and  run  using  the  Xerox  Stand  Alone  Loader  with  I/O  Handlers.  The 
loader  is  supplied  as  a  bootable  absolute  binary  deck.  SAGPAD  is 
supplied  as  a  relocatable  binary  deck.  The  deck  set  up  is  shown 
below.  All  cards  explicitly  shown  are  essential.  (7  *  blank) 


I 

I 

I 

A 

I  Word 


I 


h 

ft 

K 


| 


ti- 


I 


£ 


o-i 


4 


5 

6 

7 

8 

9 

10 
11 


12 


13 

14 

15 

16 

17 

18 

19 

20 
21 
22 


23 


24 


25 
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ACODAC  -  SIGMA  7  DIGITIZE*  OUTPUT  TAPE 

Header  Record  Format 

Tape  I.D.  any  8  chars.  (EBCDIC) 

Max.  data  record  length  in  bytes  (not  incl.  header)  (Bin.  Integer) 

No.  pts.  per  word  (1  or  2)  (Bin.  Integer) 

Gain  setting  for  digitizer  (Bin.  Real) 

Type  of  scan,  1  =  random,  0  =  seq.  w/clock 
Control,  1  =  time  code  specs,  0  =  manual  control 

Digitizing  rate  per  second  per  channel  in  hydrophone  time  (Bin.  Int.) 
No.  of  channels  being  digitized  (Bin.  Int.  1-6) 

Playback  ratio  during  digitizing  (Bin.  Real) 

Deployment  No.  4  char.  EBCDIC 
Mooring  type  "  " 

Launch  day  no.  3  char.  EBCDIC  Ut.  Just. 

"  year  4  "  " 

Start  time  of  recording  GMT  Day  (Bin.  Int.) 


«» 

tl 

tl 

•  I 

tt 

Hour 

It 

It 

tt 

II 

tt 

ft 

Min. 

It 

Stop 

I' 

tl 

tt 

tt 

Day 

tt 

♦  1 

tt 

tt 

tt 

tt 

Hour 

tt 

It 

tt 

It 

tt 

tt 

Min. 

tt 

Duty  cycle  on  time  Min.  (Bin.  Int.) 

it  it  off  "  "  *• 

Recording  tape  speed  (Bin.  Real) 

Mooring  latitude  12  char.  EBCDIC  Rt.  Just. 

•  I  «•  tt  It  tl  ft  II 

•  I  M  If  tl  ft  It  II 


ACODAC 

Word 

26-28 

29-34 

35-40 

41 

49 

99 
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SIGMA  7  DIGITIZER  OUTPUT  TAPE  (Cont.) 

Header  Record  Format 

Mooring  longitude  12  EBCDIC  chars.  Rt.  Just.  (3  words) 

Hydrophone  sta.  no.'s  in  order  they  appear  in  data  (6  words)  (Bin.  Int.) 

"  sensitivity  no.'s  in  order  they  appear  in  data  (6  wds)  " 
Reserved  for  future  use  (9  words) 

Comments  (words  49-99) 


EOR  GAP 

Data  Record  Format 

Start  of  burst  time  in  elapsed  time  (secs.)  since  recording  start 
time  (Bin.  Int.) 

Record  number  within  burst  rt.  just.  (Bin.  Int.) 

Last  record  in  burst  indicator  (B3T  31)  0  =  not  last 

1  =  last 

DATA  -  number  of  bytes  given  by  word  #2  of  header  record. 


END  OF  RECORD  GAP 
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Stand  Alone  Loader  with  I/O  Handlers 

!A  F:LO,  (DEVICE,  LL) 

!A  F:OC,  (DEVICE, OC) 

2 A  F:C,  (DEVICE,  BI) 

!L  MAP,  M100,  N 

SAGPAD  relocatable  binary  deck 

2EOD 
2EOD 
!  R 

user's  input  deck  (see  INPUT,  section  I) 

The  loading  and  running  procedure  is  as  follows: 

A.  Place  deck  in  card  reader  and  ready  reader.  Ready  teletype 
and  line  printer  also.  Mount  your  output  tape,  dial  unit 
to  0  and  ready  unit. 

B.  Set  sense  switch  1  on. 

C.  Boot  from  card  reader. 

D.  When  card  reading  stops  (this  will  be  just  after  Loader  is 
read  in),  push  console  interrupt  button.  The  teletype  will 
print . 

2 ! KEY IN 

and  wait  fo--  input.  Type  the  following  message. 

2SYST  C.JRA03  NL 

This  tells  the  loader  to  accept  all  further  commands  from  the 
card  reader. 

E.  Set  sense  switch  1  off.  The  loader  will  load  SAGPAD  and  turn 
control  over  to  it. 

F.  SAGPAD  will  type  out 

HIT  CONSOLE  INTERRUPT  TO  START 

and  will  enter  the  wait  state.  This  occurs  prior  to  reading 
user's  input  decks.  If  this  is  a  rerun,  the  user  may  place 
his  next  input  deck  in  the  reader.  The  output  tape  should  be 
ready  on  unit  ABO  since  SAGPAD  will  write  a  header  record 
during  the  next  phase. 

When  the  console  interrupt  button  is  pushed,  SAGPAD  will  read  and 
process  the  input  deck.  If  an  error  is  detected,  SAGPAD  v/ill 
flag  it  and  enter  the  wait  state,  giving  the  user  an  opportunity 
to  correct  the  bad  card.  Hitting  the  console  interrupt  will 
cause  resunption  of  the  read.  See  "ERRORS  $  DIAGNOSTICS"  for 
more  details. 


I 
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G.  When  the  input  is  read  in  satisfactorily,  SAGPAD  will  write 
the  header  record  on  the  output  tape  and  type  out 

READY  TO  DIGITIZE  -  HIT  CONS.  INT.  TO  START 

and  will  enter  the  wait  state. 

The  user  should  now  ready  his  analog  input.  When  the  console 
interrupt  is  pressed,  the  program  will  initialize  its  internal 
clock  at  0  and  enter  the  digitizing  phase.  If  time  code  input 
is  indicated,  it  will  conmence  reading  the  time  code  channel 
and  will  reset  its  clock  to  stay  in  phase  with  the  time  code. 

If  the  program  is  being  run  in  automatic  mode  it  will  digitize 
as  called  for  in  the  script  set  up  by  the  shot  specification 
cards.  If  manual  mode,  it  will  digitize  only  one  command  from 
the  interrupt  button  on  A/D  unit. 

H.  Normal  program  termination  may  occur  in  one  of  two  ways. 

1.  If  in  automatic  mode,  the  program  will  terminate 
itself  when  all  the  digitizing  specified  on  the 
shot  cards  is  done.  When  this  occurs,  the  program 
will  type: 

DIGITIZING  RUN  COMPLETED. 

2.  The  operator  may  manually  terminate  the  run  by 
setting  sense  switch  1  on.  This  is  the  normal 
method  of  termination  for  a  manual  mode  operation 
but  it  may  also  be  used  to  effect  early  termination 
of  an  automatic  mode  run.  If  manually  terminated 
the  program  will  type: 

DIGITIZING  RUN  MANUALLY  TERMINATED. 

In  both  cases  the  program  will  write  two  end-of-file  marks 
on  the  output  tape,  backspace  over  the  second  end-of-file 
and  will  go  back  to  step  F,  allowing  the  user  to  make  another 
run. 

RESTRICTIONS: 

SAGPAD  may  only  be  run  with  the  Stand  Alone  Loader  with  I/O  Handlers. 

It  cannot  be  operated  under  the  Batch  Processing  Monitor. 

STORAGE  REQUIREMENTS 

4800j6  hex  locations  including  buffers  for  SAGPAD  proper.  The  Stand 

Alone  I/O  Handlers  must  also  be  in  core  at  run  time. 


/* 


t 


SUBPROGRAMS  REQUIRED;  None. 
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OPERATION  AL  ENVIRONMENT: 


Device 


Function 


console  interrupt 
console  teletype 
A/D  unit 
line  printer 
magnetic  tape 


input 

input/output 

input 

output 

output 


Special  Requirements 

none 

none 

See  INPUT  section 
9-track  800  BPI 
Unit  must  be  set  to 
"A80",  i.e.,  selector 
turned  to  0. 


OPERATIONAL  CHARACTERISTICS: 

Nonnal  for  Stand  Alone  Program 

TIMING:  There  are  two  areas  of  time  critical  operations.  The  first  is 
the  actual  digitizing  process,  the  second  is  the  transfer  of 
data  from  the  disk  buffer  to  magnetic  tape. 

A.  Digitizing 

The  maximum  allowable  digitizing  rate  (in  Sigma  7  or  wall 
clock  time)  is  100,000  data  points  per  second.  This  rate 
assumes  two  data  points  per  word  are  transferred.  The  only 
other  consideration  is  the  capacity  of  the  disk.  At  present 
S AGP AD  uses  one  7232  disk  unit  which  has  a  capacity  of  1.5 
million  Sigma  7  words  or  3.0  million  data  points  if  two  data 
points  per  word  are  used.  Therefore,  the  length  of  a  single 
continuous  burst  must  not  exceed  this. 

B.  Transfer  to  Magnetic  Tape 

The  critical  period  here  is  the  time  between  digitizing 
bursts.  It  must  be  long  enough  to  allow  data  collected  on 
the  disk  to  be  transferred  to  magnetic  tape.  This  is  dependent 
on  two  things.  The  number  of  data  points  in  the  preceding 
burst  and  the  magnetic  tape  record  size.  Two  things  must  be 
calculated.  First  the  time  to  read  the  data  from  the  disk, 
and  second,  the  time  to  write  the  data  to  tape. 

As  an  example,  assume  that  a  burst  contains  200,000  data 
points,  collected  as  two  points  per  32 -bit  Sigma  7  word. 

This  represents  100,000  Sigma  7  words  or  400,000  Sigma  7 
bytes.  Assume  also  that  a  magnetic  tape  record  size  of  2048 
bytes  is  specified.  This  size  will  also  be  used  for  disk 
reads.  For  each  disk  access,  the  maximum  latency  time  is 
.  034  seconds  and  the  transfer  rate  is  364,000  bytes  per 
second.  Therefore  the  maximum  time  required  to  read  the 
data  from  the  disk  is: 


•r  _  r400,000  v  . 

TD  ~204~8  X  '°54]  + 


400,000 
364,000  = 


7.74  seconds. 


To  calculate  the  time  required  to  write  the  data  to  magnetic 
tape  assune  an  inter-record  delay  time  of  .010  seconds  per 
record  and  a  transfer  rate  of  60,000  bytes  per  second. 
Therefore, 


_  ,400,000  , 

TT  “  1  2048  X  '0I0J 


400,000 

60,006 


8.63  seconds 


Therefore  total  transfer  time  (worst  case)  is  16.37  seconds, 
Sigma  7  time  and  the  user  must  allow  at  least  this  between 
bursts. 


Remember  that  these  are  wallclock  times  and  that  shot 
specifications  are  given  in  hydrophone  time.  If  the  play¬ 
back  rate  is  10  to  1,  this  means  that  the  specified  burst 
interval  must  be  greater  than  10  X  16.37  plus  burst  length  or 
163.7  seconds  plus  burst  length. 

It  is  evident  from  the  example  that  Ask  latency  or  seek  times 
and  inter-record  times  on  tape  account  for  much  of  the  total 
time.  Therefore,  one  can  cut  the  time  considerably  by  using 
records  which  are  as  long  as  practicable.  For  example,  if  the 
record  size  in  the  foregoing  illustration  were  doubled  from 
2048  bytes  to  4096  bytes,  the  time  would  be  decreased  from 
16.37  seconds  to  12.06  seconds. 

ERRORS  5  DIAGNOSTICS: 


Ihe  general  procedure  with  error  and  diagnostic  messages  is  to 
output  a  full  description  on  the  line  printer  and  to  either  repeat 
or  output  a  shorter  version  on  the  teletype.  In  addition,  the 
teletype  will  give  directions  for  recovery,  if  recovery'  is  possible. 

The  following  list  will  give  both  the  line  printer  (LP)  and  teletype 
(TTY)  message  for  each  error  or  diagnostic. 

A.  Parameter  Input  and  Initialization 

During  parameter  input  and  initialization,  the  following  messages 
may  appear. 

1.  (LP)  ERROR  IN  LAST  CARD 

(TTY)  PAR.  CARD  ERR.,  HIT  CONSOLE  INTERRUPT  TO  CONTINUE 

The  program  will  enter  the  wait  state.  The  error  may  be 
a  non-recognizable  key  or  an  error  in  the  parameter  value. 

To  recover,  correct  the  bad  card,  place  it  in  front  of  the 
unread  portion  of  the  deck  and  push  the  console  interrupt. 

2.  (LP)  MUX  GAIN  ILLEGAL 

(TTY)  PAR.  VAL.  ERR.  HIT  C.I.  TO  REREAD  PARAMETER  CARDS 

The  digitizer  gain  value  (key  -  GAIN)  must  be  one  of  four 
values,  1.25,  2.5,  5.0,  or  10.0.  This  message  occurs  after 
the  :SH0TS  card  is  read,  but  before  any  shot  cards  are 
processed.  To  recover,  correct  the  offending  parameter 
card,  retrieve  the  :SH0TS  card  from  the  output  hopper  and 
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place  the  corrected  parameter  card  and  the  :SHOTS  card 
(in  that  order)  in  front  of  the  unread  deck.  Hit  the 
console  interrupt  to  continue. 

3.  (LP)  TAPE  REC.  SIZE  MUST  BE  MULTIPLE  OF  1024  BYTES 

-  MAX.  =  16384. 

(TTY)  same  as  for  error  A-2 

Specific  tape  record  si2e  is  incorrect.  Recovery  is  the 
same  as  for  error  A-2 

4.  (TTY)  READY  TAPE  A80 

SAGPAD  is  ready  to  write  header  record,  but  tape  is  not 
ready.  Ready  tape  unit.  Message  will  repeat  every  five 
seconds  until  unit  is  readied. 

5.  (LP)  ERROR  LAST  SHOT  CARD 

(TTY)  SPEC.  ERR.,  CORRECT  AND  HIT  CON.  INT.  TO  REREAD 

There  is  an  error  in  the  last  shot  specification  card. 

Either  the  values  are  illegal,  there  is  a  time  inconsistency 
within  the  card,  or  there  is  a  time  overlap  with  the  last 
card.  To  recover,  correct  the  card,  place  it  at  the  head 
of  the  unread  input  deck  and  hit  the  console  interrupt  to 
continue  reading. 

6.  (LP)  TOO  MANY  SHOT  SPECS  -  LIMIT  *  100 

(TTY)  SHOT  OVERFLOW,  CORRECT  AND  HIT  CON.  INT.  TO  REREAD 

To  recover,  the  shot  deck  must  be  reduced  to  100  or  less 
and  the  whole  shot  deck,  including  the  :SHOT  card,  must 
be  placed  in  the  input  stack.  Hit  the  console  interrupt 
to  reread. 

B.  Digitizing 

The  following  error  - jssages  may  be  output  during  the  digitizing 
phase.  Each  of  the  messages  will  be  preceded  by  line  printer  output 
giving  the  sequence  nunber  of  the  shot  specification  card  currently 
being  processed  and  internal  (hydrophone)  time  information  consisting 
of  last  time  code  value  read,  if  any,  and  the  value  of  the  internal 
clock. 

1.  (LP  6  TIT)  RUN  ABORTED  -  TOO  MANY  TIME  OVERRUNS 

A  time  overrun  occurs  whenever  the  interval  between 
digitizing  bursts  is  too  short  to  allow  data  transfer  to 
magnetic  tape  (see  TIMING  section).  This  causes  SAGPAD 
to  skip  digitizing  the  next  schedule  burst.  The  maximum 
number  of  these  time  overruns  allowed  is  set  by  the  user 
via  parameter  0VRX.  When  this  number  is  exceeded  the  run 
is  aborted  and  SAGPAD  readies  itself  for  the  next  run. 

No  recovery  of  the  aborted  run  is  possible.  An  end-of- 
file  is  written  on  the  output  tape. 
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2.  CLP  S  TTY)  RUN  ABORTED  -  TOO  MANY  TIME  CODE  INTERRUPT 

FAILURES 

if  time  code  option  is  used  (TCDE=1),  SAGPAD  expects  a 
time  code  read  out  to  be  available  once  a  minute  (hydrophone 
time) .  Whenever  it  has  time  before  the  next  burst  is  due, 
it  attempts  to  read  the  time  code  channel.  If  an  interrupt 
is  not  received  within  TCWX  seconds  (hydrophone  time)  after 
issuing  a  read  command,  it  is  counted  as  a  time  code  failure. 
When  the  number  of  failures  exceeds  TCMX,  the  run  is  aborted 
in  the  same  manner  as  for  error  B-l. 

3.  (LP  §  TTY)  RUN  ABORTED  -  RAD  OVERFLOW 


During  actual  digitizing  SAGPAD  double  buffers  data  from 
core  to  disk.  If,  for  any  reason,  it  cannot  start  disk 
I/O  during  the  period,  it  aborts  with  the  above  message. 

The  cause  of  the  failure  may  either  be  due  to  hardware 
failure  in  the  Sigma  7  or  it  may  be  due  to  a  burst  length 
which  exceeds  the  capacity  of  the  disk  (see  TIMING  section). 
Abort  procedure  is  the  same  as  for  error  B-l. 


Failure  with  no  error  message 


Because  of  the  rapid  digitizing  rate  (up  to  100  KC) ,  SAGPAD  uses  a 
Read  Direct  loop  with  no  interrupt  when  reading  digitized  data. 

If  there  is  a  failure  in  the  A/D  *jnit  or  in  the  Sigma  7  direct  I/O 
logic,  the  program  will  hang  on  the  Read  Direct  instruction. 

There  is  no  way,  given  tne  timing  constraints,  that  the  program  can 
diagnose  this  problem.  In  order  to  aid  the  user  in  diagnosing  this 
situation,  the  location  of  the  Read  Direct  instruction  is  declared 
an  external  definition,  with  the  tag  READ.  Thus,  its  core  location 
will  be  contained  in  the  load  map  which  the  Stand  Alone  Loader 
outputs  on  the  line  printer  when  it  loads  SAGPAD.  The  user  should 
note  the  location.  If  the  program  hangs  at  this  core  address,  he 
can  then  be  certain  that  there  is  a  Read  Direct  failure  and  that 
the  problem  most  likely  lies  in  the  A/D  digitizer. 
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APPENDIX  E 

PRINCIPALS  OF  RMS  MEASUREMENTS  USING  THE  GENERAL  RADTO 
1926  ANALYZER  (RMS  DETECTOR) 


A  voltage  wave  form  v(t),  has  a  mean  value  v(t) 

v(t)  =  ^  /^vtdt  (1) 

The  mean  square  value  of  v(t),  taken  about  its  mean, 
is  called  the  variance,  a2 

a2  =  /g  [v(t)  -  v(t)]2dt  (2) 

The  square  root  of  a2  is  called  the  root -mean-square,  or  ims, 
value  of  v(t).  (Note  that  o2  excludes  dc.) 

If  we  assune  that  7(t)  =  0,  equation  (2)  reduces  to 

°2  =  f;o  [v(t)i2dt  (3) 

and  the  xms  value  becomes 


°  =  %  /J[v(t)]2dt  (4) 

When  the  data  are  discrete,  the  voltage  v(t)  is  known  only  at 
discrete  instants  of  time  (t).  The  integral  of  equation  (4) 
reduces  to  a  summation  and  we  have 


o 


/!  N 

if  l  [V(t)]z 

i=l 


(5) 


when  N  equals  the  number  of  samples. 

The  tms  value,  therefore  is  equal  to  the  square  root  of  the  sum 
of  the  squares  of  v(t^)  divided  by  the  nunber  of  samples. 

The  1926  Detector  measures  the  value  v(t.)  with  a  1-mv  reference 
level.  The  rms  level  in  db  above  1  mv  il  therefore 


xms  *  20  lOg10(^-)db 
o 

*  10  log10(~)2db 
o 


(6) 
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Combining  equations  5  and  6 


rms  =  10  log 


10 


i  9  '(*«)  , 

.(s.n-r^]2)ab 

i=l  o 


o„  and  are  reference  levels, 
o  o 

Equation  7  is  computed  by  the  1926  Detector  and  avoids  the 
square  root  operation  by  using  10  log.Q.  Also  N  is  aJways 
a  binary  number  so  that  division  by  N  is  simply  a  binary 
point  shift. 

Decibel  conversion 


All  output  measurements  of  the  1926  Analyzer  are  in  db  unit. 
During  the  output  phase  of  operation,  the  digital  memory  has 
the  mean  square  measurement  values  (called  the  variance,  a2) 
in  storage,  which  are  converted  to  rms  in  db,  using  the 
operations  described  by  equation  6: 

tins  =  10  iog10(^-)2db 
o 

The  memory  format  is  floating  point  binary,  so  that  a2  is  a 
number  stored  in  two  parts,  a  mantissa,  and  a  binary  exponent, 
E; 

o2  =  m.2^ 

Substituting  (8)  into  (6) 

rms  =  10  log1()m  +  10E  log10  2  db  -  10  log1Q  aQ 

=  10  log^m  ♦  3.0 IE  db  -  10  log1Q  oq 


The  (10  log1Qm)  term  of  equation  9  is  a  component  that  lies 
between  0  and  3  db,  because  m  is  a  number  that  lies  between 
1.0  and  2.0.  The  second  half  of  equation  9  is  a  component 
that  lies  between  0  and  60.2  db  because  E  is  an  integer  number 
that  lies  between  0  and  20.  The  third  term  is  a  constant  that 
accounts  for  the  reference  level. 

The  output  data  is  quantized  in  0.25  db  steps,  making  a  total 
of  240  possible  output  steps  for  the  dynamic  range  of  o0  db. 
Between  0  and  3  db  there  are  12  steps  (the  possible  outcomes 
of  the  output  computation  using  the  mantissa,  ra) . 

For  outputs  between  3  and  6  db,  the  exponent  E  equals  1  and 
again  12  steps  are  needed.  This  process  repeats  up  to  a 
maximun  of  60  db,  where  0  *  10  (actually  for  60.2  db). 


ttW'-.ssr  j^^saw***,  "«3«Ea^ 
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Using  12  steps  for  the  mantissa  and  20  for  the  exponent, 
therefore  a  total  of  240  steps  are  possible,  using  a  very 
simple  computational  algorithm. 

The  (3.01E)  term  of  .equation  9  is  not  a  whole  nunber,  but 
the  above  discussion  assumes  it  to  be.  In  the  actual 
computation  this  fractional  difference  is  ignored. 

Computation  of  the  first  half  of  equation  9  is  achieved  with 
a  fixed  conversion  matrix.  It  accepts  a  5-bit  input  (32 
possible  inputs)  and  has  13  possible  outputs  (12  steps  +  0) . 
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